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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 求助！三个未知数 四个方程用1stopt(Anto2Fit)约束优化
北京石油化工学院2026年研究生招生接收调剂公告
111/2返回列表12下一页
查看: 2413  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fugitivesk

铁虫 (初入文坛)

[求助] 求助！三个未知数 四个方程用1stopt(Anto2Fit)约束优化 已有2人参与

计算小白，新来的没啥金币
用的Anto2Fit 5.5试用版（明明三个parameters为啥说超出四个了）
方程都比较长
不知道为什么三个未知数会出现四个结果
第四个结果好像是我给的常数变量或者变参量（偷梁换柱？）
贴上mathematica里写的表达式（这个软件里运行过三个方程三个变量的求解，方程本身没问题）


Constant  cc = 3*10^8;
                 la = cc/10^9;
                 k0 = 2*pi/la;
                 r1 =0.17532;
                 r2 =0.21216;
                 t1 =0.77139;
                 t2 =-0.3507;
                 
Parameter k1[30,60],k2[-2,0],dz[0.5,0.7];
         
Function a=sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz));
b=cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz));
c=cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz));
d=sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz));
kk=2*(k1^2+k2^2);
(r1*c)-(r2*d)+r1*(b*k1-a*k2)/(2*k0)+r1*(a*k0*k2+b*k0*k1)/kk-r2*(a*k1+b*k2)/(2*k0)-r2*(a*k0*k1-b*k0*k2)/kk-(b*k1-a*k2)/(2*k0)+(a*k0*k2+b*k0*k1)/kk;
(r2*c)+(r1*d)+r1*(a*k1+b*k2)/(2*k0)+r1*(a*k0*k1-b*k0*k2)/kk+r2*(b*k1-a*k2)/(2*k0)+r2*(a*k0*k2+b*k0*k1)/kk-(a*k1+b*k2)/(2*k0)-(b*k0*k2-a*k0*k1)/kk;
(t1*c)-(t2*d)+t1*(b*k1-a*k2)/(2*k0)+t1*(a*k0*k2+b*k0*k1)/kk-t2*(a*k1+b*k2)/(2*k0)-t2*(a*k0*k1-b*k0*k2)/kk+2;
(t2*c)+(t1*d)+t1*(a*k1+b*k2)/(2*k0)+t1*(a*k0*k1-b*k0*k2)/kk+t2*(b*k1-a*k2)/(2*k0)+t2*(a*k0*k2+b*k0*k1)/kk;

//方程我用了两种写法第二种太长了

Function
(r1*cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)))-(r2*sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz)))+r1*(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k1-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k2)/(2*k0)+r1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k2+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k1)/2*(k1^2+k2^2)-r2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k1+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k2)/(2*k0)-r2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k1-cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k2)/2*(k1^2+k2^2)-(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k1-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k2)/(2*k0)+(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k2+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k1)/2*(k1^2+k2^2);

(r2*cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)))+(r1*sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz)))+r1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k1+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k2)/(2*k0)+r1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k1-cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k2)/2*(k1^2+k2^2)+r2*(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k1-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k2)/(2*k0)+r2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k2+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k1)/2*(k1^2+k2^2)-(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k1+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k2)/(2*k0)-(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k2-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k1)/2*(k1^2+k2^2);

(t1*cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)))-(t2*sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz)))+t1*(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k1-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k2)/(2*k0)+t1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k2+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k1)/2*(k1^2+k2^2)-t2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k1+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k2)/(2*k0)-t2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k1-cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k2)/2*(k1^2+k2^2)+2;

(t2*cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)))+(t1*sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz)))+t1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k1+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k2)/(2*k0)+t1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k1-cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k2)/2*(k1^2+k2^2)+t2*(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k1-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k2)/(2*k0)+t2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k2+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k1)/2*(k1^2+k2^2);

求助！三个未知数 四个方程用1stopt(Anto2Fit)约束优化
22.JPG@月只蓝
回复此楼

» 猜你喜欢
1楼 2020-06-24 18:48:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shikang999

新虫 (小有名气)

【答案】应助回帖

