 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 744 | 回復: 1 | ||
想飛就飛1985捐助貴賓 (初入文壇)
|
[求助]
求解一個含約束的非線性參數(shù)擬合問題，方程為冪函數(shù)，請各位大神幫忙，多謝！ 已有1人參與
|
|
y=[184,107,49,15366,575,19120,5504] x=[5.15,8.21,10.39,1.07,2.37,0.94,1.50] f(x)為指數(shù)函數(shù)：f(x)=a*x^b 含約束的非線性參數(shù)擬合的最優(yōu)目標為: min{[y1-f(x1)]^2+[y2-f(x2)]^2+[y3-f(x3)]^2……[yn-f(xn)]^2} 約束條件為： 0.5yi<f(xi)<2yi 本題的目的是想求使目標函數(shù)最優(yōu)的冪函數(shù)的參數(shù)a和b。 我初學matlab，看過一些含約束的曲線擬合的例子，都是對方差等參數(shù)進行約束，也看過一些參數(shù)擬合的例子，但是怎么添加約束沒有找到，這個問題困擾我近1個月了，請各位幫幫忙，幫我解決一下這個問題，多謝！ |
回復此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有12人回復
-
一志愿西南交大，求調(diào)劑
已經(jīng)有4人回復
-
0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287
已經(jīng)有7人回復
-
297求調(diào)劑
已經(jīng)有7人回復
-
材料專碩306英一數(shù)二
已經(jīng)有13人回復
-
0703化學調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有10人回復
-
0817 化學工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有6人回復
-
收復試調(diào)劑生
已經(jīng)有4人回復
-
材料專碩274一志愿陜西師范大學求調(diào)劑
已經(jīng)有6人回復
-
331求調(diào)劑（0703有機化學
已經(jīng)有8人回復
文刀星
木蟲 (小有名氣)
- 應助: 13 (小學生)
- 金幣: 3651.2
- 散金: 500
- 紅花: 5
- 帖子: 60
- 在線: 100.9小時
- 蟲號: 3693767
- 注冊: 2015-02-22
- 性別: GG
- 專業(yè): 導航、制導與傳感技術
【答案】應助回帖
感謝參與，應助指數(shù) +1
|
擬合不太懂，不過這個問題可以用優(yōu)化算法來算。只是你的數(shù)據(jù)太奇怪，找不到滿足約束的a，b。程序如下，僅作參考。 A=[]; b=[]; Aeq=[]; beq=[]; lb=[0,0]; ub=[1000,1000]; options = gaoptimset('TolFun',1e-8,'TolCon',1e-8); X0 = ga(@my_obj,2,A,b,Aeq,beq,lb,ub,@my_con,options); %X0是用遺傳算法優(yōu)化出來的結果，可能不太好，進一步用下面的算法； X = fmincon(@my_obj,X0,A,b,Aeq,beq,lb,ub,@my_con); %將遺傳算法得出的結果作為初值進行優(yōu)化可得更精確的結果； function f=my_obj(x) % 目標函數(shù) %輸入?yún)?shù)x=[a,b]; % Y=[184,107,49,15366,575,19120,5504]; % X=[5.15,8.21,10.39,1.07,2.37,0.94,1.50]; %由于你的數(shù)據(jù)找不到滿足約束的點，所以用下面的來替代，令a=1,b=2;這樣可以看出算法是有效的。 X=1:7; Y=X.^2; fx=x(1)*(X.^x(2)); f=sum((Y-fx).^2); function [c,ceq]=my_con(x) %約束條件 % Y=[184,107,49,15366,575,19120,5504]; % X=[5.15,8.21,10.39,1.07,2.37,0.94,1.50]; X=1:7; Y=X.^2; fx=x(1)*(X.^x(2)); c1=0.5*Y'-fx'; c2=fx'-2*Y'; c=[c1;c2]; ceq=[]; |
2樓2015-09-29 10:01:30
相關版塊跳轉
我要訂閱樓主 想飛就飛1985 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +7 | 戲精丹丹丹 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿985，本科211，0817化學工程與技術319求調(diào)劑 +6 | Liwangman 2026-03-15 | 6/300 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
|
[考研] 085600材料與化工 +4 | 安全上岸！ 2026-03-16 | 4/200 |
|
|
|
[考研] 東南大學364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學獎覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 材料與化工一志愿南昌大學327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 復試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
|
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
|
[考研] 26調(diào)劑/材料科學與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|