| 查看: 745 | 回复: 7 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 zhangzg698958 的 14 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】50金币求解觉矩阵中未知数画图问题,先行谢过！！！
|
|||
|
已知： ， I=[1 0;0 1]; T1=[1 1;1 -1]; T2=[exp(-j*k2*d2) exp(j*k2*d2);Z1*exp(-j*k2*d2)/Z2 -Z1*exp(j*k2*d2)/Z2]; T3=[exp(-j*k3*d3) exp(j*k3*d3);Z2*exp(-j*k3*d3)/Z3 -Z2*exp(j*k3*d3)/Z3]; T4=[exp(-j*k4*d4) exp(j*k4*d4);Z3*exp(-j*k4*d4)/Z4 -Z3*exp(j*k4*d4)/Z4]; T5=[exp(-j*k5*d5) exp(j*k5*d5);Z4*exp(-j*k5*d5)/Z5 -Z4*exp(j*k5*d5)/Z5]; T6=[exp(-j*k6*d6) exp(j*k6*d6);Z5*exp(-j*k6*d6)/Z6 -Z5*exp(j*k6*d6)/Z6]; 其中： k1= 2*pi*f/c1; k2=2*pi*f/c2; k3=2*pi*f/c3; k4=2*pi*f/c4; k5=2*pi*f/c5; k6=2*pi*f/c6; d1=0.002;d2=0.002;d3=0.002;d4=0.002;d5=0.002;d6=0.002; c1=5960; c2=96;c3=5960;c4=96;c5=5960;c6=96; rho1=7850; rho2=1300; rho3=7850; rho4=1300; rho5=7850; rho6=1300; Z1= rho1*c1; Z2= rho2*c2; Z3= rho3*c3; Z4=rho4*c4; Z5=rho5*c5; Z6=rho6*c6; f=250:1:20000 令S=det( T6^-1*T1*T5^-1*T1*T4^-1*T1*T3^-1*T1*T2^-1*T1-exp(xi)*I)=0时，如何画出f与real(xi)的曲线。 以下是我编的程序： close all clear all syms xi for f=250:1:20000; c1=5960; c2=96;c3=5960;c4=96;c5=5960;c6=96; rho1=7850; rho2=1300; rho3=7850; rho4=1300; rho5=7850; rho6=1300; Z1= rho1*c1; Z2= rho2*c2; Z3= rho3*c3; Z4=rho4*c4; Z5=rho5*c5; Z6=rho6*c6; k1= 2*pi*f/c1; k2=2*pi*f/c2; k3=2*pi*f/c3; k4=2*pi*f/c4; k5=2*pi*f/c5; k6=2*pi*f/c6; d1=0.002;d2=0.002;d3=0.002;d4=0.002;d5=0.002;d6=0.002; T1=[1 1;1 -1];I=[1 0;0 1]; T2=[exp(-j*k2*d2) exp(j*k2*d2);Z1*exp(-j*k2*d2)/Z2 -Z1*exp(j*k2*d2)/Z2]; T3=[exp(-j*k3*d3) exp(j*k3*d3);Z2*exp(-j*k3*d3)/Z3 -Z2*exp(j*k3*d3)/Z3]; T4=[exp(-j*k4*d4) exp(j*k4*d4);Z3*exp(-j*k4*d4)/Z4 -Z3*exp(j*k4*d4)/Z4]; T5=[exp(-j*k5*d5) exp(j*k5*d5);Z4*exp(-j*k5*d5)/Z5 -Z4*exp(j*k5*d5)/Z5]; T6=[exp(-j*k6*d6) exp(j*k6*d6);Z5*exp(-j*k6*d6)/Z6 -Z5*exp(j*k6*d6)/Z6]; S= T6^-1*T1*T5^-1*T1*T4^-1*T1*T3^-1*T1*T2^-1*T1-exp(xi)*I; det(S)=0; end f=250:1:20000; plot(f,real(xi)) 未能得到希望结果，敬请高手指正。附件中为详细说明。 [ Last edited by zhangzg698958 on 2010-3-5 at 00:30 ] |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有18人回复
-
面上意见出来了
已经有17人回复
-
找导师
已经有5人回复
-
国社科又开始会评了，不知道这次命运如何
已经有3人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
高手帮忙！！！
|
又进行了一些列修改，还是存在问题！！！在下面的程序中，T1='T'; xi=solve('T');可解出两个数值，而画图只需要一个数值，到底错在哪里？是solve的问题还是plot的问题？？？高手帮忙！！！我是学材料的，matlab半路出家，易学难精啊！而且，指望这个程序毕业哪！！！ close all clear all syms xi for f=250:100:20000; c1=5960; c2=96;c3=5960;c4=96;c5=5960;c6=96; rho1=7850; rho2=1300; rho3=7850; rho4=1300; rho5=7850; rho6=1300; Z1= rho1*c1; Z2= rho2*c2; Z3= rho3*c3; Z4=rho4*c4; Z5=rho5*c5; Z6=rho6*c6; k1= 2*pi*f/c1; k2=2*pi*f/c2; k3=2*pi*f/c3; k4=2*pi*f/c4; k5=2*pi*f/c5; k6=2*pi*f/c6; d1=0.002;d2=0.002;d3=0.002;d4=0.002;d5=0.002;d6=0.002; T1=[1 1;1 -1];I=[1 0;0 1]; T2=[exp(-j*k2*d2) exp(j*k2*d2);Z1*exp(-j*k2*d2)/Z2 -Z1*exp(j*k2*d2)/Z2]; T3=[exp(-j*k3*d3) exp(j*k3*d3);Z2*exp(-j*k3*d3)/Z3 -Z2*exp(j*k3*d3)/Z3]; T4=[exp(-j*k4*d4) exp(j*k4*d4);Z3*exp(-j*k4*d4)/Z4 -Z3*exp(j*k4*d4)/Z4]; T5=[exp(-j*k5*d5) exp(j*k5*d5);Z4*exp(-j*k5*d5)/Z5 -Z4*exp(j*k5*d5)/Z5]; T6=[exp(-j*k6*d6) exp(j*k6*d6);Z5*exp(-j*k6*d6)/Z6 -Z5*exp(j*k6*d6)/Z6]; S= T6^-1*T1*T5^-1*T1*T4^-1*T1*T3^-1*T1*T2^-1*T1-exp(xi)*I; end T=det(S); T1='T'; xi=solve('T'); xi1=real(xi); xi2=imag(xi); f=250:100:20000; plot(f, xi1) |
6楼2010-03-04 21:50:59
ftp0000
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1822.1
- 帖子: 292
- 在线: 6.6小时
- 虫号: 585530
- 注册: 2008-08-08
- 专业: 凝聚态物性 II ：电子结构
2楼2010-03-03 13:34:35
3楼2010-03-03 19:17:32
4楼2010-03-03 20:25:15
30