| 查看: 690 | 回复: 1 | ||
[求助]
matlab请教 已有1人参与
|
|
初学,编了个程序,全通型微环的输入输出波,编好号运行出错,提示“In an assignment A(I) = B, the number of elements in B and I must be the same. Error in fft (line 37) ot1m(m)=Ow*exp(i*hw*t1(l))*dw/2*pi;” 将“ot1m(m)"改为"ot1m"能运行，但感觉结果完全不对啊，求大神指教。代码如下： function fft neff=1.5; R=10e-6; tao=1; t=0.9; c=3e8; TR=2*pi*R*neff./c; F=pi*t^(0.5)/(1-t); Tc=F*TR; lambdac=1.4960e-6; i=sqrt(-1); figure; for n=1:3; subplot(1,3,n); u=[0.25,0.5,1]; tao1=u*Tc; to1=3*tao1; wc=2*pi*c./lambdac; t1=linspace(1,25*4^(n-1),25*4^(n-1)*8); ft1=zeros(25*4^(n-1)*8,1); ot1=zeros(25*4^(n-1)*8,1); for l=1:25*4^(n-1)*8; ft1(l)=exp(-((t1(l)-to1)/tao1).^2)*exp(i*wc*t1(l)); lambda=linspace(1.48e-6,1.51e-6,500); w=zeros(500,1); ot1m=zeros(499,1); for m=1:499; w(m)=2*pi*c./lambda(m); w(m+1)=2*pi*c./lambda(m+1); dw=w(m+1)-w(m); hw=(w(m+1)+w(m))/2; A=exp(-(((hw-wc).*tao1)/2).^2); B=exp(-i*(hw-wc)*to1); Fw=sqrt(pi)*tao1.*A.*B; Hw=(t1(l)-tao*exp(-i*hw*TR))./(t1(l)-tao*t*exp(-i*hw*TR)); Ow=Fw*Hw; ot1m(m)=Ow*exp(i*hw*t1(l))*dw/2*pi; end; ot1(l)=sum(ot1m); end plot(t1,ft1,'r--',t1,ot1,'b-.','LineWidth',2); xlabel('时间（ps）'); legend('输入','输出','Location', 'East'); end |
回复此楼» 猜你喜欢
5个%2F，非固定段22个字母
已经有10人回复
-
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有16人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
国社科通讯评审
已经有6人回复
-
面上的水很深
已经有21人回复
2楼2014-01-08 21:24:05
10