作业帮 > 数学 > 作业

请帮忙编写一个matlab计算程序,拟合二元一次函数z=f(x,y)

来源:学生作业帮 编辑:作业帮 分类:数学作业 时间:2024/11/07 23:58:09
请帮忙编写一个matlab计算程序,拟合二元一次函数z=f(x,y)
希望得到一个z=ax+by+c的二元一次函数
其中:
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
请帮忙别写一个拟合程序
目前这个运算结果最后生成的图形是一个平面。如果我想将原始数据xyz描出一个曲面,能不能再帮忙补充两句代码,同时将给出的已知数据描点连成曲面。
clear
clc

x=[50 100 150 200 250 300 350 400 450  500];
y=[40 80 120 160 200];
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];
[y,x]=meshgrid(y,x);
p=[x(:),y(:),ones(numel(z),1)]\z(:)%p=[a;b;c]
plot3(x(:),y(:),z(:),'ro');
hold on;grid on
mesh(x,y,p(1)*x+p(2)*y+p(3))