请matlab高手根据给出的数据编写一个拟合二元函数的程序(本人完全不懂matlab)希望得到一个二元不限次的函数,一次的肯定不行,从图形上看偏差很大.建议先将数据点描出,再根据分布趋势确
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/05 18:49:54
![请matlab高手根据给出的数据编写一个拟合二元函数的程序(本人完全不懂matlab)希望得到一个二元不限次的函数,一次的肯定不行,从图形上看偏差很大.建议先将数据点描出,再根据分布趋势确](/uploads/image/z/6843704-32-4.jpg?t=%E8%AF%B7matlab%E9%AB%98%E6%89%8B%E6%A0%B9%E6%8D%AE%E7%BB%99%E5%87%BA%E7%9A%84%E6%95%B0%E6%8D%AE%E7%BC%96%E5%86%99%E4%B8%80%E4%B8%AA%E6%8B%9F%E5%90%88%E4%BA%8C%E5%85%83%E5%87%BD%E6%95%B0%E7%9A%84%E7%A8%8B%E5%BA%8F%EF%BC%88%E6%9C%AC%E4%BA%BA%E5%AE%8C%E5%85%A8%E4%B8%8D%E6%87%82matlab%EF%BC%89%E5%B8%8C%E6%9C%9B%E5%BE%97%E5%88%B0%E4%B8%80%E4%B8%AA%E4%BA%8C%E5%85%83%E4%B8%8D%E9%99%90%E6%AC%A1%E7%9A%84%E5%87%BD%E6%95%B0%2C%E4%B8%80%E6%AC%A1%E7%9A%84%E8%82%AF%E5%AE%9A%E4%B8%8D%E8%A1%8C%2C%E4%BB%8E%E5%9B%BE%E5%BD%A2%E4%B8%8A%E7%9C%8B%E5%81%8F%E5%B7%AE%E5%BE%88%E5%A4%A7.%E5%BB%BA%E8%AE%AE%E5%85%88%E5%B0%86%E6%95%B0%E6%8D%AE%E7%82%B9%E6%8F%8F%E5%87%BA%2C%E5%86%8D%E6%A0%B9%E6%8D%AE%E5%88%86%E5%B8%83%E8%B6%8B%E5%8A%BF%E7%A1%AE)
请matlab高手根据给出的数据编写一个拟合二元函数的程序(本人完全不懂matlab)希望得到一个二元不限次的函数,一次的肯定不行,从图形上看偏差很大.建议先将数据点描出,再根据分布趋势确
请matlab高手根据给出的数据编写一个拟合二元函数的程序
(本人完全不懂matlab)
希望得到一个二元不限次的函数,一次的肯定不行,从图形上看偏差很大.建议先将数据点描出,再根据分布趋势确定函数类型,这需要根据经验了.最后希望得到完整的代码,运行之后可以得到拟合的方程、原始数据点图和拟合出的函数曲面图.
其中:
x:50 100 150 200 250 300 350 400 450 500.
y:40 80 120 160 200
z是一个10X5的矩阵(其实就是X行和Y列得到的数据点):
0.05 0.05 0.05 0.05 0.05
0.25 0.15 0.11 0.05 0.05
0.5 0.2 0.2 0.04 0.04
1.6 1 0.5 0.11 0.07
2.5 2.4 1.14 0.34 0.12
3.2 2.7 1.54 0.7 0.21
3.4 3.1 2.1 1.15 0.3
3.75 3.5 2.56 1.4 0.6
4 3.8 2.95 2 0.9
4.2 4 3.2 2.2 1.1
请matlab高手根据给出的数据编写一个拟合二元函数的程序(本人完全不懂matlab)希望得到一个二元不限次的函数,一次的肯定不行,从图形上看偏差很大.建议先将数据点描出,再根据分布趋势确
在不知道原理的情况下,函数的选择本身就有些随意性.我选了一种函数来拟合,程序和结果你运行后就能看见.
clear all
clc
x=[50 100 150 200 250 300 350 400 450 500];
y=[40 80 120 160 200];
[X,Y]=meshgrid(x,y);
z=[0.05 0.05 0.05 0.05 0.05
0.25 0.15 0.11 0.05 0.05
0.5 0.2 0.2 0.04 0.04
1.6 1 0.5 0.11 0.07
2.5 2.4 1.14 0.34 0.12
3.2 2.7 1.54 0.7 0.21
3.4 3.1 2.1 1.15 0.3
3.75 3.5 2.56 1.4 0.6
4 3.8 2.95 2 0.9
4.2 4 3.2 2.2 1.1]';
DF=@(x) sum(sum((z-(x(1)*X.^2./(x(2)*X.^2+x(3)*Y.^3+x(4)))).^2));
[coef,fval]=fminsearch(DF,[1,1,1,1]);
[X1,Y1]=meshgrid(0:20:500,0:20:200);
Fs=coef(1)*X1.^2./(coef(2)*X1.^2+coef(3)*Y1.^3+coef(4));
surf(X,Y,z,'FaceColor','none','Marker','.','MarkerSize',30,'MarkerFaceColor',[0 0 0])
hold on
surf(X1,Y1,Fs,'FaceAlpha',0.5,'FaceColor','interp')
xlabel('x')
ylabel('y')
zlabel('z')
syms x y a b c d
Fitfun=a*x^2/(b*x^2+c*y^3+d);
disp('The fitting function is ')
pretty(Fitfun)
disp(['where a=',num2str(coef(1)),', b=',num2str(coef(2)),', c=',num2str(coef(3)),', d=',num2str(coef(4))])
disp(['The sum of least-square errors :',num2str(fval)])