帮看个matlab程序

来源:百度知道 编辑:UC知道 时间:2024/09/28 15:40:07
a=2;b=2;
x=-a:0.2:a;
y=-b:0.2:b;
for i=1:length(x)
for j=1:length(y);
if x(i)+y(j)>1
z(i,j)=0.5457*exp(-0.75*y(j)^2-3.75*x(i)^2-1.5*x(i));
elseif x(i)+y(j)>-1 & x(i)+y(j)<=1
z(i,j)=0.7575*exp(-y(j)^2-6*x(i)^2);
else z(i,j)=0.5457*exp(-0.75*y(j)^2-3.75*x(i)^2+1.5*x(i));
end
end
end
axis([-a,a,-b,b,min(min(z)),max(max(z))]);
[X,Y]=meshgrid(x,y);
colormap(flipud(winter));
surf(X,Y,z)
运行完后
size(z)=21 126
我就无语了

a=2;b=2; 

x=-a:0.2:a; 

y=-b:0.2:b; 

for i=1:length(x);

    for j=1:length(y); 

        if x(i)+y(j)>1 

            z(i,j)=0.5457*exp(-0.75*y(j)^2-3.75*x(i)^2-1.5*x(i)); 

        elseif x(i)+y(j)>-1 & x(i)+y(j)<=1 

            z(i,j)=0.7575*exp(-y(j)^2-6*x(i)^2); 

        else z(i,j)=0.5457*exp(-0.75*y(j)^2-3.75*x(i)^2+1.5*x(i)); 

        end 

    end 

end 

axis([-a,a,-b,b,min(min(z)),max(max(z))]); 

[X,Y]=meshgrid(x,y); 

colormap(flipud(winter)); 

sur