24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 zhangzg698958 的 14 个金币，回帖就立即获得 1 个金币，每人有 1 次机会

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhangzg698958

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 1088
帖子: 711
在线: 112.1小时
虫号: 722666
注册: 2009-03-14
专业: 人工晶体

[交流] 【求助】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 ]

回复此楼

» 猜你喜欢

体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有9人回复
过年走亲戚时感受到了所开私家车的鄙视链已经有9人回复
今年春晚有几个节目很不错，点赞！已经有10人回复
情人节自我反思：在爱情中有过遗憾吗？已经有10人回复
基金正文30页指的是报告正文还是整个申请书已经有5人回复

1楼 2010-03-02 12:14:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

我仔细检查了lz的程序，发现主要问题不是逆矩阵的问题，虽然有些情况下逆矩阵求不出来

其实lz就是想解一个方程，是xi的函数，lz编程的思路是这样的

S=f(xi);
det(S)=0;

想让程序根据第二句直接给出xi的值，这样是行不通的，因为matlab中det求解是不知道xi是什么的

所以我的建议就是lz花点时间把这个方程显式形式写出来，然后求解，这个好像是个二阶的方程，不难用求根公式求