应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】50金币求解觉矩阵中未知数画图问题,先行谢过！！！
51/1返回列表
查看: 742  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 zhangzg698958 的 14 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhangzg698958

金虫 (正式写手)

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

» 猜你喜欢
1楼 2010-03-02 12:14:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
★ ★
zhangzg698958(金币+30):谢谢！！！ 2010-03-05 00:17
zhangzg698958(金币+20):再次感谢！！！ 2010-03-05 00:27
adu886886(金币+2):谢谢交流，欢迎常来仿真编程版 2010-03-05 08:50
function A111
clc
clear all
syms x
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;
d1=0.002;d2=0.002;d3=0.002;d4=0.002;d5=0.002;d6=0.002;
f=250:100:20000;
for n=1:length(f)
    k1=2*pi*f(n)/c1; k2=2*pi*f(n)/c2; k3=2*pi*f(n)/c3; k4=2*pi*f(n)/c4; k5=2*pi*f(n)/c5; k6=2*pi*f(n)/c6;
    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(x).*I;
    T=det(S);
    xi=solve(T);
    xi=double(xi);
    xi1(n,: )=real(xi);
    xi2(n,: )=imag(xi);
end
plot(f, xi1)
赞一下(1人) 回复此楼
8楼2010-03-04 23:53:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

ftp0000

木虫 (小有名气)
kuhailangyu:欢迎参与交流 2010-03-03 13:47
zhangzg698958(金币+1):抛砖引玉，谢谢！！！ 2010-03-03 19:13
lz最好不用用这样的附件形式，还要流量呢，如果能用图片形式就好了
赞一下(1人) 回复此楼
2楼2010-03-03 13:34:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangzg698958

金虫 (正式写手)

高手帮忙！！！

我连续修改数十次未有结果，此程序涉及循环、指数、复数、矩阵、隐函数等诸多方面，我实在无能为力！
对于高手而言，此程序有可能过于简单，但对于我，目前实在无办法！
以前在此得到很多帮助，希望高手不吝出手！！！
些许金币，不成敬意！！！

谢谢谢谢-------!!!!!!
赞一下 回复此楼
3楼2010-03-03 19:17:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangzg698958

金虫 (正式写手)

自查！！！

问题有可能出现在逆矩阵中，逆矩阵的运算用到行列式，而行列式无法进行未知数的运算！！！

请高手解决此问题！！！
赞一下 回复此楼
4楼2010-03-03 20:25:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 03:26 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:29 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:18 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-30 01:58 by ZPa0EcMwuECS
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:14 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 01:04 by ZPa0EcMwuECS
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 00:52 by ZPa0EcMwuECS
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
信息提示
关闭
请填处理意见
关闭