| 查看: 619 | 回复: 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 |
» 猜你喜欢
280求调剂
已经有3人回复
085400 328分 求调剂
已经有7人回复
材料专硕322
已经有12人回复
一志愿华中农微生物,288分,三年实验经历
已经有3人回复
314求调剂
已经有14人回复
调剂 化学 307
已经有11人回复
化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文
已经有10人回复
材料299专硕求调剂
已经有5人回复
083200 初试305分 求调剂
已经有10人回复
332,085601求调剂
已经有3人回复
2楼2014-01-08 21:24:05













回复此楼