Question

matlab三维图

程序：
[X,Y]=meshgrid([200:2:300],[-pi/4:pi/100:pi/4]);
M=-pi/4:pi/100:pi/4;
f=sqrt(1-((1600+(X.*cos(M+pi/4)-0.707*X+40*sin(Y)).^2+(X.*sin(M+pi/4)-0.707*X-40*cos(Y)).^2-(0.707*X-40*sin(Y)).^2-(0.707*X+40*cos(Y)).^2)./(80*X)).^2);
Z=(pi/2)/(trapz(M,f));
mesh(X,Y,Z)
运行错误
??? Error using ==> times
Matrix dimensions must agree.
应该是矩阵不对什么的，望大神帮忙改下
函数如下

Answer 1

简单说明两点：

1、表达式中涉及到三个变量：X、Y和M，似乎没有太好的方法避免使用循环（即使有，可能也会把程序搞得更难理解），所以程序改用循环结构。

2、按照你给的表达式计算，结果Z会出现复数（我直接使用你的f表达式，没有仔细检查有没有问题，请你再自己核实一下），所以绘图分别绘制实部和虚部。

代码如下：

[X,Y] = meshgrid([200:2:300],[-pi/4:pi/100:pi/4]);
M=-pi/4:pi/100:pi/4;
for i = 1 : size(X,1)
    for j = 1 : size(X,2)
        x = X(i, j);
        y = Y(i, j);
        f = sqrt(1-((1600+(x.*cos(M+pi/4)-0.707*x+40*sin(y)).^2+...
            (x.*sin(M+pi/4)-0.707*x-40*cos(y)).^2-(0.707*x-...
            40*sin(y)).^2-(0.707*x+40*cos(y)).^2)./(80*x)).^2);
        Z(i, j) = (pi/2)/(trapz(M,f));
    end
end
subplot 121
mesh(X,Y,real(Z))
xlabel('{\itl}_2'); ylabel('\theta');
zlabel('Z'); title('Real part'); 
subplot 122
mesh(X,Y,imag(Z))
xlabel('{\itl}_2'); ylabel('\theta');
zlabel('Z'); title('Imaginary part');

 

关于 计算Z会出现复数的问题，我们取x=200, y=0，看一下f表达式中分母的根号里面1减去的那部分的值：

x=200;
y=0;
W=((1600+(x.*cos(M+pi/4)-0.707*x+40*sin(y)).^2+(x.*sin(M+pi/4)-...
    0.707*x-40*cos(y)).^2-(0.707*x-40*sin(y)).^2-...
    (0.707*x+40*cos(y)).^2)./(80*x));
plot(M,W)

从画出的图很容易看到，W绝对值多数时候大于1，拿1减去W的平方再开方，当然就会得到复数了。

Answer 2

确实是你的矩阵维数不对，稍作修改如下：
有个问题你需要解决，怎么保证开根号里的值是正的，否则你会开根号得复数！
clear;clc;
x = [200:2:300];
y = [-pi/4:pi/100:pi/4];
h = pi/100;
M=-pi/4:pi/100:pi/4;
[X,Y] = meshgrid([200:2:300],[-pi/4:pi/100:pi/4]);
for ii = 1:length(x)
for jj = 1:length(y)
temp1 = 1600+(x(ii).*cos(M+pi/4)-0.707*x(ii)+40*sin(y(jj))).^2;
temp2 = (x(ii).*sin(M+pi/4)-0.707*x(ii)-40*cos(y(jj))).^2;
temp3 = (0.707*x(ii)-40*sin(y(jj))).^2*ones(1,51);
temp4 = (0.707*x(ii)+40*cos(y(jj))).^2*ones(1,51);
temp5 = 80*x(ii)*ones(1,51);
f = sqrt(abs(1-((temp1+temp2-temp3-temp4)./temp5).^2));
%%%% 这里开根号出现复数，我采用加绝对值，只为调试用，你自己修改
z1(ii,jj)=trapz(M,f);
end
end
Z=(pi/2)./z1;
mesh(X,Y,Z)