| 查看: 1953 | 回复: 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 |
回复此楼» 猜你喜欢
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有164人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt能进行隐函数非线性拟合吗？谁能帮拟合一下新型液体的粘度混合计算公式
已经有5人回复
- 求大神用1stOpt帮忙模拟一个复数拟合。
已经有3人回复
- 急求，常微分方程的参数计算问题，有1stopt的帮忙解答一下
已经有5人回复
- 求大神帮忙用1stOpt 拟合复数数据
已经有13人回复
- 求助，用1stopt解方程组，两次得到的解不唯一？什么原因，大神帮忙
已经有6人回复
- sos!!1stopt解非线性方程组，一直出现常数定义错误的问题！求大神帮助！
已经有12人回复
- 求大侠帮助跑一下1stopt
已经有12人回复
- 哪位朋友有注册版1stOpt软件？急需
已经有5人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- 求高版本1stopt，解决一个问题，微分方程组参数拟合
已经有3人回复
- 1stopt运行程序，总是出现字符函数错误的提示，请高手帮助
已经有5人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 1stopt4阶非线性常微分方程，帮小弟看一下，跪求。
已经有7人回复
- 1stOpt高版本求助！求高版本的帮忙走一下 谢了！
已经有16人回复
- 急切求助！一个mathematica求解一个简单一元非线性方程的问题
已经有4人回复
- 1stOpt求解常微分方程边值问题
已经有5人回复
- 哪位大哥能帮我用1stopt3.0以上版本跑一个拟合程序
已经有11人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3880
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
wh1125: 金币+5, ★有帮助, 谢谢你！按照你的办法能够求解，但这样得出的是实数解，没有复数。结果不是我期望的。 2013-10-07 19:45:33
dingd: 金币+1, 谢谢交流。 2013-10-07 21:45:05
感谢参与，应助指数 +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
wh1125
木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 3078.9
- 红花: 5
- 帖子: 499
- 在线: 432.2小时
- 虫号: 417863
- 注册: 2007-07-02
- 性别: GG
- 专业: 集成电路设计与测试
3楼2013-10-07 19:47:58
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2013-10-07 21:49:52
wh1125
木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 3078.9
- 红花: 5
- 帖子: 499
- 在线: 432.2小时
- 虫号: 417863
- 注册: 2007-07-02
- 性别: GG
- 专业: 集成电路设计与测试
5楼2013-10-07 22:32:19
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★
wh1125: 金币+5, ★有帮助, 谢谢！我的解要求K01和K03不为零，因为后面还有计算，他们要做分母。 2013-10-07 23:13:11
csgt0: 金币+1, 谢谢 2013-10-09 10:30:42
wh1125: 金币+5, ★有帮助, 谢谢！我的解要求K01和K03不为零，因为后面还有计算，他们要做分母。 2013-10-07 23:13:11
csgt0: 金币+1, 谢谢 2013-10-09 10:30:42
|
下面这组参考下： k00.实部: 0 k00.虚部: 7290129.24733499 k01.实部: 4269089.10229991 k01.虚部: 0 k02.实部: 0 k02.虚部: -7290129.24733499 k03.实部: 0 k03.虚部: 0 beta0.实部: 25346940.7278267 beta0.虚部: 0 |
6楼2013-10-07 23:09:46