急求，为何运行后surf不显示图像

gh784478335 · 发表于 2019-8-17 09:18

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

syms x  z q m t c
m=0.5;
R=40;
q=200;
e=1;f=1;
Z= []
for X=[-55:2:55];
  for Y=[-55:2:55]
syms r y
if ((X.*X)./3600+(Y.*Y)./400)<1;

Z(e,f)=nan;

f=f+1;
continue
   else
   a=X.*X
   b=Y.*Y

   end
   [r]=solve((a./(R.^2*(1./r+m.*r).^2)+b./(R.^2.*(1./r-m.*r).^2))==1,r);
   [y]=solve((a./(4.*R.^2.*m.*(cos(y)).^2)-b./(4.*R.^2.*m.*(sin(y)).^2))==1,y);

   y=vpa(y);
J=y(y>0);
K=J(J<2*pi);
L=K(imag(K)==0);

r=vpa(r);
G=r(r>0);
H=G(G<1);
I=H(imag(H)==0);

x=0.5*pi;

A=q*real(((2.*exp(2.*i.*x)-m).*(I.*exp(i.*L)).^2-1)./(m.*(I.*exp(i.*L)).^2-1));
B=((q*(m*I^4+(I*exp(i*L))^2)*(I*exp(i*L))^2)./(I^4*(m-(I*exp(i*L))^2./I^4)*(m*(I*exp(i*L))^2-1))*((2*exp(2*i*x)-m+m*(1+m*(I*exp(i*L))^2-2*exp(2*i*x)*(I*exp(i*L))^2)/(m*(I*exp(i*L))^2-1)))+q/(I^2*(m-(I*exp(i*L))^2/I^4))*(exp(-2*i*x)-((3*exp(2*i*x)*(I*exp(i*L))^2+m*exp(2*i*x)-m^2-1)*(I*exp(i*L))^2)/(m*(I*exp(i*L))^2-1)+((exp(2*i*x)*(I*exp(i*L))^2+m*exp(2*i*x)-m^2-1)*(I*exp(i*L))^4)*2*m/(m*(I*exp(i*L))^2-1)^2));

O=B./((I.*exp(i.*L)).^2.*(R.*(m-1./(I.*exp(i.*L)).^2))).*(I.^2.*R.*(m-(I.*exp(i.*L)).^2./I.^4));

D=(A+real(O))./2;
E=(A-real(O))./2;
F=sqrt(D.^2+E.^2);
Z(e,f)=E

f=f+1;
  end
  e=e+1;
  f=1;
end
[d]=[-55:2:55];
[c]=[-55:2:55];
%idy =d.^2./3600+c.^2./400< 1;
%d(idy) = nan;
%c(idy) = nan;

figure(3)
suRF(d,c,Z)
%surf(d,c,Z)
%view([0,1,0]);
shading interp;
colorbar;
%axis equal;
这个函数运行后可以产生Z的矩阵值，但是surf不显示图像，但我将xy的范围定义在-45到45的时候会显示图像，请问是什么问题

Liberallh21 · 发表于 2019-8-19 18:10

没看懂

帐号		自动登录	找回密码
密码			注册

急求，为何运行后surf不显示图像

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1