 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 839 | 回復(fù): 3 | ||||
zw073916木蟲(chóng) (著名寫手)
|
[求助]
求助Matlab編程，萬(wàn)分感謝 已有1人參與
|
|
有2組數(shù)： x: 0 10 40 110 155 215 310 410 505 605 700 800 900 y: 0 0.00468 0.00704 0.00996 0.01132 0.01308 0.01568 0.01868 0.02148 0.02468 0.02828 0.03292 0.03888 用方程 y’=(p(1)*sinh((p(2)*130*(1-p(5)*(1-exp(-(p(4)*y/(p(5)*130))))))/(((1+p(6)*x)^(-1/3))*(1-p(3)*y)))) 進(jìn)行擬合 求解參數(shù)P(1)-p(6)？ y'為y的導(dǎo)數(shù)，上面方程可以表示為y'=f(x,y)由于這是個(gè)隱函數(shù)，用origin擬合不了，據(jù)說(shuō)可以用matlab用最小二乘法進(jìn)行非線性擬合，可是我不會(huì)。哪位大牛能幫忙把程序?qū)懸幌拢? |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有6人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有193人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有65人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助各位大神，幫忙用MATLAB編程數(shù)出一個(gè)立方體內(nèi)所有顆粒的個(gè)數(shù)
已經(jīng)有0人回復(fù)
- 求matlab wolfe線搜索下DY共軛梯度法程序，新人求助，萬(wàn)分感謝
已經(jīng)有1人回復(fù)
- 求助matlab牛人隨機(jī)生成序列相鄰兩數(shù)不相等
已經(jīng)有2人回復(fù)
- matlab動(dòng)力學(xué)模型編程求助
已經(jīng)有6人回復(fù)
- 求助matlab編程問(wèn)題~
已經(jīng)有1人回復(fù)
- 求助高手幫忙運(yùn)用matlab分析數(shù)據(jù)，得出回歸系數(shù)
已經(jīng)有2人回復(fù)
- 求助代碼
已經(jīng)有0人回復(fù)
- 【求助】求助 matlab編程：可否實(shí)現(xiàn)清除文本文檔數(shù)據(jù)
已經(jīng)有5人回復(fù)
- matcom混合編程求助
已經(jīng)有1人回復(fù)
- 【求助】向高手求助matlab編碼
已經(jīng)有2人回復(fù)
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲(chóng)號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
zw073916: 金幣+100, ★★★★★最佳答案 2016-05-16 08:41:48
感謝參與，應(yīng)助指數(shù) +1
zw073916: 金幣+100, ★★★★★最佳答案 2016-05-16 08:41:48
| MATLAB代碼參見(jiàn)：http://www.gaoyang168.com/bbs/viewthread.php?tid=7603645 |
2樓2016-05-14 23:15:42
zw073916
木蟲(chóng) (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1475.9
- 散金: 8397
- 紅花: 1
- 帖子: 1236
- 在線: 144.2小時(shí)
- 蟲(chóng)號(hào): 3433901
- 注冊(cè): 2014-09-22
- 性別: GG
- 專業(yè): 金屬材料的力學(xué)行為
|
我按照你給我的代碼修改后，怎么計(jì)算不了 clear all;clc format long tspan=[20 50 120 165 225 320 420 515 615 710 810 910]-10; yexp=[0.00468 0.00704 0.00996 0.01132 0.01308 0.01568 0.01868 0.02148 0.02468 0.02828 0.03292 0.03888]'; k0=[1 1 1 1 1 1]; %%%請(qǐng)注意這里，初值的選取 y0=0; lb=-[1 1 1 1 1 1]*1e3; ub=[1e12 1e9 1e6 1000 10 1]; yy=[y0 yexp']; [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp); fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為:\n') fprintf('\t待擬合參數(shù) k1 = %.6f\n',k(1)) fprintf('\t待擬合參數(shù) k2 = %.6f\n',k(2)) fprintf('\t待擬合參數(shù) k3 = %.6f\n',k(3)) fprintf('\t待擬合參數(shù) k4 = %.6f\n',k(4)) fprintf('\t待擬合參數(shù) k5 = %.6f\n',k(5)) fprintf('\t待擬合參數(shù) k6 = %.6f\n',k(6)) fprintf(' \t殘差平方和= %.6f\n\n',resnorm) ts=0:1:max(tspan); [ts ys]=ode45(@KineticsEqs,ts,y0,[],k); [ttt XXsim] = ode45(@KineticsEqs,tspan,y0,[],k); y=XXsim(2:end); xexp=yexp; R2=1-sum((xexp-y).^2)./sum((xexp-mean(y)).^2); fprintf('\n\t決定系數(shù)R-Square = %.6f',R2); figure(1) plot(ts,ys,'b',tspan,yy,'or'),legend('計(jì)算值','實(shí)驗(yàn)值','Location','best'); yr=y-yexp; figure(2) plot(tspan(2:end),yr,'r*',[-100 1000],[0 0]),axis([-100 1000 -50 50]); figure(3) plot(yexp,y,'ro',[0 0.05],[0 0.05],'b-'); %--------------------------------------------------------- function f = ObjFunc(k,tspan,y0,yexp) [t Xsim] = ode45(@KineticsEqs,tspan,y0,[],k) ; ysim = Xsim(2:end); size(ysim); size(yexp); f=ysim-yexp; %---------------------------------------------------------- function dydt = KineticsEqs(t,y,k) beta(1)=k(1); beta(2)=k(2); beta(3)=k(3); beta(4)=k(4); beta(5)=k(5); beta(6)=k(6); dydt = (beta(1)*sinh((beta(2)*130*(1-beta(5)*(1-exp(-(beta(4)*y/(beta(5)*130))))))/(((1+beta(6)*t)^(-1/3))*(1-beta(3)*y)))); |
3樓2016-05-16 18:12:58
zw073916
木蟲(chóng) (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1475.9
- 散金: 8397
- 紅花: 1
- 帖子: 1236
- 在線: 144.2小時(shí)
- 蟲(chóng)號(hào): 3433901
- 注冊(cè): 2014-09-22
- 性別: GG
- 專業(yè): 金屬材料的力學(xué)行為
4樓2016-05-16 18:14:03
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 學(xué)員yMdbWV 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-17 | 8/400 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +11 | chaojifeixia 2026-03-19 | 12/600 |
|
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +5 | 77chaselx 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +5 | 陽(yáng)陽(yáng)哇塞 2026-03-14 | 10/500 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|