斑竹请不要删！用MATALB编了一个计算程序，有一个问题，请大家帮忙！

mam1231 · 发表于 2006-6-20 23:02

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

斑竹和高手帮忙啊！也可给点其他程序解决这样问题的建议！已经标注错误的地方！ 数据a0R，alR，a0G，alG，是随频率变化得出的一系列值，我想在下面求a,b等时，是用上面相同频率点下的值一一对应相乘.(例如a=a0G*alR，让a0G，alR在相同频率下的值相乘，附给a，这样就是把一列值附给a）。肯定不能把一连串的乘积结果只附给一个a，请问该如何解决啊？非常感谢！ function jisuan % 导体参数 L=10;%导线长度 LG=0.8e-6;LR=0.8e-6;Lm=0.2e-6; CG=5e-13;CR=5e-13;Cm=0.3e-13; %输入阻抗和负载阻抗 Zs=50;Z0R=50;ZLR=50;ZLG=50; %源电压 Us=1;f=1e2:1:1e5; %求解 j=sqrt(-1);w=2*pi*f;p=w*sqrt(LG*CG); v=1./sqrt(LG*CG); q=cos(p*L);S=sin(p*L)./(p*L); k=Lm./sqrt(LG*LR); %耦合系数 ZCR=v*LR*w*sqrt(1-k^2);ZCG=v*LG*w*sqrt(1-k^2); a0R=Z0R./ZCR;alR=ZLR./ZCR;a0G=Zs./ZCG;alG=ZLG./ZCG; %由于角频率w，得出的数据随w变化 a=a0G*alR;b=alG*a0R;c=a0R*alR;d=a0G*alG; %这里开始出错，说是矩阵维数不匹配，但不知道该如何处理，请大家帮忙！ n=(1-a)*(1-b)./((1+c)*(1-d)); tG=LG*L./(Zs+ZLG)+(CG+Cm)*L*Zs*ZLG./(Zs+ZLG); tR=LR*L./(Z0R+ZLR)+(CG+Cm)*L*Z0R*ZLR./(Z0R+ZLR); D=q^2-(S*w).^2*tR*tG*(1-k^2*n)+j*w*q*S*(tR+tG); g=S./D;h=j*w*ZLR*Lm*L./((Z0R+ZLR)*(ZLG+Zs)); i=j*w*Z0R*ZLR*Cm*L*ZLG./((Z0R+ZLR)*(ZLG+Zs)); URLUs=g*(h-i); plot(f,abs(URLUs)); grid; xlabel('MHZ');ylabel('UR(L)/Us');

[此贴子已经被作者于2006-6-20 23:13:57编辑过]

lvyehuaigu · 发表于 2006-6-20 23:41

我看了你的程序，变量比较多，错误的原因主要是矩阵乘法的维数不对，我觉得你对matlab语言的.*和*及.^和^、./ 和/的区别没有搞清楚，凡是加点的是向量的相应运算，就是对应元素进行相关运算，不加点就是按照矩阵法则进行的运算，和高等代数（线形代数）的矩阵运算法则对应。我都把你的程序错误的语句改为了.*或./，结果就没有错误了，但是可能和你的本意不是很一样，因为我不是很懂你的程序要实现的具体实际意义，但是你的错误就是矩阵*或^的维数不对，你自己理解了matlab的关于向量及矩阵的四则运算法则后应该很容易解决的。 function jisuan % 导体参数 L=10;%导线长度 LG=0.8e-6;LR=0.8e-6;Lm=0.2e-6; CG=5e-13;CR=5e-13;Cm=0.3e-13; %输入阻抗和负载阻抗 Zs=50;Z0R=50;ZLR=50;ZLG=50; %源电压 Us=1;f=1e2:1:1e5; %求解 j=sqrt(-1);w=2*pi*f;p=w*sqrt(LG*CG); v=1./sqrt(LG*CG); q=cos(p*L);S=sin(p*L)./(p*L); k=Lm./sqrt(LG*LR); %耦合系数 ZCR=v*LR*w*sqrt(1-k^2);ZCG=v*LG*w*sqrt(1-k^2); a0R=Z0R./ZCR;alR=ZLR./ZCR;a0G=Zs./ZCG;alG=ZLG./ZCG; %由于角频率w，得出的数据随w变化 a=a0G*alR'; b=alG*a0R'; c=a0R*alR'; d=a0G*alG'; %这里开始出错，说是矩阵维数不匹配，但不知道该如何处理，请大家帮忙！ n=(1-a)*(1-b)./((1+c)*(1-d)); tG=LG*L./(Zs+ZLG)+(CG+Cm)*L*Zs*ZLG./(Zs+ZLG); tR=LR*L./(Z0R+ZLR)+(CG+Cm)*L*Z0R*ZLR./(Z0R+ZLR); D=q.^2-(S.*w).^2.*tR.*tG.*(1-k.^2.*n)+j.*w.*q.*S.*(tR+tG); g=S./D;h=j*w*ZLR*Lm*L./((Z0R+ZLR)*(ZLG+Zs)); i=j*w*Z0R*ZLR*Cm*L*ZLG./((Z0R+ZLR)*(ZLG+Zs)); URLUs=g.*(h-i); plot(f,abs(URLUs)); grid; xlabel('MHZ');ylabel('UR(L)/Us');

cc800 · 发表于 2006-6-21 22:37

我也看了同意楼上的观点!~~~~

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 斑竹请不要删！用MATALB编了一个计算程序，有一个问题，请大家帮忙！

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。