版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

津津愛睡覺

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1100.8
散金: 33
紅花: 1
帖子: 145
在線: 60小時
蟲號: 2029058
注冊: 2012-09-25
性別: GG
專業(yè): 信號理論與信號處理

[求助] matlab解方程已有2人參與

x=[0,0.25,0.5,0.75,1,1.25,1.5,1.75,2,2.25,2.5,3,3.5,4,4.5,5,5.5,6,6.5]
x=1/(15.3109*49.005^2).*[log((y./86.4241).*((86.4241-49.005)./(y-49.005)))-49.005.*(1./y-1./86.4241)]我已知x值，相求y值。。使用了fzero、fslove都各種出錯，卡在這好幾天了

。。新手，處理數(shù)據(jù)臨時學(xué)的，還請大神幫忙。。感激不盡。。希望把程序告訴我啊。。還有好多要處理。。

回復(fù)此樓

» 猜你喜歡

274求調(diào)劑已經(jīng)有7人回復(fù)
北科281學(xué)碩材料求調(diào)劑已經(jīng)有11人回復(fù)
085600 材料與化工 329分求調(diào)劑已經(jīng)有3人回復(fù)
302求調(diào)劑已經(jīng)有3人回復(fù)
312求調(diào)劑已經(jīng)有5人回復(fù)
網(wǎng)絡(luò)空間安全0839招調(diào)劑已經(jīng)有4人回復(fù)
329求調(diào)劑已經(jīng)有3人回復(fù)
【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額已經(jīng)有8人回復(fù)
085602 289分求調(diào)劑已經(jīng)有7人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有4人回復(fù)

有夢一起追

1樓 2015-07-20 09:14:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

引用回帖:

8樓: Originally posted by 津津愛睡覺 at 2015-07-20 14:09:49
想問下擬合的代碼是什么啊。。我擬合出的結(jié)果都不對。。太感謝了，你讓我覺得世界都亮了...

1stOpt軟件代碼：

CODE:

Variable  x,y ;

Function x=1/(x2*49.005^2)*(ln((y/x1)*((x1-49.005)/(y-49.005)))-49.005*(1/y-1./x1));

Data;

         0   83.9499

    0.2500   69.1806

    0.5000   63.7801

    0.7500   62.1649

    1.0000   60.6538

    1.2500   59.8982

    1.5000   58.4652

    1.7500   57.7097

    2.0000   56.9541

    2.5000   55.7696

    3.0000   54.9544

    3.5000   53.8983

    4.0000   53.1428

    4.5000   52.1404

    5.0000   51.6460

    5.5000   51.3437

    6.0000   51.1179

    6.5000   50.8551

計算結(jié)果：
相關(guān)系數(shù)(R): 0.994536826184908
相關(guān)系數(shù)之平方(R^2): 0.98910349863795
決定系數(shù)(DC): 0.98910349863795
卡方系數(shù)(Chi-Square): -0.475114996414728
F統(tǒng)計(F-Statistic): 1452.36121690615

參數(shù) 最佳估算
---------- -------------
x2 0.000124355240207823
x1 71.0613092892056

1stopt軟件獲得鏈接：http://www.gaoyang168.com/bbs/viewthread.php?tid=4586018

如果你想用MATLAB Origin等軟件，具體代碼或者操作步驟，可參見：
http://www.gaoyang168.com/bbs/viewthread.php?tid=6572963

附圖1.png

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

9樓2015-07-20 14:23:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
ben_ladeng: 金幣+2, 辛苦了 2015-07-20 12:52:08

用1stOpt很方便：

CODE:

LoopConstant x=[0,0.25,0.5,0.75,1,1.25,1.5,1.75,2,2.25,2.5,3,3.5,4,4.5,5,5.5,6,6.5];

Function

x=1/(15.3109*49.005^2)*(ln((y/86.4241)*((86.4241-49.005)/(y-49.005)))-49.005*(1/y-1/86.4241));

x y
0 -52.980531685766
0.25 -0.005325680664864
0.5 -0.00266410044168563
0.75 -0.0017763656641459
1 -0.00133239055348768
1.25 -0.00106596977280202
1.5 -0.000888340653964983
1.75 -0.000761455084693858
2 -0.000666286667945962
2.25 -0.000592264240768995
2.5 -0.000533044674782525
3 -0.000444212582842599
3.5 -0.000380759044626785
4 -0.000333167756019921
4.5 -0.000296151628028883
5 -0.00026653829158802
5.5 -0.000242308908315198
6 -0.000222117554561104
6.5 -0.000205032419060659

