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

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

24小時(shí)熱門版塊排行榜

返回列表

yan_fei

銀蟲 (小有名氣)

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

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

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

方程

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

學(xué)習(xí)知識(shí)

» 猜你喜歡

0703化學(xué)調(diào)劑，求導(dǎo)師收已經(jīng)有3人回復(fù)
生物學(xué)學(xué)碩求調(diào)劑已經(jīng)有6人回復(fù)
305分求調(diào)劑（食品工程）已經(jīng)有7人回復(fù)
材料與化工考研調(diào)劑已經(jīng)有5人回復(fù)
299求調(diào)劑已經(jīng)有6人回復(fù)
材料292調(diào)劑已經(jīng)有5人回復(fù)
085404電子信息284分求調(diào)劑已經(jīng)有4人回復(fù)
收08調(diào)劑生已經(jīng)有3人回復(fù)
081700 調(diào)劑 267分已經(jīng)有9人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)

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

matlab求解方程中的參數(shù) 已經(jīng)有21人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
高手進(jìn)來用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
求matlab求解一個(gè)方程組的問題已經(jīng)有6人回復(fù)
求各位前輩幫忙！matlab解復(fù)數(shù)方程，出現(xiàn)圖中的警告，怎來修改呢？謝謝了啊已經(jīng)有5人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問題已經(jīng)有21人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
【求助】如何用matlab求解這個(gè)微分方程？【已解決】已經(jīng)有3人回復(fù)
【資源】用matlab求解超越方程已經(jīng)有13人回復(fù)

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

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

yan_fei

銀蟲 (小有名氣)

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

沒有朋友回復(fù)啊，是金幣少嗎？小弟金幣有限，是每天辛苦簽到積攢的，再追加100個(gè)金幣，請(qǐng)知道的朋友幫忙，怎么用matlab編程把這個(gè)方程解出來，公式求和符號(hào)上面的【】是取整數(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í)
蟲號(hào): 149791
注冊(cè): 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的回帖

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í)
蟲號(hào): 149791
注冊(cè): 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í)
蟲號(hào): 896603
注冊(cè): 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ù)，我剛試了下您的程序，好像對(duì)初始值的依賴性很強(qiáng)，沒次給不同的初始值求的結(jié)果都不一樣，而且結(jié)果幾乎和給的初始值差不多。我是想改變不同的alpha的值求出對(duì)應(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的回帖

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í)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

★ ★
xiegangmai: 金幣+2, 專家考核, 辛苦了！ 2012-08-02 23:03:57

能給下參數(shù)嗎，有可能alpha值不合適，你可以看下這個(gè)例子隨t的改變 alpha是收斂的
function yan
global s phi lam alpha
s=3.2;
phi=1;
lam=3;
alpha=0.6;
%[X,FVAL,EXITFLAG,OUTPUT]=fsolve(@sumt,10)
t=0.1:1:10000;
for i=1:length(t)
ya(i)=sumt(t(i));
end
plot(t,ya,'*-')

function y=sumt(t)
global s phi lam alpha
y=0;
N=0:floor(s*t-phi);
yN= poisspdf(N,lam*t);
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.

6樓2012-07-30 22:59:15

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

yan_fei

銀蟲 (小有名氣)

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

引用回帖:

6樓: Originally posted by dbb627 at 2012-07-30 22:59:15
能給下參數(shù)嗎，有可能alpha值不合適，你可以看下這個(gè)例子隨t的改變 alpha是收斂的
function yan
global s phi lam alpha
s=3.2;
phi=1;
lam=3;
alpha=0.6;
%=fsolve(@sumt,10)
t=0.1:1:10000;
for i=1:len ...

原始數(shù)據(jù)是這樣的：phi=9;lam=0.1389;alpha=0,7%,21%,38%,54%,68%,78%,86%91%,94%,97%;
phi和lam值不變，alpha依次取0~97%對(duì)應(yīng)計(jì)算出t的值是：t<20,t=21,t=23;t=25;t=27;t=29;t=31;t=33;t=35;t=37;t=39.

