{"id":5298,"date":"2014-01-08T09:31:10","date_gmt":"2014-01-08T00:31:10","guid":{"rendered":"http:\/\/www.moonmile.net\/blog\/?p=5298"},"modified":"2014-01-06T23:32:12","modified_gmt":"2014-01-06T14:32:12","slug":"f-%e9%80%a3%e7%ab%8b%e4%b8%80%e6%ac%a1%e6%96%b9%e7%a8%8b%e5%bc%8f%e3%82%92%e8%a7%a3%e3%81%8f","status":"publish","type":"post","link":"http:\/\/www.moonmile.net\/blog\/archives\/5298","title":{"rendered":"[F#] \u9023\u7acb\u4e00\u6b21\u65b9\u7a0b\u5f0f\u3092\u89e3\u304f"},"content":{"rendered":"<p>\n\u9006\u884c\u5217\u3001\u884c\u5217\u5f0f\u3068\u3067\u304d\u305f\u306e\u3067\u3001\u9023\u7acb\u4e00\u6b21\u65b9\u7a0b\u5f0f\u3092\u89e3\u304b\u305b\u308b\u3002\u6709\u9650\u8981\u7d20\u6cd5\u306e\u3044\u308f\u3086\u308b\u300csolve\u300d\u3068\u3044\u3046\u3084\u3064\u3067\u3001\u89e3\u6790\u306e\u624b\u9806\u3068\u3057\u3066\u3001\n<\/p>\n<ul>\n<li>pre : \u69cb\u9020\u3092\u8a2d\u5b9a\u3001\u8981\u7d20\u306b\u5206\u5272<\/li>\n<li>solve : \u9023\u7acb\u4e00\u6b21\u65b9\u7a0b\u5f0f\u3092\u89e3\u304f\uff08\u305d\u308c\u3060\u3051\u3058\u3083\u306a\u3044\uff1f\uff09<\/li>\n<li>post : \u3072\u305a\u307f\u306e\u7d50\u679c\u3092\u8868\u793a\u3059\u308b<\/li>\n<\/ul>\n<p>\n\u3068\u3044\u3046\u624b\u9806\u306b\u306a\u3063\u3066\u3044\u308b\u3002<br \/>\nGauss \u306e\u6d88\u53bb\u6cd5\u3092\u4f7f\u3063\u3066\u76f4\u63a5\u7684\u306b\u89e3\u304f\u308f\u3051\u3060\u304c\u3001\u3053\u308c\u3063\u3066\u3001\u305d\u306e\u307e\u3093\u307e LU \u5206\u89e3\u306a\u308f\u3051\u3067\u3002\u306a\u306e\u3067\u3001solve \u81ea\u4f53\u306f\u53cd\u5fa9\u6cd5\u306a\u3069\u3082\u542b\u3081\u305f\u610f\u5473\u3067\u4f7f\u308f\u306a\u3044\u3068\u30c0\u30e1\u3063\u307d\u3044\u3002\n<\/p>\n<pre class=\"brush: csharp; title: ; notranslate\" title=\"\">\nlet solve (A&#039; : matrix) (F&#039;: vector) = \n    let A = A&#039;.Copy()\n    let F = F&#039;.Copy()\n    let n = A.NumCols-1\n    \/\/ \u524d\u9032\u6d88\u53bb\n    for i in 0..n-1 do\n        for j in i+1..n do\n            let p = A.&#x5B;j,i]\/A.&#x5B;i,i]\n            F.&#x5B;j] &lt;- F.&amp;#91;j&amp;#93; - p * F.&amp;#91;i&amp;#93;\n            for k in 0..n do\n                A.&amp;#91;j,k&amp;#93; &lt;- A.&amp;#91;j,k&amp;#93; - p * A.&amp;#91;i,k&amp;#93;\n    \/\/ \u5f8c\u9000\u4ee3\u5165\n    F.&amp;#91;n&amp;#93; &lt;- F.&amp;#91;n&amp;#93;\/A.&amp;#91;n,n&amp;#93;\n    for i in n-1..-1..0 do\n        for j in i+1..n do\n            F.&amp;#91;i&amp;#93; &lt;- F.&amp;#91;i&amp;#93; - A.&amp;#91;i,j&amp;#93;*F.&amp;#91;j&amp;#93;\n        F.&amp;#91;i&amp;#93; &lt;- F.&amp;#91;i&amp;#93;\/A.&amp;#91;i,i&amp;#93;\n    F\n&amp;#91;\/code&amp;#93;\n&lt;p&gt;\n\u300c\u6709\u9650\u8981\u7d20\u6cd5\u6982\u8aac\u300d\u306b\u3042\u3063\u305f\u524d\u9032\u6d88\u53bb\u300c\u4e0a\u4e09\u89d2\u5316\u300d\u3057\u305f\u3042\u3068\u306b\u5f8c\u9000\u4ee3\u5165\u3067\u3072\u3068\u3064\u305a\u3064\u8a08\u7b97\u3059\u308b\u3002\n&lt;\/p&gt;\n&#x5B;code lang=&quot;csharp&quot;]\n\/\/ 2x -y = 2\n\/\/ -x +2y= 5 \u3092\u89e3\u304f\n\nlet M = matrix&#x5B;&#x5B;2.;-1.];\n               &#x5B;-1.;2.]]\nlet F = vector&#x5B;2.0;5.0]\n\nsolve M F \n<\/pre>\n<p>\n\u3053\u3093\u306a\u98a8\u306b\u3044\u308c\u3066\u304a\u3044\u3066\u3001\u5b9f\u884c\u3059\u308b\u3068\n<\/p>\n<pre class=\"brush: csharp; title: ; notranslate\" title=\"\">\nval it : Vector&lt;float&gt; = vector &#x5B;|3.0; 4.0|]\n<\/pre>\n<p>\n\u3053\u3093\u306a\u98a8\u306b\u7d50\u679c\u304c\u51fa\u308b\u3002\n<\/p>\n<p>\n\u5b9f\u306f\u3001F# MathProvider \u306b\u3082\u540c\u3058\u3082\u306e\u304c\u3042\u3063\u3066\u3001\n<\/p>\n<pre class=\"brush: csharp; title: ; notranslate\" title=\"\">\nlet M = matrix&#x5B;&#x5B;2.;-1.];\n               &#x5B;-1.;2.]]\nlet F = vector&#x5B;2.0;5.0]\n\nlet sol = MathProvider.LinearAlgebra.solve M F\n<\/pre>\n<p>\n\u3068\u3057\u3066\u89e3\u304f\u3068\u540c\u3058\u4e8b\u304c\u3067\u304d\u308b\u3002\n<\/p>\n<p>\n\u6709\u9650\u8981\u7d20\u6cd5\u306e\u5834\u5408\u3001\u3072\u3068\u3064\u306e\u7bc0\u70b9\u306e\u81ea\u7531\u5ea6\u306f6\u306b\u306a\u308b\u306e\u3067\u3001\u7bc0\u70b9\u306e\u6570\u306e6\u500d\u306e\u8a08\u7b97\u91cf\u306b\u306a\u308b\u3002\u3055\u3089\u306b\u3001\u8981\u7d20\u304c\u6301\u3064\u7bc0\u70b9\uff08\u9802\u70b9\uff09\u306f\u4e09\u89d2\u5f62\u4e8c\u6b21\u8981\u7d20\u5834\u5408\u306f6\u70b9\u306b\u306a\u308b\u306e\u3067\u3001\u8981\u7d20\u6570\u306e36\u500d\u306e\u8a08\u7b97\u91cf\u306a\u308b\u3002\u304c\u3001\u8a08\u7b97\u91cf\u304c\u591a\u304f\u306a\u308b\u3068\u3044\u3063\u3066\u3082\u3001\u591a\u304f\u306a\u3063\u3066\u30822\u6841\u5897\u3048\u308b\u3050\u3089\u3044\u3060\u304b\u3089\u5927\u3057\u305f\u3053\u3068\u306f\u306a\u3044\uff1f\u3053\u3068\u306f\u306a\u3044\u304b\u3002\u524d\u9032\u6d88\u53bb\u3067\u3001n^3\/2 \u3050\u3089\u3044\u306e\u30aa\u30fc\u30c0\u30fc\u306b\u306a\u308b\u306e\u3067\u3001n \u304c 100\u500d\u306b\u306a\u308b\u3068\u3001100\u4e07\u500d\u3050\u3089\u3044\u306e\u8a08\u7b97\u91cf\u304c\u5897\u3048\u308b\u3002\u3068\u306a\u308b\u3068\u3001\u3053\u306e\u7d20\u76f4\u306a\u65b9\u6cd5\u3067\u89e3\u3044\u3066\u3057\u307e\u3046\u3068\u3001\u30e1\u30c3\u30b7\u30e5\u6570\u30922\u500d\u7d30\u304b\u304f\u3059\u308b\u3068\u7d048\u500d\u306e\u8a08\u7b97\u91cf\u304c\u5fc5\u8981\u306b\u306a\u308b\u3002\u78ba\u304b\u306b\u3001\u9ad8\u901f\u5316\u304c\u5fc5\u8981\u306a\u5206\u91ce\u3067\u3042\u308b\u304b\u3082\u3002<br \/>\n\u3053\u306e\u3042\u305f\u308a\u306f\u5f8c\u3067\u5b9f\u6e2c\u3057\u3066\u8a66\u3057\u3066\u307f\u3088\u3046\u3002<\/p>\n","protected":false},"excerpt":{"rendered":"<p>\u9006\u884c\u5217\u3001\u884c\u5217\u5f0f\u3068\u3067\u304d\u305f\u306e\u3067\u3001\u9023\u7acb\u4e00\u6b21\u65b9\u7a0b\u5f0f\u3092\u89e3\u304b\u305b\u308b\u3002\u6709\u9650\u8981\u7d20\u6cd5\u306e\u3044\u308f\u3086\u308b\u300csolve\u300d\u3068\u3044\u3046\u3084\u3064\u3067\u3001\u89e3\u6790\u306e\u624b\u9806\u3068\u3057\u3066\u3001 pre : \u69cb\u9020\u3092\u8a2d\u5b9a\u3001\u8981\u7d20\u306b\u5206\u5272 solve : \u9023\u7acb\u4e00\u6b21\u65b9\u7a0b\u5f0f\u3092\u89e3\u304f\uff08\u305d\u308c\u3060\u3051\u3058\u3083\u306a\u3044\uff1f\uff09 p &hellip; <a href=\"http:\/\/www.moonmile.net\/blog\/archives\/5298\">\u7d9a\u304d\u3092\u8aad\u3080 <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[59],"tags":[],"class_list":["post-5298","post","type-post","status-publish","format-standard","hentry","category-f"],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"http:\/\/www.moonmile.net\/blog\/wp-json\/wp\/v2\/posts\/5298","targetHints":{"allow":["GET"]}}],"collection":[{"href":"http:\/\/www.moonmile.net\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"http:\/\/www.moonmile.net\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"http:\/\/www.moonmile.net\/blog\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"http:\/\/www.moonmile.net\/blog\/wp-json\/wp\/v2\/comments?post=5298"}],"version-history":[{"count":1,"href":"http:\/\/www.moonmile.net\/blog\/wp-json\/wp\/v2\/posts\/5298\/revisions"}],"predecessor-version":[{"id":5299,"href":"http:\/\/www.moonmile.net\/blog\/wp-json\/wp\/v2\/posts\/5298\/revisions\/5299"}],"wp:attachment":[{"href":"http:\/\/www.moonmile.net\/blog\/wp-json\/wp\/v2\/media?parent=5298"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"http:\/\/www.moonmile.net\/blog\/wp-json\/wp\/v2\/categories?post=5298"},{"taxonomy":"post_tag","embeddable":true,"href":"http:\/\/www.moonmile.net\/blog\/wp-json\/wp\/v2\/tags?post=5298"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}