★
fugitivesk: 金币+1, ★有帮助 2020-06-28 19:26:45
按你第一种写法, 如果最后四个方程是等于0的约束的话，求得一个如下的解(其中四个方程的误差平方和为2.07045657911154E-05）（注意，下面不是1stopt的结果，所以不保证和1stopt的结果一致）
k1 = 48.5939462992421,  k2 = -0.161774554371271,  dz = 0.588775267394145
赞一下 回复此楼
我还是那个我，过去如此，未来亦如此！
2楼2020-06-28 15:36:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fugitivesk

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by shikang999 at 2020-06-28 15:36:14
按你第一种写法, 如果最后四个方程是等于0的约束的话，求得一个如下的解(其中四个方程的误差平方和为2.07045657911154E-05）（注意，下面不是1stopt的结果，所以不保证和1stopt的结果一致）
k1 = 48.5939462992421 ...

请问求解思路是什么，用matlab还是mathemetica求解的吗？
赞一下 回复此楼
3楼2020-06-28 19:11:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fugitivesk

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by shikang999 at 2020-06-28 15:36:14
按你第一种写法, 如果最后四个方程是等于0的约束的话，求得一个如下的解(其中四个方程的误差平方和为2.07045657911154E-05）（注意，下面不是1stopt的结果，所以不保证和1stopt的结果一致）
k1 = 48.5939462992421 ...

我原先就只是四选三作为方程组求解，感觉不能保证误差。四个方程的约束用什么方法比较合适呢?

发自小木虫Android客户端
赞一下 回复此楼
4楼2020-06-28 19:24:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shikang999

新虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
独孤神宇: 金币+5, 鼓励交流 2020-06-28 20:41:18
fugitivesk: 金币+7, ★★★很有帮助 2020-06-28 21:14:32
引用回帖:
3楼: Originally posted by fugitivesk at 2020-06-28 19:11:15
请问求解思路是什么，用matlab还是mathemetica求解的吗？...

1. 主要使用自己写的MathSword软件进行求解。（PS：附件可以下载MathSword，在打开程序界面最上端选择【高级优化】，在优化模块右键选择上传的附件 solve.vb文件，然后点击求解，就能看到你问题的答案了。）

2. 至于你说的方程用什么约束比较合适，不太明白你的意思。因为一个方程本身就是一个约束，有约束就可以构建误差，一般的优化器内部会根据你构建的1个或多个误差进行优化。
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : MathSword.2020.06.27.rar
    • 2020-06-28 20:08:56, 2.97 M
  • 附件 2 : solve.vb
    • 2020-06-28 20:09:09, 6.46 K
我还是那个我，过去如此，未来亦如此！
5楼2020-06-28 20:26:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fugitivesk: 金币+12, ★★★★★最佳答案 2020-06-28 21:13:55
代码写的有问题。

k1: 48.5939435307918
k2: -0.161774549680272
dz: 0.588775301298235

------------------------------------------------

Constant  cc = 3*10^8,
r1 =0.17532,
r2 =0.21216,
t1 =0.77139,
t2 =-0.3507,
la = cc/10^9,
k0 = 2*pi/la;
Parameter k1=[30,60],k2=[-2,0],dz=[0.5,0.7];
ConstStr  kk=2*(k1^2+k2^2),
a=sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz)),
b=cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz)),
c=cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)),
d=sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz));
Function
(r1*c)-(r2*d)+r1*(b*k1-a*k2)/(2*k0)+r1*(a*k0*k2+b*k0*k1)/kk-r2*(a*k1+b*k2)/(2*k0)-r2*(a*k0*k1-b*k0*k2)/kk-(b*k1-a*k2)/(2*k0)+(a*k0*k2+b*k0*k1)/kk=0;
(r2*c)+(r1*d)+r1*(a*k1+b*k2)/(2*k0)+r1*(a*k0*k1-b*k0*k2)/kk+r2*(b*k1-a*k2)/(2*k0)+r2*(a*k0*k2+b*k0*k1)/kk-(a*k1+b*k2)/(2*k0)-(b*k0*k2-a*k0*k1)/kk=0;
(t1*c)-(t2*d)+t1*(b*k1-a*k2)/(2*k0)+t1*(a*k0*k2+b*k0*k1)/kk-t2*(a*k1+b*k2)/(2*k0)-t2*(a*k0*k1-b*k0*k2)/kk+2=0;
(t2*c)+(t1*d)+t1*(a*k1+b*k2)/(2*k0)+t1*(a*k0*k1-b*k0*k2)/kk+t2*(b*k1-a*k2)/(2*k0)+t2*(a*k0*k2+b*k0*k1)/kk=0;
赞一下 回复此楼
数值计算
6楼2020-06-28 20:40:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fugitivesk

铁虫 (初入文坛)
感谢大神

发自小木虫Android客户端
回复此楼
7楼2020-06-28 21:15:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
1stOpt正确代码如下，4个方程3个未知数，只能有近似解。
CODE:
Constant  cc = 3*10^8,
         la = cc/10^9,
         k0 = 2*pi/la,
         r1 =0.17532,
         r2 =0.21216,
         t1 =0.77139,
         t2 =-0.3507;

Parameter k1=[30,60],k2=[-2,0],dz=[0.5,0.7];

ConstStr  a=sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz)),
        b=cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz)),
        c=cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)),
        d=sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz)),
        kk=2*(k1^2+k2^2);
