| 查看: 696 | 回复: 8 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】请各位帮忙看看这个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:1  number_theta+1)% E(j)=abs(sum(j)); % PolarXY=10*log10(E(j)); %end %for j=1:1 number_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 ] |
回复此楼» 猜你喜欢
咱们一起用铁证分析2026国家社科基金中标与否
已经有17人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有47人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有6人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有4人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
5楼2009-05-25 08:35:36
2楼2009-05-23 16:28:31
3楼2009-05-24 13:09:23
4楼2009-05-24 23:04:25
