24小时热门版块排行榜    

查看: 2330  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wh1125

木虫 (正式写手)

[求助] 一个1stOpt的初级问题:复数方程组计算问题(有代码)

初学1stOpt,以下1stOpt代码中有5个方程,求解5个变量K00,K01,K02,K03,beta0,其中i已定义为虚数单位,但是运行结果中,怎么总把i当成一个变量,有结果i=0。怎么样才能正确计算我要求解的那5个变量呢?请大神指点啊!

ComplexStr = i;
Constant K=2*pi/(0.88e-6),n0=3.4,n1=3.6,n2=3.4,n3=3.55,t=1e-6,g=0.2e-6;
Function
K00^2=K^2*n0^2-beta0^2;
K01^2=K^2*n1^2-beta0^2;
K02^2=K^2*n2^2-beta0^2;
K03^2=K^2*n3^2-beta0^2;
K03*(K01^2+K00*K02)*tan(K01*(t-g))+K01*(K03^2+K00*K02)*tan(K03*g)+i*(K01*K03*(K00+K02)-(K00*K03^2+K01^2*K02)*tan(K03*g)*tan(K01*(t-g)))=0;

运行结果:
k00: 3.28705929208081
beta0: 4.83100198907778
k01: 2.24062444176525
k02: 4.41804389934987
k03: 3.40456867241301
i: 0
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

wh1125

木虫 (正式写手)

引用回帖:
2楼: Originally posted by feixiaolin at 2013-10-07 19:01:58
删除ComplexStr = i;
拆分
K03*(K01^2+K00*K02)*tan(K01*(t-g))+K01*(K03^2+K00*K02)*tan(K03*g)+i*(K01*K03*(K00+K02)-(K00*K03^2+K01^2*K02)*tan(K03*g)*tan(K01*(t-g)))=0;
成两个式子:
K03*(K01^2+K00*K0 ...

虽然能得出结果,全部是实数解,但是明显解的数量级都不对。
3楼2013-10-07 19:47:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wh1125: 金币+5, ★有帮助, 谢谢你!按照你的办法能够求解,但这样得出的是实数解,没有复数。结果不是我期望的。 2013-10-07 19:45:33
dingd: 金币+1, 谢谢交流。 2013-10-07 21:45:05
删除ComplexStr = i;
拆分
K03*(K01^2+K00*K02)*tan(K01*(t-g))+K01*(K03^2+K00*K02)*tan(K03*g)+i*(K01*K03*(K00+K02)-(K00*K03^2+K01^2*K02)*tan(K03*g)*tan(K01*(t-g)))=0;
成两个式子:
K03*(K01^2+K00*K02)*tan(K01*(t-g))+K01*(K03^2+K00*K02)*tan(K03*g)=0;
(K01*K03*(K00+K02)-(K00*K03^2+K01^2*K02)*tan(K03*g)*tan(K01*(t-g)))=0;
因为负数为零等价于实部为0 & 虚部为0
2楼2013-10-07 19:01:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★
感谢参与,应助指数 +1
csgt0: 金币+1, 谢谢 2013-10-09 10:30:26
不知道一楼的问题中待求参数是实数还是虚数?如果都是实数,2楼的方法很好啊,只是似乎无解。方程是否正确?
4楼2013-10-07 21:49:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wh1125

木虫 (正式写手)

引用回帖:
4楼: Originally posted by dingd at 2013-10-07 21:49:52
不知道一楼的问题中待求参数是实数还是虚数?如果都是实数,2楼的方法很好啊,只是似乎无解。方程是否正确?

方程是正确的。解没有限制,可以是实数解也可以是复数解。
5楼2013-10-07 22:32:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:26 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见