| syms k n a L deta m0 m1 s zs l aj x y ax ay az u0 k=1;n=6;a=0.25;L=0.5;deta=0.5136;m0=-0.229;m1=0.618;s=0.85;zs=1.25;l=6.5227;aj=0.175; f1=sin(x)*(l/2); f2=1.25-cos(x)*cos(y)*(l/2); f3=-f1; f4=1.25-cos(x)*cos(y)*(-l/2); A=f1^2+f2^2;B=f3^2+f4^2; Q11=[f1*B^(3/2)-f3*A^(3/2) f2*B^(3/2)-f4*A^(3/2)] Q12=[ax;az] Q1=Q11*Q12 Q2=(A*B)^(3/2) q0=Q1/Q2 b=20;c=20; [X,Y]=meshgrid(linspace(-pi/9,pi/9,b),linspace(-pi/9,pi/9,c)); Z=-m0*u0/(a*pi).*D; surf(X,Y,Z); xlabel('\beta(rad)'),ylabel('\alf(rad)'),zlabel('B_{1r}(T)') 这个程序哪出错了,帮忙改一下,怎么把Z换成一个矩阵 |
| 欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) | Powered by Discuz! X2.5 |