當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

yan_fei

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1406.5
帖子: 81
在線: 218.6小時(shí)
蟲號: 896603
注冊: 2009-11-07
專業(yè): 交通工程

[求助] 各位朋友，如何用matlab編寫這個(gè)方程，并求解？

各位朋友，最近遇到這么一個(gè)方程，不知道如何用matlab解決，方程如下：
公式中的未知數(shù)是t，其他的都是已知的。怎么用matlab編寫這個(gè)方程，并求解？請各位朋友幫忙，謝謝！

方程

回復(fù)此樓

1樓 2012-07-30 11:06:37

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+3, 謝謝應(yīng)助！ 2012-07-30 21:50:21
yan_fei: 金幣+130, 感謝這位朋友的回復(fù)，金幣全部給你！不過不知道為什么初值不同，結(jié)果差別很大，我想只改變alph的值求相應(yīng)的t，隨便一個(gè)初值都能求出準(zhǔn)確的解，如知道望告知！感謝您的幫助~~~ 2012-07-30 22:39:56

CODE:

function yan

global s phi lam alpha

s=3.2;

phi=1;

lam=3;

alpha=0.6;

[X,FVAL,EXITFLAG,OUTPUT]=fsolve(@sumt,10)

function y=sumt(t)

global s phi lam alpha

y=0;

for N=0:floor(s*t-phi)

    y=y+exp(-lam*t)*(lam*t)^N/gamma(N+1);

end

y=y-alpha;

Equation solved.

fsolve completed because the vector of function values is near zero
as measured by the default value of the function tolerance, and
the problem appears regular as measured by the gradient.

X =

10.0881

FVAL =

-9.6716e-010

EXITFLAG =

   1

OUTPUT =

   iterations: 2
      funcCount: 6
      algorithm: 'trust-region dogleg'
firstorderopt: 2.0551e-010
      message: [1x695 char]

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

3樓2012-07-30 21:09:48

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

查看全部 16 個(gè)回答

yan_fei

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1406.5
帖子: 81
在線: 218.6小時(shí)
蟲號: 896603
注冊: 2009-11-07
專業(yè): 交通工程

沒有朋友回復(fù)啊，是金幣少嗎？小弟金幣有限，是每天辛苦簽到積攢的，再追加100個(gè)金幣，請知道的朋友幫忙，怎么用matlab編程把這個(gè)方程解出來，公式求和符號上面的【】是取整數(shù)的意思。

贊一下

回復(fù)此樓

2樓2012-07-30 20:39:08

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

★ ★
xiegangmai: 金幣+2, 專家考核, 辛苦了！ 2012-07-30 21:50:33

或者這么寫也可以
function yan
global s phi lam alpha
s=3.2;
phi=1;
lam=3;
alpha=0.6;
[X,FVAL,EXITFLAG,OUTPUT]=fsolve(@sumt,10)

function y=sumt(t)
global s phi lam alpha
y=0;
% for N=0:floor(s*t-phi)
% y=y+exp(-lam*t)*(lam*t)^N/gamma(N+1);
% end
N=0:floor(s*t-phi);
yN=exp(-lam*t)*(lam*t).^N./gamma(N+1);
y=sum(yN);
y=y-alpha;

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

4樓2012-07-30 21:15:00

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

yan_fei

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1406.5
帖子: 81
在線: 218.6小時(shí)
蟲號: 896603
注冊: 2009-11-07
專業(yè): 交通工程

引用回帖:

4樓: Originally posted by dbb627 at 2012-07-30 21:15:00
或者這么寫也可以
function yan
global s phi lam alpha
s=3.2;
phi=1;
lam=3;
alpha=0.6;
=fsolve(@sumt,10)

function y=sumt(t)
global s phi lam alpha
y=0;
% for N=0:floor(s*t-phi)
% y=y+ ...

首先感謝您的回復(fù)，我剛試了下您的程序，好像對初始值的依賴性很強(qiáng)，沒次給不同的初始值求的結(jié)果都不一樣，而且結(jié)果幾乎和給的初始值差不多。我是想改變不同的alpha的值求出對應(yīng)的t. 程序運(yùn)行后的提示為：

No solution found.

fsolve stopped because the relative size of the current step is less than the
default value of the step size tolerance squared, but the vector of function values
is not near zero as measured by the default value of the function tolerance.

fsolve stopped because the relative norm of the current step, 4.694166e-013, is less than
max(options.TolX^2,eps) = 1.000000e-012. However, the sum of squared function values,
r = 9.759502e-003, exceeds sqrt(options.TolFun) = 1.000000e-003.

Optimization Metric                                              Options
relative norm(step) =  4.69e-013             max(TolX^2,eps) =  1e-012 (default)
r =  9.76e-003                                  sqrt(TolFun) = 1.0e-003 (default)

fsolve stopped because the relative norm of the current step, 4.694166e-013, is less than
max(options.TolX^2,eps) = 1.000000e-012. However, the sum of squared function values,
r = 9.759502e-003, exceeds sqrt(options.TolFun) = 1.000000e-003.

Optimization Metric                                              Options
relative norm(step) =  4.69e-013             max(TolX^2,eps) =  1e-012 (default)
r =  9.76e-003                                  sqrt(TolFun) = 1.0e-003 (default)

贊一下

回復(fù)此樓

5樓2012-07-30 21:53:23

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

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

[考研] 292求調(diào)劑 +9	是妍子也是研子 2026-03-30	9/450	2026-03-30 13:34 by chemdavid
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6	王一一依依 2026-03-30	8/400	2026-03-30 11:52 by yjolah
[考研] 318求調(diào)劑 +7	陳晨79 2026-03-30	7/350	2026-03-30 10:49 by 探123
[考研] 290求調(diào)劑 +3	dfffsar 2026-03-29	3/150	2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 281求調(diào)劑 +4	亞克西good 2026-03-26	6/300	2026-03-29 21:06 by 啊李999
[考研] 總分293求調(diào)劑 +8	加一一九 2026-03-25	11/550	2026-03-29 19:53 by 無際的草原
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級。 +4	XX小鄧 2026-03-29	4/200	2026-03-29 17:34 by 無際的草原
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +4	12.27 2026-03-27	4/200	2026-03-29 14:40 by zhshch
[考研] 329求調(diào)劑 +10	鈕恩雪 2026-03-25	10/500	2026-03-29 13:32 by peike
[考研] 343求調(diào)劑 +6	愛羈絆 2026-03-29	6/300	2026-03-29 12:00 by 無際的草原
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-27	6/300	2026-03-28 14:10 by 唐沐兒
[考研] 322求調(diào)劑 +5	舊吢 2026-03-24	5/250	2026-03-28 13:26 by Iveryant
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7	hanamiko 2026-03-27	7/350	2026-03-28 09:56 by 李上岸0921
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚的貓
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 求調(diào)劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325
[有機(jī)交流] 有機(jī)合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753