yuyi2009

銅蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 60.6
散金: 29
紅花: 4
帖子: 289
在線: 174.3小時
蟲號: 1103053
注冊: 2010-09-19
性別: GG
專業(yè): 計算機網(wǎng)絡(luò)

[求助] matlab程序問題求解答

function[c,err,yc]=bisect(f,a,b,delta)
% Input  - f is the function input as a string 'f'
%       - a and b are the left and right endpoints
%       - delta is the tolerance
% Output - c is the zero
%       - yc = f(c)
%       - err is the error estimate for c
ya = feval(f,a);
yb = feval(f,b);
if ya * yb > 0, end
max1 = 1+ round((log(b-a)-log(delta))/log(2));
for k =1:max1
c =(a+b)/2;
yc = feval(f,c);
if yc==0;
      a=c;
      b=c;
elseif yb * yc > 0
      b=c;
      yb=yc;
else
      a=c;
      ya=yc;
end
if b - a < delta, break, end
end
c=(a+b)/2;
err=abs(b-a);
yc=feval(f,c);    % 求函數(shù)在c處的函數(shù)值

for T=200:10:300
K = exp(-9.3932+3204.71/T+0.835931*log(T)...
   +2.35267*10^-3*T-1.8736*10^-6*T^2+5.1606*10^-10*T^3)
f=inline('(4*K-1)*x^2-8*K*x+4*K');
bisect(f,0,1,10^-5)
end
主程序調(diào)用函數(shù)出錯

回復(fù)此樓

» 猜你喜歡

284求調(diào)劑已經(jīng)有10人回復(fù)
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
07化學(xué)280分求調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有12人回復(fù)
求材料，環(huán)境專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
335求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
一志愿吉大化學(xué)322求調(diào)劑已經(jīng)有4人回復(fù)
環(huán)境學(xué)碩288求調(diào)劑已經(jīng)有8人回復(fù)
341求調(diào)劑(一志愿湖南大學(xué)070300) 已經(jīng)有6人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

請教MATLAB中的LMI求解問題，急！��！已經(jīng)有7人回復(fù)
求高手解答Matlab求解微分方程的初始值和邊界值問題已經(jīng)有5人回復(fù)
matlab出錯。。。求大神解答。。。已經(jīng)有3人回復(fù)
matlab稀疏矩陣問題求教已經(jīng)有3人回復(fù)
在matlab里運行SVM的問題，期待著您的解答...... 已經(jīng)有5人回復(fù)
關(guān)于matlab線性規(guī)劃背包問題的解法求救呀~~方程知道但不知道如何解答···· 已經(jīng)有4人回復(fù)
求助如何用matlab取出直線上的點已經(jīng)有10人回復(fù)
請幫忙用matlab求矩陣特征值已經(jīng)有8人回復(fù)
matlab用多元線性最小二乘法求解多元線性方程的問題已經(jīng)有11人回復(fù)
求助��！求一程序，用matlab程序做，用蒙特卡洛方法模擬已經(jīng)有10人回復(fù)
matlab最小二乘問題求助怎么定義多變量已經(jīng)有3人回復(fù)
matlab出錯求助高手解答已經(jīng)有5人回復(fù)
求matlab求解一個方程組的問題已經(jīng)有6人回復(fù)
求助matlab軟件壓縮包已經(jīng)有4人回復(fù)
【求助】求MAtlab求解程序！��！已經(jīng)有20人回復(fù)
如何利用MATLAB求最小值已經(jīng)有4人回復(fù)
【求助】求HMM的Matlab程序已經(jīng)有4人回復(fù)

1樓 2013-11-24 16:50:01

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

zhg19

金蟲 (小有名氣)

應(yīng)助: 33 (小學(xué)生)
金幣: 894.6
紅花: 3
帖子: 212
在線: 114.5小時
蟲號: 1895765
注冊: 2012-07-16
性別: GG
專業(yè): 半導(dǎo)體物理

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

感謝參與，應(yīng)助指數(shù) +1

LZ詳細描述一下問題。。。

