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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 1stOpt编程求助
41/1返回列表
查看: 749  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

CarrieQh

银虫 (正式写手)

[求助] 1stOpt编程求助 已有1人参与

各位虫友，小妹初学1stOpt，发现这个软件上手很快，也很好用，准备用其处理我的拟合问题。但是我的拟合条件非常复杂，我又是初学，所以不知道该如何处理了。希望在这里得到大家的帮助，因为之前确实用MATLAB处理了好久，也没有满意的结果，不是满足不了约束条件，就是拟合精度不够。。小妹先谢谢大家了。。
我的问题是这样的：
我有一组800行2列的数据，要拟合的就是这组数（变量是k和S）
拟合公式是确定的，就是Function S=((1+p1*i*3.14*(k)/2500+p2*(i*3.14*(k)/2500)^2+p3*(i*3.14*(k)/2500)^3+p4*(i*3.14*(k))/2500)^4+p5*(i*3.14*(k)/2500)^5+p6*(i*3.14*(k)/2500)^6)/(1+p7*i*3.14*(k)/2500+p8*(i*3.14*(k)/2500)^2+p9*(i*3.14*(k)/2500)^3+p10*(i*3.14*(k)/2500)^4+p11*(i*3.14*(k)/2500)^5);
参数是p，共11个，即p1,p2,p3...p11
我的限制条件是分母这个多项式等于0时，这个方程的根的实部全部小于0.
我知道在MATLAB里面可以用roots来求多项式的根，不知道在1stOpt里面，有没有相关的命令呢？

