在matlab中使用矩阵拟合2维函数
来源:互联网 发布:美洽电脑软件下载 编辑:程序博客网 时间:2024/04/29 22:03
最近要解决用2维矩阵拟合2维函数的问题,貌似在中文网站中没有找到方法,最后在stack overflow才问到了答案,解答如下。
假设有5*5数组:
d= [0.0177104427823448,0.00246661459209512,0.0399831543374395,0.0615494164555707,0.0476204124707652;0.0275276152854314,0.0219153841813084,0.0581144391404502,0.144890028400954,0.157839631316098;0.0622883972729130,0.0716157303159909,0.245482781674067,0.123999612575059,0.177495187746408;0.0200735764542146,0.0573087934038160,0.0636451189717613,0.0160810084568415,0.0484992279558924;0.0185180386159227,0.00841167700273800,0.0372017422726281,0.0173721095082637,0.0459520362441099]
想用最小均方法拟合2维函数:
r = alfa*sin(pi*(n1+delta1))*sin(pi*(n2+delta2)) / (25*sin(pi/5*(n1+delta1))*sin(pi/5*(n2+delta2)))
[n,m]=size(d);%assumes that d is a n x m matrix[X,Y]=meshgrid(1:n,1:m);%your x-y coordinatesx(:,1)=X(:); % x= first columnx(:,2)=Y(:); % y= second columnf=d(:); % your data f(x,y) (in column vector)%--- now define the function in terms of x%--- where you use x(:,1)=X and x(:,2)=Yfun = @(c,x) c(1)*sin(pi*(x(:,1)+c(2))).*sin(pi*(x(:,2)+c(3))) ./ (25*sin(pi/5*(x(:,1)+c(2))).*sin(pi/5*(x(:,2)+c(3))));%--- now solve with lsqcurvefitoptions=optimset('TolX',1e-6);c0=[1 0 0];%start-guess herecc=lsqcurvefit(fun,c0,x,f,[],[],options);Ifit=fun(cc,x); Ifit=reshape(Ifit,[n m]);%fitting data reshaped as matrixsurf(X,Y,Ifit);hold on;plot3(X, Y, dataArray);
- 在matlab中使用矩阵拟合2维函数
- 数据拟合---使用自定义函数进行非线性拟合 -在Origin。matlab拟合工具箱cftool
- 数据拟合---使用自定义函数进行非线性拟合 -在Origin。matlab拟合工具箱cftool
- 数据拟合---使用自定义函数进行非线性拟合 -在Origin。matlab拟合工具箱cftool
- matlab函数公式拟合
- Matlab 自定义函数拟合
- MATLAB拟合函数使用说明
- matlab 给定函数 拟合
- MATLAB 自定义函数拟合
- matlab 函数拟合polyfit
- MATLAB拟合函数使用说明
- matlab中矩阵运算函数
- 在C#应用中使用Matlab函数
- 关于matlab中数据拟合的一些函数polyfit()
- Matlab三维拟合函数griddata
- Matlab 使用squeeze函数提取多维矩阵的某一维
- MATLAB中数据数据拟合
- Matlab拟合工具箱的使用
- do while(FLASE)的妙用【转】
- hdu 4602题解
- 针对乐视网的主页无法打开的解决办法
- HDU--杭电--2191--悼念512汶川大地震遇难同胞——珍惜现在,感恩生活--背包
- uva839 Not so Mobile (重建二叉树)
- 在matlab中使用矩阵拟合2维函数
- ListView item设置间距无效问题
- android基础篇------------java基础(4) (面向对象程序设计)
- 认识js中的function和this
- 飞跃的心
- poj3667 线段树+区间合并
- ACM 动态规划初步 Partitioning by Palindromes
- hdu 2013 蟠桃记
- Aho-Corasick算法