求助:积分级数综合画图问题

yewei317

大家帮忙看看怎么画图,我迷糊了~~~

I(x,f,s)...的图像..

c(f,s)的图像...

1. n = 1:8;
   
2. fmax = 70*pi/180;
   
3. alp = @(s,f) cell2mat(cellfun(@(x)reshape(x,1,1,8),arrayfun(@(s,f)acos(cos(f)*cos(50*pi/180)+sin(f)....
   
4. *cos(s)*sin(50*pi/180).*cos(n*45*pi/180)+sin(f).*sin(s)*...
   
5. sin(50*pi/180).*sin(n*45*pi/180)),s,f,'Uni',0),'Uni',0));
   
6. fu = @(s,f) sin(f).*sum(298*cos(alp(s,f))-27.84*sqrt(1.455^2-sin(alp(s,f).^2)).*(1-alp(s,f)/fmax).^2,3);
   
7. fd = @(s,f) sin(f).*sum((1-alp(s,f)/fmax).^2,3);
   
8. q = quad2d(fu,0,22.5*pi/180,0,fmax)/quad2d(fd,0,22.5*pi/180,0,fmax)

复制代码