为了方便，把我的数据也写上吧（空间有限，写150行）
Data;
1        0.279479639197136 + 0.0384417832209621i
2        0.281727699282176 + 0.0480552970486635i
3        0.285786220529823 + 0.0563204916752128i
4        0.289050452923830 + 0.0638514374708868i
5        0.294537097176165 + 0.0699183664359744i
6        0.298876314491868 + 0.0749050754867195i
7        0.301760869716538 + 0.0791306803416643i
8        0.303687224737348 + 0.0829270503941455i
9        0.307079398270538 + 0.0859237681025002i
10        0.309997440590737 + 0.0886875536805402i
11        0.310508772572096 + 0.0912725857065100i
12        0.313152809226795 + 0.0934490256951712i
13        0.313109302543649 + 0.0956300178119395i
14        0.314883975187872 + 0.0981752336734337i
15        0.313372638498908 + 0.100284800783650i
16        0.314173716343903 + 0.103059850285706i
17        0.312993578407729 + 0.104883314837608i
18        0.312871615363587 + 0.107843241112237i
19        0.311584953669249 + 0.110556990263263i
20        0.312378452768188 + 0.113118190574925i
21        0.311606870083760 + 0.115796823293006i
22        0.308493066799768 + 0.119154490884026i
23        0.307815074075734 + 0.122138321663898i
24        0.306756945265756 + 0.125472886855160i
25        0.304004301363881 + 0.129284440862274i
26        0.302466303815432 + 0.132723513548496i
27        0.301657925476687 + 0.136782510234684i
28        0.299724992503746 + 0.141025766593126i
29        0.298741447214539 + 0.145334827625565i
30        0.294511722296633 + 0.149448981026675i
31        0.292140932344386 + 0.154378189423179i
32        0.291731953424778 + 0.159432426112257i
33        0.290700381154457 + 0.163997096516899i
34        0.289933953186610 + 0.169110402164468i
35        0.285759234004963 + 0.174515468750946i
36        0.285301850651697 + 0.180226689402304i
37        0.285318778978834 + 0.186243354943311i
38        0.282780641586369 + 0.191700451606933i
39        0.282224631865583 + 0.197410777269511i
40        0.282357794303744 + 0.203542493119098i
41        0.281573705026058 + 0.209433408920143i
42        0.280583905140873 + 0.215254780847297i
43        0.279992345740081 + 0.221340486264525i
44        0.279482701307005 + 0.227669759794188i
45        0.278772864533912 + 0.233616891411661i
46        0.281313625702324 + 0.239533890241768i
47        0.279742222322056 + 0.245538085531388i
48        0.281383078474079 + 0.250792563746291i
49        0.281713036994699 + 0.256671941670209i
50        0.282327173845756 + 0.261924190914206i
51        0.283038674259945 + 0.267391149378866i
52        0.285103789430794 + 0.271830478080321i
53        0.286143322610054 + 0.276844159095643i
54        0.287325442040975 + 0.281501531551064i
55        0.287796710509510 + 0.286214519189313i
56        0.289349569061817 + 0.290191669956809i
57        0.289986384207662 + 0.294160869712012i
58        0.289241719036391 + 0.298050277467399i
59        0.291670816495331 + 0.301824335846840i
60        0.290139194320732 + 0.305102677733350i
61        0.290246071733169 + 0.308383318810720i
62        0.290936893443977 + 0.311939648600401i
63        0.291768489044532 + 0.315249899450593i
64        0.290785248754759 + 0.318469217233985i
65        0.290748913170460 + 0.321352077786297i
66        0.288876259469215 + 0.324445790461322i
67        0.287297708288495 + 0.327316590471862i
68        0.288672793174997 + 0.330522590222016i
69        0.285286312738354 + 0.334297055803381i
70        0.284786735119206 + 0.337486330606448i
71        0.283261587169566 + 0.340898371593846i
72        0.281855144174891 + 0.344137158185844i
73        0.279275048948165 + 0.348103895289608i
74        0.276953429939229 + 0.352077208732792i
75        0.275146050222064 + 0.355610731512224i
76        0.272093392023027 + 0.359756362104508i
77        0.270321033989090 + 0.363707794177817i
78        0.266251284881468 + 0.368463101281507i
79        0.263749790572167 + 0.373039930114154i
80        0.262281167386285 + 0.378372079435673i
81        0.257992150155298 + 0.383360901355126i
82        0.256342865829986 + 0.388786068105767i
83        0.251932446613127 + 0.394282215227473i
84        0.249038565513940 + 0.399980561938111i
85        0.245536966637245 + 0.406751453905034i
86        0.242040781305733 + 0.412383941089901i
87        0.240236587177524 + 0.419165494060697i
88        0.236725517958843 + 0.426049082895119i
89        0.234229828347182 + 0.433097693315253i
90        0.230886649401522 + 0.440396797687712i
91        0.229340490505599 + 0.447316846362464i
92        0.226732352666274 + 0.455291418348478i
93        0.224222721708530 + 0.463288173876788i
94        0.222933667238278 + 0.471379622262633i
95        0.219876444836137 + 0.479137141406710i
96        0.219249888496767 + 0.487314509598339i
97        0.217154125258385 + 0.496123219206257i
98        0.216474406231333 + 0.504631452945457i
99        0.214987970032758 + 0.512540347392450i
100        0.215601625119156 + 0.521472989180999i
101        0.215739974837045 + 0.530346471042275i
102        0.214628763199637 + 0.538884158539275i
103        0.214895958798120 + 0.547066734761426i
104        0.216132012172760 + 0.555456223114273i
105        0.215838804994722 + 0.563927903621178i
106        0.216645575593618 + 0.571894151863945i
107        0.217653040899196 + 0.579597438262141i
108        0.219726401068226 + 0.588179286484852i
109        0.219253244812545 + 0.596355605323546i
110        0.221933074908222 + 0.604119809404095i
111        0.222940549904489 + 0.611602065494655i
112        0.225029909689336 + 0.619178610510521i
113        0.225266940012791 + 0.626614481150402i
114        0.227512925306853 + 0.633558465767751i
115        0.227713880793591 + 0.640462452883964i
116        0.230137547963400 + 0.647559868377860i
117        0.231491530862348 + 0.654186518079246i
118        0.234540652019009 + 0.661126461715605i
119        0.235361363960209 + 0.667833750559685i
120        0.234636020582275 + 0.674324994324235i
121        0.237101073029435 + 0.681254960930132i
122        0.237326478871147 + 0.687608247749951i
123        0.238140486296109 + 0.694414323516892i
124        0.238795997891099 + 0.701358995796543i
125        0.241473088980627 + 0.707913104541403i
126        0.242789669379477 + 0.714902631067352i
127        0.245040997945817 + 0.722222363912149i
128        0.243478752370918 + 0.728703014598619i
129        0.245506106405337 + 0.735335567657973i
130        0.247818240262639 + 0.742612819046929i
131        0.248515954462911 + 0.749746929454615i
132        0.249466148527073 + 0.756437169939372i
133        0.251607893014301 + 0.763609496951999i
134        0.251093367841386 + 0.770362197322573i
135        0.252969234579990 + 0.777769822523277i
136        0.256114169676954 + 0.784458598707879i
137        0.256322428807093 + 0.791874247446914i
138        0.257206939514029 + 0.799015935876225i
139        0.259827822327856 + 0.806380424267288i
140        0.260343170446372 + 0.813275530218118i
141        0.262486979832761 + 0.819695579015996i
142        0.263090102309425 + 0.827243572543503i
143        0.265140910441813 + 0.834259386145045i
144        0.266016640969459 + 0.841221815077815i
145        0.267762998565223 + 0.848413320476126i
146        0.269096303762242 + 0.855552512071508i
147        0.270170203071007 + 0.862557731100249i
148        0.273775092224720 + 0.870043869487210i
149        0.276785362298543 + 0.876882071477134i
150        0.277724798345201 + 0.884111232753154i
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-10-14 19:08:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
先试试不加任何约束这个复数拟合效果如何，可以接受在慢慢加约束；反之不加约束拟合效果都不好的话，是否得先检查下公式和数据是否正确。