贊一下

回復(fù)此樓

2樓2015-07-20 10:55:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
津津愛睡覺: 金幣+100 2015-07-20 11:48:52

CODE:

function solve_eqs

clear all;clc

format long

global x

y0=-0.1;

x_data=[0,0.25,0.5,0.75,1,1.25,1.5,1.75,2,2.25,2.5,3,3.5,4,4.5,5,5.5,6,6.5];

for i=1:length(x_data)

    x=x_data(i);

[y(i) fval(i)]=fsolve(@eqs,y0);

end

[x_data ;y; fval]'

計算結(jié)果：
               x                               y                   方程殘差
               0             -1.506414574649781 0.000781864700102
0.250000000000000  -0.005325680664661 0.000000000000122
0.500000000000000  -0.002664100440992 0.000000000130176
0.750000000000000  -0.001776365663977 0.000000000043270
1.000000000000000  -0.001332390553487 0.000000000000432
1.250000000000000  -0.001065969772692 0.000000000128479
1.500000000000000  -0.000888340653914 0.000000000029346
1.750000000000000  -0.000761455084694 0.000000000000222
2.000000000000000  -0.000666286667945 0.000000000003717
2.250000000000000  -0.000592264240769 0.000000000000011
2.500000000000000  -0.000533044674760 0.000000000010490
3.000000000000000  -0.000444212582825 0.000000000004250
3.500000000000000  -0.000380759044612 0.000000000000032
4.000000000000000  -0.000333167756007 0.000000000004873
4.500000000000000  -0.000296151628018 0.000000000000042
5.000000000000000  -0.000266538291578 0.000000000001760
5.500000000000000  -0.000242308908306 0.000000000009863
6.000000000000000  -0.000222117554561 0.000000000000224
6.500000000000000  -0.000205032419053 0.000000000000051

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

3樓2015-07-20 11:42:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

津津愛睡覺

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1100.8
散金: 33
紅花: 1
帖子: 145
在線: 60小時
蟲號: 2029058
注冊: 2012-09-25
性別: GG
專業(yè): 信號理論與信號處理

引用回帖:

2樓: Originally posted by dingd at 2015-07-20 10:55:11
用1stOpt很方便：

LoopConstant x=;
Function
x=1/(15.3109*49.005^2)*(ln((y/86.4241)*((86.4241-49.005)/(y-49.005)))-49.005*(1/y-1/86.4241));

x y
0 -52.980531685766
0.25 -0.005325680664864
0 ...

哇哇

