Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

mam1231 发表于 2006-6-20 13:47

[求助]我编了一个程序，运行附值出错，大家帮忙！谢谢啊

谢谢大家，程序根据大家的建议做了修改！但红色标出的地方，自己不知道如何处理！数据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:27:04编辑过]

ericlin 发表于 2006-6-20 16:25

a=a0G.*alR;b=alG.*a0R;c=a0R.*alR;d=a0G.*alG;

happy 发表于 2006-6-20 20:32

回复：（mam1231）[求助]我编了一个程序，运行有一个地...

<DIV class=quote>以下是引用mam1231在2006-6-20 13:47:12的发言： a0R=Z0R./ZCR;alR=ZLR./ZCR;a0G=Zs./ZCG;a1G=ZLG./ZCG; %由于角频率w，</DIV> 另一个问题，红色标记部分变量写错了，其中l你写成了数字的1了

wei343 发表于 2006-6-20 22:02

你的ZOG在哪

wei343 发表于 2006-6-20 22:07

你的第23中的ZOG在前面好像没值呀。

mam1231 发表于 2006-6-20 22:37

谢谢各位，根据大家自己又将帖子修正了一下，可还有一个地方，红色标注了！大家再帮忙啊！谢谢了

branch 发表于 2006-6-20 23:22

你把function去掉运行看看a=a0G*alR;b=alG*a0R;c=a0R*alR;d=a0G*alG;%这里开始出错，说是矩阵维数不匹配，但不知道该如何处理，请大家帮忙！ 这样在工作数据空间可以看到数据的维数，好几个是1*999001的数据，相乘用*是不允许的 可以用.*或者转置后乘。

mam1231 发表于 2006-6-20 23:29

嗯用转置！谢谢

mam1231 发表于 2006-6-20 23:33

用转置求取后，该如何提取数据用于下面计算呢！能给一个小例子不？

mam1231 发表于 2006-6-23 12:56

谢谢大家

页: [1]

声振论坛's Archiver

[求助]我编了一个程序，运行附值出错，大家帮忙！谢谢啊

回复：（mam1231）[求助]我编了一个程序，运行有一个地...