另外：“限制条件是分母这个多项式等于0时，这个方程的根的实部全部小于0”

“方程的根”具体指谁？每一步的K？K是已知实数，求哪个根？
赞一下 回复此楼
2楼2015-10-14 21:37:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CarrieQh

银虫 (正式写手)
引用回帖:
2楼: Originally posted by dingd at 2015-10-14 21:37:13
先试试不加任何约束这个复数拟合效果如何，可以接受在慢慢加约束；反之不加约束拟合效果都不好的话，是否得先检查下公式和数据是否正确。

另外：“限制条件是分母这个多项式等于0时，这个方程的根的实部全部小于 ...

谢谢您的回复！我尝试不加约束拟合了，但是不知道是不是版本的问题，复数拟合的话，只有实部可以拟合上，而且拟合效果很好。
至于多项式的根。。上面那个式子是我自己为了处理方便写的，实际的公式就是
1+p7*w+p8*w^2+p9*w^3+p10*w^4+p11*w^5
这个根就是w，只要求得的w的实部都小于0就行
赞一下 回复此楼
3楼2015-10-15 09:28:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CarrieQh

银虫 (正式写手)
求助呀
回复此楼
4楼2015-10-23 11:18:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 CarrieQh 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程085602 305分求调剂 +3 RichLi_ 2026-03-25 3/150 2026-03-25 19:11 by qingfeng258
[考研] 考研调剂 +5 呼呼？~+123456 2026-03-24 5/250 2026-03-25 18:15 by xcjcqu
[考研] 335分 | 材料与化工专硕 | GPA 4.07 | 有科研经历 +5 cccchenso 2026-03-23 5/250 2026-03-25 18:08 by xcjcqu
[考研] 274求调剂 +7 顾九笙要谦虚 2026-03-24 7/350 2026-03-25 15:18 by baoball
[考研] 材料学硕333求调剂 +4 北道巷 2026-03-24 4/200 2026-03-25 14:16 by mapenggao
[考研] 求调剂 +3 李李不服输 2026-03-25 3/150 2026-03-25 13:03 by cmz0325
[考研] 303求调剂 +6 元夕元 2026-03-20 7/350 2026-03-25 12:00 by edmund7
[考研] 086003食品工程求调剂 +6 淼淼111 2026-03-24 6/300 2026-03-25 10:29 by 3Strings
[考研] 292求调剂 +4 鹅鹅鹅额额额额� 2026-03-24 4/200 2026-03-24 16:41 by peike
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 305分求调剂（食品工程） +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 一志愿国科过程所081700，274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 一志愿陕师大生物学071000，298分，求调剂 +3 SYA！ 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 一志愿070300浙大化学358分，求调剂！ +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 285求调剂 +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by 星空星月
[考研] 材料与化工（0856）304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 261求B区调剂，科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
信息提示
关闭
请填处理意见
关闭