Function
(r1*c)-(r2*d)+r1*(b*k1-a*k2)/(2*k0)+r1*(a*k0*k2+b*k0*k1)/kk-r2*(a*k1+b*k2)/(2*k0)-r2*(a*k0*k1-b*k0*k2)/kk-(b*k1-a*k2)/(2*k0)+(a*k0*k2+b*k0*k1)/kk;
(r2*c)+(r1*d)+r1*(a*k1+b*k2)/(2*k0)+r1*(a*k0*k1-b*k0*k2)/kk+r2*(b*k1-a*k2)/(2*k0)+r2*(a*k0*k2+b*k0*k1)/kk-(a*k1+b*k2)/(2*k0)-(b*k0*k2-a*k0*k1)/kk;
(t1*c)-(t2*d)+t1*(b*k1-a*k2)/(2*k0)+t1*(a*k0*k2+b*k0*k1)/kk-t2*(a*k1+b*k2)/(2*k0)-t2*(a*k0*k1-b*k0*k2)/kk+2;
(t2*c)+(t1*d)+t1*(a*k1+b*k2)/(2*k0)+t1*(a*k0*k1-b*k0*k2)/kk+t2*(b*k1-a*k2)/(2*k0)+t2*(a*k0*k2+b*k0*k1)/kk;

结果：
CODE:
k1: 48.5939431819858
k2: -0.161774549476569
dz: 0.588775305571065
赞一下 回复此楼
8楼2020-06-28 22:22:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fugitivesk

铁虫 (初入文坛)
引用回帖:
5楼: Originally posted by shikang999 at 2020-06-28 20:26:02
1. 主要使用自己写的MathSword软件进行求解。（PS：附件可以下载MathSword，在打开程序界面最上端选择【高级优化】，在优化模块右键选择上传的附件 solve.vb文件，然后点击求解，就能看到你问题的答案了。）

2. ...

感谢，回头参考参考

发自小木虫Android客户端
赞一下 回复此楼
9楼2020-06-28 23:01:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fugitivesk

铁虫 (初入文坛)
引用回帖:
8楼: Originally posted by dingd at 2020-06-28 22:22:03
1stOpt正确代码如下，4个方程3个未知数，只能有近似解。

Constant  cc = 3*10^8,
         la = cc/10^9,
         k0 = 2*pi/la,
         r1 =0.17532,
         r2 =0.21216,
         t1 =0.77139,
   ...

感谢

发自小木虫Android客户端
回复此楼
10楼2020-06-28 23:01:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fugitivesk 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学学硕311分求调剂 +11 梁富贵险中求 2026-04-04 13/650 2026-04-06 07:24 by houyaoxu
[考研] 机械专硕274求调剂，不挑专业学校 +6 泛泛2333 2026-04-05 7/350 2026-04-05 23:20 by chyhaha
[考研] （调剂）一志愿报考哈尔滨工业大学0857资源与环境专业378分考生 +7 狠狠加油 2026-04-05 7/350 2026-04-05 22:31 by dongzh2009
[考研] 306分材料与化工求调剂 +7 黎吧啦啦你很有� 2026-04-03 7/350 2026-04-05 17:18 by Hdyxbekcb
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[考研] 材料专硕322分 +11 哈哈哈吼吼吼哈 2026-04-02 11/550 2026-04-04 23:37 by 永字号
[考研] 调剂 +11 JLLLLLLLLLL 2026-04-03 11/550 2026-04-04 22:21 by hemengdong
[考研] 298求调剂 +5 zzz，，r 2026-04-02 8/400 2026-04-04 19:55 by 蓝云思雨
[考研] 一志愿南农090401，268，求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[考研] 考研调剂 +4 zybz冲冲冲 2026-04-03 6/300 2026-04-04 13:08 by zybz冲冲冲
[考研] 0710生物学336分求调剂 +6 kiyy 2026-04-01 8/400 2026-04-04 10:10 by kiyy
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[考研] 一志愿北京交通大学材料工程总分358 +4 cs0106 2026-04-03 4/200 2026-04-03 13:41 by 百灵童888
[考研] 一志愿华东理工大学，080500学硕，317分，求调剂 +13 s1145 2026-03-31 15/750 2026-04-03 11:44 by msi123
[考研] 348求调剂 +11 zzzzyk123 2026-04-01 11/550 2026-04-02 16:52 by Wang200018
[考研] 272求调剂，接受跨专业调剂！ +4 闲鱼卢 2026-03-31 4/200 2026-04-02 11:18 by guyan1000
[考研] 07生物学求调剂 一志愿同济大学359分 +3 LAMC. 2026-03-30 3/150 2026-04-02 10:26 by 18828373951
[考研] 0855机械初试280求调剂 +3 kazenotori 2026-03-31 3/150 2026-04-01 10:08 by JourneyLucky
[考研] 本2一志愿C9-333分，材料科学与工程，求调剂 +9 升升不降 2026-03-31 9/450 2026-03-31 18:01 by 无际的草原
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
信息提示
关闭
请填处理意见
关闭