rwhn.net
当前位置:首页 >> 请帮忙编写一个mAtlAB计算程序,拟合二元一次函数z... >>

请帮忙编写一个mAtlAB计算程序,拟合二元一次函数z...

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....

楼主说的这个问题是线性拟合问题,不必用专门的拟合命令: 推导: z = a0+a1*x+a2*y ==> z = [1 x y]*[a0 a1 a2] 所以,要求未知参数列向量a,直接用矩阵除法就可以了: a = xy\z 其中 xy = [1 1 0 1 2 0 1 1 1 1 2 1 1 0 1 1 0 2] %第2列为x,...

[x,y]=meshgrid(-9:0.5:9);z = 2*sin(x).*sin(y)./(x.*y); % z的表达式sinyy是什么?mesh(x,y,z);

>> x=[3,4,5,6,7,8,9,10]; y=[4,5,6,7,8,9,10,11]; z=21.89-9.87*x+8.09*y-7.98*x.*y-1.87*x.^2-7.98*y.^2 z = 1.0e+003 * -0.2156 -0.3662 -0.5524 -0.7742 -1.0317 -1.3249 -1.6537 -2.0182 >>

x=linspace(-5,5,50); y=linspace(-5,5,50); [x,y]=meshgrid(x,y); z=sin(x).*sin(y); figure(1) meshc(x,y,z); figure(2) [C, h] = contourf(z); clabel(C, h, 'LabelSpacing', 72); colorbar; tm1 = colormap(hot(9)); tm1(1, :) = [1 1 1]; c...

>> z1=@(x1,y1) exp(-1./(x1.^2+y1.^2)); >> z1(1,3) ans = 0.9048 >> ezsurf(z1) 第二个尝试了下: >> z2=@(x2,y2) dblquad(@(x1,y1) z1(x1,y1).*exp(-2*pi.*(x1.*x2+y1.*y2)),-100,100,-100,100); >>z2(1,2) 貌似算不出来,不知怎么弄啊

举个例子:已知z=a*(exp(y)+1)-sin(x)*b且有多组(x,y,z)的值,现在求最佳系数a,b >> f=@(a,x)a(1)*(exp(x(2,:))+1)-sin(x(1,:))*a(2);>> x=2:10;>> y=10*sin(x)./log(x);>> z=4.5*(exp(y)+1)-sin(x)*13.8;>> lsqcurvefit(f,[1 2],[x;y],z) %注意[...

可以转化成三重积分:int(int(int(f,z,z1(x,y),z2(x,y)),y,y1(x),y2(x)),x,a,b)

你看看程序,不懂得可以问我。 %程序 xgrid = -2:0.1:2; ygrid = -2:0.1:2; [x,y] = meshgrid(xgrid,ygrid); z_p = sqrt(x.^2+y.^2); z_n = -sqrt(x.^2+y.^2); mesh(x,y,z_p);%锥面上半部分 hold on; mesh(x,y,z_n);%锥面下半部分 hold on; t=0:...

参考代码: % 生成测试数据X = rand(100,1);Y = rand(100,1);k = 2;a = 1.5;b = 2.5;Z = k*(X.^a).*(Y.^b); % 数据加入噪声Z = Z + randn(size(Z))*0.1; % 数据拟合f=inline('c(1)*x(:,1).^c(2).*x(:,2).^c(3)','c','x');c=lsqcurvefit(f,[1 1 1...

网站首页 | 网站地图
All rights reserved Powered by www.rwhn.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com