马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
谢谢大家,程序根据大家的建议做了修改!但红色标出的地方,自己不知道如何处理!数据a0R,alR,a0G,alG,是随频率变化得出的一系列值,我想在下面求a,b等时,将上面相同频率下的值一一对应相乘.(例如a=a0G*alR,让a0G,alR在相同频率下的值相乘,附给a,这样就是把一列值附给a)。肯定不能把一连串的乘积结果只附给一个a,请问该如何解决啊?<br>function jisuan<br>% 导体参数<br>L=10;%导线长度<br>LG=0.8e-6;LR=0.8e-6;Lm=0.2e-6;<br>CG=5e-13;CR=5e-13;Cm=0.3e-13;<br>%输入阻抗和负载阻抗<br>Zs=50;Z0R=50;ZLR=50;ZLG=50;<br>%源电压 <br>Us=1;f=1e2:1:1e5;<br>%求解<br>j=sqrt(-1);w=2*pi*f;p=w*sqrt(LG*CG);<br>v=1./sqrt(LG*CG);<br>q=cos(p*L);S=sin(p*L)./(p*L);<br>k=Lm./sqrt(LG*LR); %耦合系数<br>ZCR=v*LR*w*sqrt(1-k^2);ZCG=v*LG*w*sqrt(1-k^2);<br>a0R=Z0R./ZCR;alR=ZLR./ZCR;a0G=Zs./ZCG;alG=ZLG./ZCG; %由于角频率w,得出的数据随w变化<br>a=a0G*alR;b=alG*a0R;c=a0R*alR;d=a0G*alG; <FONT color=#dd2222>%这里开始出错,说是矩阵维数不匹配,但不知道该如何处理,请大家帮忙!<br></FONT>n=(1-a)*(1-b)./((1+c)*(1-d));<br>tG=LG*L./(Zs+ZLG)+(CG+Cm)*L*Zs*ZLG./(Zs+ZLG);<br>tR=LR*L./(Z0R+ZLR)+(CG+Cm)*L*Z0R*ZLR./(Z0R+ZLR);<br>D=q^2-(S*w).^2*tR*tG*(1-k^2*n)+j*w*q*S*(tR+tG);<br><br>g=S./D;h=j*w*ZLR*Lm*L./((Z0R+ZLR)*(ZLG+Zs));<br>i=j*w*Z0R*ZLR*Cm*L*ZLG./((Z0R+ZLR)*(ZLG+Zs));<br>URLUs=g*(h-i);<br>plot(f,abs(URLUs));<br>grid;<br>xlabel('MHZ');ylabel('UR(L)/Us');<br>
[此贴子已经被作者于2006-6-20 23:27:04编辑过]
|