贊一下

回復(fù)此樓

7樓2012-07-31 08:27:41

已閱回復(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í)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

引用回帖:

7樓: Originally posted by yan_fei at 2012-07-31 08:27:41
原始數(shù)據(jù)是這樣的：phi=9;lam=0.1389;alpha=0,7%,21%,38%,54%,68%,78%,86%91%,94%,97%;
phi和lam值不變，alpha依次取0~97%對(duì)應(yīng)計(jì)算出t的值是：t<20,t=21,t=23;t=25;t=27;t=29;t=31;t=33;t=35;t=37;t=39....

s的取值是多少呢？，如果alpha依次取0~97%對(duì)應(yīng)計(jì)算出t的值是：t<20,t=21,t=23;t=25;t=27;t=29;t=31;t=33;t=35;t=37;t=39，那么計(jì)算結(jié)果應(yīng)該不會(huì)解不出來

贊一下

回復(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.

8樓2012-07-31 08:52:14

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

yan_fei

銀蟲 (小有名氣)

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

引用回帖:

8樓: Originally posted by dbb627 at 2012-07-31 08:52:14
s的取值是多少呢？，如果alpha依次取0~97%對(duì)應(yīng)計(jì)算出t的值是：t<20,t=21,t=23;t=25;t=27;t=29;t=31;t=33;t=35;t=37;t=39，那么計(jì)算結(jié)果應(yīng)該不會(huì)解不出來...

不好意思忘了寫s的值了，s=0.5.

贊一下

回復(fù)此樓

9樓2012-07-31 08:55: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í)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

★ ★
xiegangmai: 金幣+2, 謝謝參與！ 2012-08-02 23:04:08

這個(gè)確實(shí)要估計(jì)下t，可以做圖估計(jì)
function yan
global s phi lam alpha
s=0.5;
phi=9;
lam=0.1389;
A=[0,0.07,0.21,0.38,0.54,0.68,0.78,0.86 0.91,0.94,0.97];
for i=1:11
alpha=A(i) ;
[X(i),FVAL(i),EXITFLAG,OUTPUT]=fsolve(@sumt,16+i*1.9);
end
X
FVAL
t=0.1:0.1:50;
for i=1:length(t)
ya(i)=sumt(t(i))+alpha;
end
plot(t,ya,'-')
hold on;plot(X,A,'ro')

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=poisspdf(N,lam*t);
y=sum(yN)-alpha;
%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.

10樓2012-07-31 09:27:24

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

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

返回列表

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

[考研] 086003食品工程求調(diào)劑 +4	淼淼111 2026-03-24	4/200	2026-03-24 11:53 by 544594351
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4	Higraduate 2026-03-22	7/350	2026-03-24 11:23 by 種大樹
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +12	吃吃吃才有意義 2026-03-19	13/650	2026-03-24 10:06 by NCUT_Yu
[考研] 一志愿南航材料專317分求調(diào)劑 +4	炸呀炸呀炸薯?xiàng)l 2026-03-23	4/200	2026-03-24 07:32 by wangy0907
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8	皮皮皮123456 2026-03-22	8/400	2026-03-23 23:47 by 熱情沙漠
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑 +7	vv迷 2026-03-22	7/350	2026-03-23 23:44 by Txy@872106
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	6/300	2026-03-23 21:45 by 子系博
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt
[考研] 323求調(diào)劑 +6	洼小桶 2026-03-18	6/300	2026-03-23 00:29 by king123！
[考研] 0854電子信息求調(diào)劑 +3	α____ 2026-03-22	3/150	2026-03-22 21:28 by zhq0425
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	4/200	2026-03-22 15:12 by 朱芷琳
[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無際的草原
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo

24小時(shí)熱門版塊排行榜

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

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

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

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

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

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