版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

cdg234

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 6.4
散金: 30
帖子: 24
在線: 22.6小時
蟲號: 970387
注冊: 2010-03-14
專業(yè): 激光

[求助] matlab請教已有1人參與

初學(xué),編了個程序,全通型微環(huán)的輸入輸出波,編好號運行出錯,提示“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"能運行，但感覺結(jié)果完全不對啊，求大神指教。代碼如下：

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

回復(fù)此樓

» 猜你喜歡

0856材料專碩353求調(diào)劑已經(jīng)有4人回復(fù)
山東省面上項目限額評審已經(jīng)有4人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有5人回復(fù)
生物學(xué)071000 329分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿華中科技大學(xué)071000，求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)一志愿985，分數(shù)349求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)

1樓 2014-01-08 16:24:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bill110100

銅蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 117.5
帖子: 7
在線: 13.1小時
蟲號: 2799633
注冊: 2013-11-14
專業(yè): 基礎(chǔ)物理學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
cdg234: 金幣+5 2014-01-09 09:42:37

問題是你的OW是1*3的矩陣，
后面的ot1m（m）里面是矩陣乘法，最后的結(jié)果也是一個矩陣，不能保存在矩陣的一個元里！

贊一下

回復(fù)此樓

2樓2014-01-08 21:24:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 cdg234 的主題更新