NET与Matlab结合 —— 最小二乘法直线拟合(C#)

转载 2007年09月19日 16:35:00

首先是一个.m文件drawgraph.m,确保它能够在Matlab里运行。

我这里是最小二乘法直线拟合程序。

%最小二乘法直线拟合
%Created by Safirst C. Ke 2007.8.29 Wed 14:51

function drawgraph(coords)
%传入的参数为两行向量,第一行为x坐标,第二行为坐标。

%axis ([0 100 0 100]);
grid on;
hold on;

%显示欲拟合的点的位置
plot(coords(1,:), coords(2,:), '*');

%分解x,y坐标
x = coords(1,:)
y = coords(2,:)'

b = size(coords);
c = ones(1, b(2));
 
 MT = [c; x];
 M = MT';
 
 %f为直线函数,f = mx + b;
 f = inv(MT * M) * MT * y
 ['y = ', num2str(f(2)), 'x + ', num2str(f(1))]

 
 %显示最终拟合的直线
 x = -max(x):max(x);
 y =  f(1) + f(2) * x;
 
 
 plot(x, y);
 xlabel('X轴');
 ylabel('Y轴');
 title('最小二乘法直线拟合 by Safirst C. Ke');
 
 legend(['y = ', num2str(f(2)), 'x + ', num2str(f(1))]);

然后将这个文件包含在.NET的类库工程中,并进行编译。

这里需要理解它的过程,毕竟.NET不能编译.m文件。怎么做到的呢?

通过设置这个工程的生成事件属性,添加为

call PlotDemoBuild.bat

然后在PlotDemoBuild.bat这个文件里面写好用Matlab编译器mcc编译的命令行,最重要的部分就是

mcc -M -silentsetup -vg -B "dotnet:PlotDemoComp,Plotter,2.0,private" -d ../../src ../../drawgraph.m

这样的话,点击生成,就会通过mcc产生dll,即我们需要的类库。

然后建立我们真正的C#工程,添加引用为刚才的类库,并开始写程序program.cs

using System;
using System.Collections.Generic;
using System.Text;

using MathWorks.MATLAB.NET.Utility;    

using MathWorks.MATLAB.NET.Arrays;  

//这两个引用显然要添加,不过好在这两个命名空间属于一个库MWArray.dll

//C:/Program Files/MATLAB/R2007a/toolbox/dotnetbuilder/bin/win32/v2.0/MWArray.dll

using PlotDemoComp;

namespace ConsoleApplication2
{
    class Program
    {
        //[STAThread]
        static void Main(string[] args)
        {
            try
            {
                Console.WriteLine("Please Input the points you want to fit:");

                string[] y = Console.ReadLine().Trim().Split();

                int size = y.Length;
                double[] x = new double[size];
                for(int i = 0; i < size; i++)
                {
                    x[i] = Convert.ToDouble(y[i]);
                }

                double[,] pointValues = new double[2, size / 2];

               //从开头算起,相邻的两个数为一个点,所以x和y都是间隔一个的。如1,2,3,4代表两点(1,2),(3,4)

                for(int i = 0; i < size; i += 2)   

                {
                    int index = i / 2;
                    pointValues[0, index] = x[i];
                }
                for(int i = 1; i < size; i += 2)
                {
                    int index = (i - 1) / 2;
                    pointValues[1, index] = x[i];
                }
                Plotter plotter = new Plotter();
                plotter.drawgraph((MWNumericArray)pointValues);
                Console.ReadLine();
            }
            catch(Exception exception)
            {
                Console.WriteLine("Error: {0}", exception);
            }
        }
    }
}

运行结果如下:

Please Input the points you want to fit:
1 2 3 4 5 6 -1 -2 -3 -4 -5 -6

*号标记欲拟合的点,直线为拟合直线!

就写这么多,以后再加入一个曲线拟合的程序。

特别声明:在下学习.NET C#和Matlab的时间均不超过1个月,

望高手不要批评指责,照顾一下新手学习钻研的积极性! 

直线拟合之最小二乘法 Ax+By+C=0

最小二乘法拟合直线,对于大多数人都不是很陌生,直线方程 y=kx+b,令sum((y-kx-b)^2),分别对k和b求导便可求出最佳的参数,但是如果遇到诸如“X=5”这样的直线方程,又该如何呢?这时候...
  • linxue968
  • linxue968
  • 2014年02月25日 14:49
  • 3393

最小二乘法(c语言实现线性,matlab进行拟合)及相关系数的求解

现在给定n个点,(x1,y1),(x2,y2),(x3,y3),(x4,y4),(x5,y5)..(xn,yn).现在希望得到一条最好的曲线(也就是求一个函数关系式~~)能尽可能的描述这n个点(不一定...
  • fengsigaoju
  • fengsigaoju
  • 2016年03月31日 13:01
  • 4662

关于Matlab中的线性与非线性最小二乘拟合

1、线性最小二乘拟合 最小二乘法(又称最小平方法)是一种数学优化技术,其通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据,并使得这些求得的数据与实际数据之间误差的...
  • u012117071
  • u012117071
  • 2015年08月24日 11:17
  • 3508

最小二乘法直线拟合及其Matlab实现

最小二乘法,通常用在我们已知数学模型,但是不知道模型参数的情况下,通过实测数据,计算数学模型,例如,在题目中,数学模型就是直线方程y=ax+b,但是不知道直线方程的a和b。     本来呢,我们只需...
  • cp32212116
  • cp32212116
  • 2014年12月10日 16:47
  • 10017

最小二乘法直线拟合程序(VB6.0代码编写)

  • 2009年04月24日 10:48
  • 2.53MB
  • 下载

最小二乘法拟合直线公式的初等推导

  • 2015年11月28日 18:53
  • 133KB
  • 下载

最小二乘法直线拟合及C++实现

  • 2010年08月04日 21:29
  • 21KB
  • 下载

最小二乘法拟合直线程序

  • 2010年12月19日 09:18
  • 4KB
  • 下载

【算法研究与实现】最小二乘法直线拟合

1.原理  在现实中经常遇到这样的问题,一个函数并不是以某个数学表达式的形式给出,而是以一些自变量与因变量的对应表给出,老师讲课的时候举的个例子是犯罪人的身高和留下的脚印长,可以测出一些人的数据然后得...
  • u012374012
  • u012374012
  • 2015年03月18日 14:36
  • 700

最小二乘法拟合直线c++代码

最近公司的一个项目需要计算TVDI(Temperature Vegetation Dryness Index ,温度植被干旱指数) ,TVDI的计算公式如下(具体原理自行百度): ...
  • liuxiangxxl
  • liuxiangxxl
  • 2017年12月18日 20:41
  • 48
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:NET与Matlab结合 —— 最小二乘法直线拟合(C#)
举报原因:
原因补充:

(最多只允许输入30个字)