贊一下

回復(fù)此樓

2樓2013-11-24 17:21:08

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

yuyi2009

銅蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 60.6
散金: 29
紅花: 4
帖子: 289
在線: 174.3小時
蟲號: 1103053
注冊: 2010-09-19
性別: GG
專業(yè): 計算機網(wǎng)絡(luò)

引用回帖:

2樓: Originally posted by zhg19 at 2013-11-24 17:21:08
LZ詳細描述一下問題。。。

給了一個與溫度相關(guān)的平衡方程，求200-300每隔10的平衡常數(shù)，以及平衡轉(zhuǎn)化率。。。要求用二分法求解~

贊一下

回復(fù)此樓

3樓2013-11-24 17:46:06

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

zhg19

金蟲 (小有名氣)

應(yīng)助: 33 (小學(xué)生)
金幣: 894.6
紅花: 3
帖子: 212
在線: 114.5小時
蟲號: 1895765
注冊: 2012-07-16
性別: GG
專業(yè): 半導(dǎo)體物理

引用回帖:

3樓: Originally posted by yuyi2009 at 2013-11-24 17:46:06
給了一個與溫度相關(guān)的平衡方程，求200-300每隔10的平衡常數(shù)，以及平衡轉(zhuǎn)化率。。。要求用二分法求解~...

主要是程序出現(xiàn)的問題是什么呢

贊一下

回復(fù)此樓

4樓2013-11-24 22:09:26

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

yuyi2009

銅蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 60.6
散金: 29
紅花: 4
帖子: 289
在線: 174.3小時
蟲號: 1103053
注冊: 2010-09-19
性別: GG
專業(yè): 計算機網(wǎng)絡(luò)

引用回帖:

4樓: Originally posted by zhg19 at 2013-11-24 22:09:26
主要是程序出現(xiàn)的問題是什么呢...

main

K =

26.5468

Error using inline/feval (line 23)
Not enough inputs to inline function.

Error in bisect (line 8)
ya = feval(f,a);

Error in main (line 5)
bisect(f,0,1,10^-5)

贊一下

回復(fù)此樓

5樓2013-11-24 22:21:01

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

yuyi2009

銅蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 60.6
散金: 29
紅花: 4
帖子: 289
在線: 174.3小時
蟲號: 1103053
注冊: 2010-09-19
性別: GG
專業(yè): 計算機網(wǎng)絡(luò)

引用回帖:

4樓: Originally posted by zhg19 at 2013-11-24 22:09:26
主要是程序出現(xiàn)的問題是什么呢...

就是我聲明的函數(shù)有K，被認為是變量，所以就少參數(shù)了。。。但是前面的K是可以求出來的~這個就不知道怎么弄了，它為什么不代進去啊~

贊一下

回復(fù)此樓

6樓2013-11-24 22:22:50

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

Bruce_lee

銅蟲 (著名寫手)

應(yīng)助: 5 (幼兒園)
金幣: 2632.7
散金: 424
紅花: 4
帖子: 1277
在線: 317小時
蟲號: 948103
注冊: 2010-01-25
性別: GG
專業(yè): 聚合物共混與復(fù)合材料

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
yuyi2009: 金幣+50, ★★★★★最佳答案, 再次感謝 2013-11-24 23:47:56

function problem_solve
%--------------------------------------------
disp('----------------分割線-----------------');
%--------------------------------------------
function[c,err,yc]=bisect(f,a,b,delta)
% Input  - f is the function input as a string 'f'
%       - a and b are the left and right endpoints
%       - delta is the tolerance
% Output - c is the zero
%       - yc = f(c)
%       - err is the error estimate for c
ya = feval(f,a);
yb = feval(f,b);
if ya * yb > 0, end
max1 = 1+ round((log(b-a)-log(delta))/log(2));
for k =1:max1
c =(a+b)/2;
yc = feval(f,c);
if yc==0;
      a=c;
      b=c;
elseif yb * yc > 0
      b=c;
      yb=yc;
else
      a=c;
      ya=yc;
