| 查看: 621 | 回复: 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 |
回复此楼» 猜你喜欢
0703化学
已经有33人回复
-
296求调剂
已经有7人回复
-
化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文
已经有13人回复
-
材料290调剂
已经有15人回复
-
368化学求调剂
已经有14人回复
-
调剂
已经有5人回复
-
一志愿中国科学院上海有机所，有机化学356分找调剂
已经有3人回复
-
一志愿矿大，材料工程专硕314分，0856可调都可以
已经有6人回复
-
材料专硕初试分332一志愿西北工业大学，
已经有6人回复
-
新一代电子信息294求调剂 不挑学校
已经有6人回复
2楼2014-01-08 21:24:05
