24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zzucbx

新虫 (初入文坛)

应助: 2 (幼儿园)
金币: 165
散金: 60
帖子: 23
在线: 9.4小时
虫号: 1814752
注册: 2012-05-13
专业: 光学

[求助] MATLAB 程序优化

这里一个程序，算了两天了，还是没画出来图形，高手能不能优化一下，让我得到图形，万分感谢！
format long
warning off
f=zeros(9,9);
m=-4:1:4;
n=-4:1:4;
a=2.5;
b=2.5;
c=0.5;
d=20;
l=60;
e=(d/c)^2;
o=c/d;
for p=1:9;
for q=1:9;
f(p,q)=(-(a^3*b)/4-(a*b^3)/4+2*a*b*l^2-3*a^3*b*m(p)^2-3*a*b^3*n(q)^2)/(2*l^4);
end
end;
syms x y s k;
z=e*f(1,1)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(1,2)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(1,3)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(1,4)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(1,5)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(1,6)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(1,7)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(1,8)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(1,9)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(2,1)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(2,2)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(2,3)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(2,4)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(2,5)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(2,6)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(2,7)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(2,8)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(2,9)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(3,1)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(3,2)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(3,3)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(3,4)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(3,5)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(3,6)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(3,7)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(3,8)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(3,9)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(4,1)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(4,2)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(4,3)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(4,4)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(4,5)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(4,6)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(4,7)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(4,8)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(4,9)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(5,1)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(5,2)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(5,3)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(5,4)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(5,5)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(5,6)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(5,7)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(5,8)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(5,9)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(6,1)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(6,2)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(6,3)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(6,4)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(6,5)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(6,6)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(6,7)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(6,8)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(6,9)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(7,1)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(7,2)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(7,3)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(7,4)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(7,5)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(7,6)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(7,7)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(7,8)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(7,9)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(8,1)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(8,2)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(8,3)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(8,4)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(8,5)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(8,6)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(8,7)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(8,8)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(8,9)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(9,1)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(9,2)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(9,3)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(9,4)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(9,5)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(9,6)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(9,7)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(9,8)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(9,9)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o));
h=int(z^2*cos(1i*k*x),x,-sqrt(1-y^2)+s/2,sqrt(1-y^2)-s/2);
j=subs(h,y,-sqrt(1-(s/2)^2)/100:sqrt(1-(s/2)^2)/10000:sqrt(1-(s/2)^2)/100);
g=sum(j(