[交流] 【求助】请各位帮忙看看这个Matlab程序的问题出在哪里？谢谢！【已完成】

主要是实现E=exp^i*psai(j)+exp^i*psai(j)+exp^i*psai(j)+……+exp^i*psai(j)其中psai(j)=dr*cos(theta)+delta(j);其中i表示虚数；j的值为1到n；dr是一个常数；theta是0：360；delta(j)也是一个常数，不过此常数对于不同的j可能不同（所以对应不同的j，psai(j)值也不同），此常数通过load导入，就是通过文件e导入。 按以下程序算出来的结果中，E向量里的值是相同的，但此结果是不对的。问题可能在绿色部分，但我搞不懂（也许下面这个程序的思路本来就是错的！），请各位帮忙分析下，非常感谢！ n=9; %e文件中的数据个数； c=3e8; f=8.58e8; %频率 number_theta=10; %将2*pi分为100段； a=2*pi/number_theta; lambda=c/f; %波长 d=0.29; %间距 dr=2*pi*d/lambda; %电距离 load C:\MATLAB7\work\e.txt %导入初相位，e文件的内容就是n个数据单； delta=e sum=zeros(number_theta+1,1) M=1 for num=1:1:number_theta+1     for j=1:n         psai(j)=dr*cos((num-1)*a)+delta(j)         sum(M)=sum(M)+exp(psai(j)*sqrt(-1))     end     M=M+1; end E=abs(sum) %后面的程序暂时不看； %for j=1:1number_theta+1) %    E(j)=abs(sum(j)); %    PolarXY=10*log10(E(j)); %end %for j=1:1number_theta+1)     %theta=0:a:2*pi     %plot(a*(j-1),PolarXY(j))     %hold on %end %polar(theta,E,'-'); %polar(theta,max(PolarXY,0),'-'); %grid on; %grid minor; [ Last edited by nono2009 on 2009-9-25 at 07:19 ]

回复此楼