。。怎么都是負(fù)的啊。。說明我參數(shù)擬合錯了。。我原方程是x=1/(x(2).*49.005^2).*[log((y./x(1)).*((x(1)-49.005)./(y-49.005))-49.005.*(1./y-1./x(1))]..
xdata=[0,0.25,0.5,0.75,1,1.25,1.5,1.75,2,2.5,3,3.5,4,4.5,5,5.5,6,6.5]
ydata=[83.94986152,69.18056238,63.78006779,62.16487246,60.65377495,59.89822619,58.46520206,57.7096533,56.95410455,55.769572,54.95441885,53.89832959,53.14278084,52.14041949,51.64595481,51.34373531,51.11791018,50.85514711]..我擬合出了方程中的參數(shù)x1,x2然后根據(jù)xdata求y值

贊一下

回復(fù)此樓

有夢一起追

4樓2015-07-20 11:45:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

津津愛睡覺

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1100.8
散金: 33
紅花: 1
帖子: 145
在線: 60小時
蟲號: 2029058
注冊: 2012-09-25
性別: GG
專業(yè): 信號理論與信號處理

引用回帖:

3樓: Originally posted by 月只藍(lán) at 2015-07-20 11:42:53
function solve_eqs
clear all;clc
format long

global x
y0=-0.1;
x_data=;
for i=1:length(x_data)
x=x_data(i);
=fsolve(@eqs,y0);
end

'

計算結(jié)果：
x ...

哇哇。。怎么都是負(fù)的啊。。說明我參數(shù)擬合錯了。。我原方程是x=1/(x(2).*49.005^2).*[log((y./x(1)).*((x(1)-49.005)./(y-49.005))-49.005.*(1./y-1./x(1))]..
xdata=[0,0.25,0.5,0.75,1,1.25,1.5,1.75,2,2.5,3,3.5,4,4.5,5,5.5,6,6.5]
ydata=[83.94986152,69.18056238,63.78006779,62.16487246,60.65377495,59.89822619,58.46520206,57.7096533,56.95410455,55.769572,54.95441885,53.89832959,53.14278084,52.14041949,51.64595481,51.34373531,51.11791018,50.85514711]..我擬合出了方程中的參數(shù)x1,x2然后根據(jù)xdata求y值。。能不能再幫幫我啊。。感謝。。急求

贊一下

回復(fù)此樓

有夢一起追

5樓2015-07-20 11:48:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

引用回帖:

5樓: Originally posted by 津津愛睡覺 at 2015-07-20 11:48:09
哇哇。。怎么都是負(fù)的啊。。說明我參數(shù)擬合錯了。。我原方程是x=1/(x(2).*49.005^2).*..
xdata=
ydata=..我擬合出了方程中的參數(shù)x1,x2然后根據(jù)xdata求y值。。能不能再幫幫我啊。。感謝。。急求...

不好意思，3樓代碼沒給出完全，完整的如下：

CODE:

function solve_eqs

clear all;clc

format long

global x

y0=-0.1;

x_data=[0,0.25,0.5,0.75,1,1.25,1.5,1.75,2,2.25,2.5,3,3.5,4,4.5,5,5.5,6,6.5];

for i=1:length(x_data)

    x=x_data(i);

[y(i) fval(i)]=fsolve(@eqs,y0);

end

[x_data ;y; fval]'

function f=eqs(y)

global x

f=1/(15.3109*49.005^2).*(log((y./86.4241).*((86.4241-49.005)./(y-49.005)))-49.005.*(1./y-1./86.4241))-x;

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

6樓2015-07-20 12:30:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

引用回帖:

4樓: Originally posted by 津津愛睡覺 at 2015-07-20 11:45:48
哇哇

。。怎么都是負(fù)的啊。。說明我參數(shù)擬合錯了。。我原方程是x=1/(x(2).*49.005^2).*..
xdata=
ydata=..我擬合出了方程中的參數(shù)x1,x2然后根據(jù)xdata求y值...

先擬合出：
x2 = 0.000124354608284474
x1 = 71.0611809557734

相關(guān)系數(shù)(R): 0.99453642822085
相關(guān)系數(shù)之平方(R^2): 0.989102707058286
決定系數(shù)(DC): 0.989102707058287

代入之后用MATLAB解方程：

CODE:

function solve_eqs

clear all;clc

format long

global x

y0=50;

x_data=[0,0.25,0.5,0.75,1,1.25,1.5,1.75,2,2.25,2.5,3,3.5,4,4.5,5,5.5,6,6.5];

for i=1:length(x_data)

    x=x_data(i);

[y(i) fval(i)]=fsolve(@eqs,y0);

end

[x_data ;y; fval]'

function f=eqs(y)

global x

f=1/(0.000124354608284474*49.005^2).*(log((y./71.0611809557734).*((71.0611809557734-49.005)./(y-49.005)))-49.005.*(1./y-1./71.0611809557734))-x;

計算結(jié)果：
x                                              y                      方程殘差
               0             71.061023088663632 0.000011398222637
0.250000000000000  67.988292938180209 0.000000000096047
0.500000000000000  65.534689215008385 0.000000331981567
0.750000000000000  63.529133478234002 0.000000000137998
1.000000000000000  61.859615079896656 0.000000036989439
1.250000000000000  60.449405535313012 0.000003469065679
1.500000000000000  59.244034566066560 0.000000108193081
1.750000000000000  58.203506004930126 0.000000000309269
2.000000000000000  57.297823894942496 0.000000002283940
2.250000000000000  56.503949371087856 0.000000000099167
2.500000000000000  55.803867449155170 0.000000000001618
3.000000000000000  54.630610121986017 0.000000002101571
3.500000000000000  53.693164548810202 0.000000017349405
4.000000000000000  52.934038863155749 0.000000762365175
4.500000000000000  52.312777591278476 0.000000000526119
5.000000000000000  51.800007824280115 0.000000000001291
5.500000000000000  51.373858765788839 0.000000363302177
6.000000000000000  51.017693708877445 0.000000025319006
6.500000000000000  50.718627743671938 0.000000000780623

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

7樓2015-07-20 12:42:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

津津愛睡覺

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1100.8
散金: 33
紅花: 1
帖子: 145
在線: 60小時
蟲號: 2029058
注冊: 2012-09-25
性別: GG
專業(yè): 信號理論與信號處理

引用回帖:

7樓: Originally posted by 月只藍(lán) at 2015-07-20 12:42:34
先擬合出：
x2 = 0.000124354608284474
x1 = 71.0611809557734

相關(guān)系數(shù)(R): 0.99453642822085
相關(guān)系數(shù)之平方(R^2): 0.989102707058286
決定系數(shù)(DC): 0.989102707058287

代入之后用MA ...

想問下擬合的代碼是什么啊。。我擬合出的結(jié)果都不對。。太感謝了，你讓我覺得世界都亮了

贊一下

回復(fù)此樓

有夢一起追

8樓2015-07-20 14:09:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

津津愛睡覺

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1100.8
散金: 33
紅花: 1
帖子: 145
在線: 60小時
蟲號: 2029058
注冊: 2012-09-25
性別: GG
專業(yè): 信號理論與信號處理

引用回帖:

9樓: Originally posted by 月只藍(lán) at 2015-07-20 14:23:09
1stOpt軟件代碼：
Variable  x,y ;
Function x=1/(x2*49.005^2)*(ln((y/x1)*((x1-49.005)/(y-49.005)))-49.005*(1/y-1./x1));
Data;
      0 83.9499
0.2500 69.1806
0.5000 63.7801
   ...

真的很謝謝你。。授之以魚不如授之以漁。。版主兩樣都告訴了我。。版主棒棒噠！

贊一下

回復(fù)此樓

有夢一起追

10樓2015-07-20 14:34:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主津津愛睡覺的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 材料學(xué)碩333求調(diào)劑 +4	北道巷 2026-03-24	4/200	2026-03-25 14:16 by mapenggao
[考研] 347求調(diào)劑 +4	L when 2026-03-25	4/200	2026-03-25 13:37 by cocolv
[考研] 085600材料與化工調(diào)劑 +9	A-哆啦Z夢 2026-03-23	15/750	2026-03-25 11:18 by Ainin_
[考研] 278求調(diào)劑 +4	我可以上岸的對?/a> 2026-03-19	4/200	2026-03-25 11:01 by userper
[考研] 292求調(diào)劑 +4	鵝鵝鵝額額額額?/a> 2026-03-25	4/200	2026-03-25 10:37 by barlinike
[考研] 289求調(diào)劑 +10	碩星赴 2026-03-23	10/500	2026-03-25 10:11 by userper
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +4	13831862839 2026-03-24	5/250	2026-03-25 09:47 by 無際的草原
[考研] 07化學(xué)280分求調(diào)劑 +7	722865 2026-03-23	7/350	2026-03-25 09:29 by aa331100
[考研] 材料專碩找調(diào)劑 +5	哈哈哈吼吼吼哈 2026-03-23	5/250	2026-03-24 19:07 by 了了了了。。
[考研] 材料考研調(diào)劑生 +3	黃粱一夢千年 2026-03-24	3/150	2026-03-24 17:00 by barlinike
[考研] 279分求調(diào)劑一志愿211 +18	chaojifeixia 2026-03-19	20/1000	2026-03-24 10:34 by dolphin_ycj
[考研] 291求調(diào)劑 +8	hhhhxn.. 2026-03-23	8/400	2026-03-23 23:15 by peike
[考研] 306求調(diào)劑 +5	來好運(yùn)來來來 2026-03-22	5/250	2026-03-22 16:17 by BruceLiu320
[考研] 319求調(diào)劑 +4	小力氣珂珂 2026-03-20	4/200	2026-03-22 15:53 by ColorlessPI
[考研] 311求調(diào)劑 +3	26研0 2026-03-20	3/150	2026-03-22 14:46 by ColorlessPI
[考研] 材料學(xué)碩301分求調(diào)劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5	Charlieyq 2026-03-19	5/250	2026-03-20 20:35 by JourneyLucky

24小時熱門版塊排行榜

[求助] matlab解方程 已有2人參與

» 猜你喜歡

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] matlab解方程已有2人參與