| 查看: 744 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
4楼2010-03-03 20:25:15
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
ftp0000
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1822.1
- 帖子: 292
- 在线: 6.6小时
- 虫号: 585530
- 注册: 2008-08-08
- 专业: 凝聚态物性 II ：电子结构
★
zhangzg698958(金币+10):谢谢帮忙！xi是个复数，又在指数中，方程显示不好写！想找到它的实部，然后实部与f画图！高手再帮助看看，帮到底，先行谢过！！！ 2010-03-04 10:42
adu886886(金币+1):谢谢提供意见 2010-03-05 08:50
zhangzg698958(金币+10):谢谢帮忙！xi是个复数，又在指数中，方程显示不好写！想找到它的实部，然后实部与f画图！高手再帮助看看，帮到底，先行谢过！！！ 2010-03-04 10:42
adu886886(金币+1):谢谢提供意见 2010-03-05 08:50
|
我仔细检查了lz的程序，发现主要问题不是逆矩阵的问题，虽然有些情况下逆矩阵求不出来 其实lz就是想解一个方程，是xi的函数，lz编程的思路是这样的 S=f(xi); det(S)=0; 想让程序根据第二句直接给出xi的值，这样是行不通的，因为matlab中det求解是不知道xi是什么的 所以我的建议就是lz花点时间把这个方程显式形式写出来，然后求解，这个好像是个二阶的方程，不难用求根公式求 |
5楼2010-03-04 10:11:25
30