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

>> 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) 貌似算不出来,不知怎么弄啊

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

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

参考代码: % 生成测试数据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...

你看看程序,不懂得可以问我。 %程序 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=[1 2 3 4 5 6]'; y=[7 8 9 6 5 4]'; z=[12 32 25 34 54 21]'; X=[ones(size(x)) x x.^2 y y.^2 x.*y]; Y=z; [B,BINT,R,RINT,STATS] = REGRESS(Y,X) %B(1)-B(6) 相当于a,b,c,d,e,f的值 rcoplot(R, RINT) syms x y F=B'*[1 x x.^2 y y....

function z = f(x, y) if y == 1 z = sin(x); else if y == 2 z = cos(x); else z = sin(x)*cos(x); endend

曲面拟合可以按下列步骤进行: 1、根据x,y,z数据,可以用cftool拟合工具箱的常用函数去拟合,判断(x,y,z)大概符合那个曲面方程。(注:cftool拟合工具箱有时候不一定与实际相符) 2、建立曲面方程模型函数,如func=@(a,x) a(1)*x(1)^2+a(2)*...

再已知x y z的情况下, 比如 你已经把代码列出 X(:,1)=x; X(:,2)=y; Y(:,1)=z; 接下来可以用非线性拟合 para =[1,1];%假定一个a 和b的初值,越接近真实的越好 Func = @(para,t)(para(1)./t(:,1)+para(2).*t(:,2)); result= nlinfit(X,Y,Func,par...

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