end
if b - a < delta, break, end
end
c=(a+b)/2;
err=abs(b-a);
yc=feval(f,c);    % 求函數(shù)在c處的函數(shù)值
end
%---------------------------------------------
function f=subfun(x,K)
f=(4*K-1).*x.^2-8*K.*x+4*K;
end
%---------------------------------------------
for T=200:10:300
K = exp(-9.3932+3204.71/T+0.835931*log(T)...
   +2.35267*10^-3*T-1.8736*10^-6*T^2+5.1606*10^-10*T^3)
[c,err,yc]=bisect(@(x)subfun(x,K),0,1,10^-5)
disp('----------------分割線-----------------');
end
end

樓主，希望能幫助到你。

贊一下

回復(fù)此樓

努力就有收獲！

7樓2013-11-24 23:37:40

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

Bruce_lee

銅蟲 (著名寫手)

應(yīng)助: 5 (幼兒園)
金幣: 2632.7
散金: 424
紅花: 4
帖子: 1277
在線: 317小時
蟲號: 948103
注冊: 2010-01-25
性別: GG
專業(yè): 聚合物共混與復(fù)合材料

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

★
xzhdty: 金幣+1, 謝謝參與 2013-11-25 06:46:02

將這個代碼保存為m文件，直接在命令窗口輸入problem_solve即可。

贊一下

回復(fù)此樓

努力就有收獲！

8樓2013-11-24 23:41:39

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

yuyi2009

銅蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 60.6
散金: 29
紅花: 4
帖子: 289
在線: 174.3小時
蟲號: 1103053
注冊: 2010-09-19
性別: GG
專業(yè): 計算機網(wǎng)絡(luò)

引用回帖:

7樓: Originally posted by Bruce_lee at 2013-11-24 23:37:40
function problem_solve
%--------------------------------------------
disp('----------------分割線-----------------');
%--------------------------------------------
function=bisect(f,a,b,delta)
% ...

多謝~能順便解釋下為什么嘛~

贊一下

回復(fù)此樓

9樓2013-11-24 23:47:44

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

Bruce_lee

銅蟲 (著名寫手)

應(yīng)助: 5 (幼兒園)
金幣: 2632.7
散金: 424
紅花: 4
帖子: 1277
在線: 317小時
蟲號: 948103
注冊: 2010-01-25
性別: GG
專業(yè): 聚合物共混與復(fù)合材料

像這樣有傳遞參數(shù)的情況，一般都是將需要求解的函數(shù)定義為“function f=函數(shù)名(x,k)”的形式，這里k為傳遞參數(shù)。用inlne()函數(shù)，給出的是函數(shù)的符號形式，有時就會出錯。

贊一下

回復(fù)此樓

努力就有收獲！

10樓2013-11-25 15:34:54

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

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

[考研] 07化學(xué)280分求調(diào)劑 +4	722865 2026-03-23	4/200	2026-03-24 00:01 by chixmc
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	6/300	2026-03-23 21:45 by 子系博
[考研] 070300化學(xué)求調(diào)劑 +8	苑豆豆 2026-03-20	8/400	2026-03-23 20:57 by baobaoye
[考研] 材料專碩英一數(shù)二306 +8	z1z2z3879 2026-03-18	8/400	2026-03-23 20:49 by baobaoye
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 0854電子信息求調(diào)劑 324 +3	Promise-jyl 2026-03-23	3/150	2026-03-23 13:43 by wangkm
[考研] 北科281學(xué)碩材料求調(diào)劑 +8	tcxiaoxx 2026-03-20	9/450	2026-03-23 12:16 by tcxiaoxx
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +4	酥酥魚.. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	4/200	2026-03-22 15:12 by 朱芷琳
[考研] 化學(xué)調(diào)劑 +5	yzysaa 2026-03-21	5/250	2026-03-21 22:12 by peike
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 268求調(diào)劑 +9	簡單點0 2026-03-17	9/450	2026-03-21 15:37 by lature00
[考研] 301求調(diào)劑 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518