| 查看: 700 | 回复: 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 |
回复此楼» 猜你喜欢
这样的filecode谁见过
已经有13人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
2楼2014-01-08 21:24:05
