YunPengWang

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 192.6
帖子: 25
在線: 3.8小時(shí)
蟲號(hào): 5343104
注冊(cè): 2016-12-11
性別: GG
專業(yè): 激光

[求助] 求助已有2人參與

各位大神，怎么用matlab解這個(gè)方程(求S)

發(fā)自小木蟲Android客戶端

回復(fù)此樓

» 猜你喜歡

材料考研調(diào)劑已經(jīng)有3人回復(fù)
材料調(diào)劑已經(jīng)有12人回復(fù)
英一數(shù)一408，總分284，二戰(zhàn)真誠(chéng)求調(diào)劑已經(jīng)有14人回復(fù)
085410 一志愿211 22408分?jǐn)?shù)359求調(diào)劑已經(jīng)有4人回復(fù)
271求調(diào)劑已經(jīng)有19人回復(fù)
385分生物學(xué)（071000）求調(diào)劑已經(jīng)有3人回復(fù)
一志愿安徽大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)碩，331分求調(diào)劑已經(jīng)有3人回復(fù)
318求調(diào)劑，計(jì)算材料方向已經(jīng)有8人回復(fù)
291求調(diào)劑已經(jīng)有25人回復(fù)
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑已經(jīng)有6人回復(fù)

1樓 2016-12-12 19:01:53

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

月只藍(lán)

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

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

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

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

m文件如下：
簡(jiǎn)便起見，我用的窮舉法（肯定還有其他效率更高的方法）
參數(shù)f_number表示在區(qū)間內(nèi)計(jì)算的次數(shù)，取整數(shù)，數(shù)值越大，計(jì)算精度越高，以下代碼中取為10 000。

CODE:

function YunPengWang

clear;clc

format long

f_number=1e4;

B=2;

a=0.5;

f=0.5584;

F=0.0585E9;

s_all=linspace(4.001498e6,4.0015e6,f_number);

n=length(s_all);

for i=1:n

    s=s_all(i);

y(i)=F*f- (1+B/f/s)*log(1+f*s)/(  quadgk( @(x) exp(-x*(a^2+1))./(1+f*s*exp(-a^2*x))        ,0,Inf    ))  ;

end

plot(s_all,y,'b-',[min(s_all) max(s_all)],[0 0],'r-');

[y_min index]=min(abs(y));

residual_errory=y_min

s_ans=s_all(index)

計(jì)算結(jié)果：
方程殘差 residual_errory =5.464069545269013e-04

待求參數(shù)s的數(shù)值 s_ans = 4.001498796679668e+06

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

YunPengWang

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

2樓2016-12-12 20:38:18

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

YunPengWang

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 192.6
帖子: 25
在線: 3.8小時(shí)
蟲號(hào): 5343104
注冊(cè): 2016-12-11
性別: GG
專業(yè): 激光

送紅花一朵

引用回帖:

2樓: Originally posted by 月只藍(lán) at 2016-12-12 20:38:18
m文件如下：
簡(jiǎn)便起見，我用的窮舉法（肯定還有其他效率更高的方法）
參數(shù)f_number表示在區(qū)間內(nèi)計(jì)算的次數(shù)，取整數(shù)，數(shù)值越大，計(jì)算精度越高，以下代碼中取為10 000。

function YunPengWang
clear;clc
forma ...

如果F是一個(gè)范圍呢？就是已知F范圍(0.585e9~5.834e9)，然后畫出F與S關(guān)系圖，窮舉法是否可行？謝謝大神！

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

3樓2016-12-12 21:43:23

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

tuburfransth

金蟲 (正式寫手)

應(yīng)助: 9 (幼兒園)
金幣: 2380.1
散金: 68
紅花: 5
帖子: 384
在線: 124.8小時(shí)
蟲號(hào): 3425987
注冊(cè): 2014-09-18
性別: GG
專業(yè): 模式識(shí)別

引用回帖:

你好，看到了你的回答我想問(wèn)一下你怎么知道s_all的取值范圍是在4.001498e6,4.0015e6之間？s_all=linspace(4.001498e6,4.0015e6,f_number);

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

4樓2016-12-12 22:25:36

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

pdl9527

專家顧問(wèn) (小有名氣)

專家經(jīng)驗(yàn): +8
應(yīng)助: 100 (初中生)
金幣: 2110.3
紅花: 23
帖子: 282
在線: 112小時(shí)
蟲號(hào): 1227333
注冊(cè): 2011-03-09
性別: GG
專業(yè): 分離過(guò)程
管轄: 計(jì)算模擬

引用回帖:

4樓: Originally posted by tuburfransth at 2016-12-12 22:25:36
你好，看到了你的回答我想問(wèn)一下你怎么知道s_all的取值范圍是在4.001498e6,4.0015e6之間？s_all=linspace(4.001498e6,4.0015e6,f_number);
...

首先應(yīng)該是個(gè)較大的范圍，然后根據(jù)畫圖的結(jié)果人為進(jìn)行縮小。

贊一下

回復(fù)此樓

5樓2016-12-13 06:49:26

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

pdl9527

專家顧問(wèn) (小有名氣)

專家經(jīng)驗(yàn): +8
應(yīng)助: 100 (初中生)
金幣: 2110.3
紅花: 23
帖子: 282
在線: 112小時(shí)
蟲號(hào): 1227333
注冊(cè): 2011-03-09
性別: GG
專業(yè): 分離過(guò)程
管轄: 計(jì)算模擬

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

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

引用回帖:

3樓: Originally posted by YunPengWang at 2016-12-12 21:43:23
如果F是一個(gè)范圍呢？就是已知F范圍(0.585e9~5.834e9)，然后畫出F與S關(guān)系圖，窮舉法是否可行？謝謝大神！
...

CODE:

function question_12 

%2016-12-13 

clear;clc

B=2; 

a=0.5; 

f=0.5584; 

%在值域內(nèi)取20個(gè)F的值來(lái)求s

FF=linspace(0.585e9,5.834e9,20);

y=@(F) @(s) F*f- (1+B/f/s)*log(1+f*s)/(quadgk( @(x) exp(-x*(a^2+1))./(1+f*s*exp(-a^2*x)),0,Inf));

[s,fval,exitflag]=arrayfun(@(F) fsolve(y(F),1e8),FF);

%F和s的關(guān)系圖像

plot(FF,s,'bo--');

xlabel('F');

ylabel('s')

%最大誤差

max(abs(fval))

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

YunPengWang

6樓2016-12-13 07:23:54

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

YunPengWang

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 192.6
帖子: 25
在線: 3.8小時(shí)
蟲號(hào): 5343104
注冊(cè): 2016-12-11
性別: GG
專業(yè): 激光

送紅花一朵

引用回帖:

6樓: Originally posted by pdl9527 at 2016-12-13 07:23:54
function question_12
%2016-12-13
clear;clc
B=2;
a=0.5;
f=0.5584;
%在值域內(nèi)取20個(gè)F的值來(lái)求s
FF=linspace(0.585e9,5.834e9,20);
y=@(F) @(s) F*f- (1+B/f/s)*log(1+f*s)/(quadgk( @(x) exp(-x*(a ...

好的，我回去試試！謝謝

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

7樓2016-12-13 08:01:14

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

YunPengWang

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 192.6
帖子: 25
在線: 3.8小時(shí)
蟲號(hào): 5343104
注冊(cè): 2016-12-11
性別: GG
專業(yè): 激光

引用回帖:

你好，我還想問(wèn)一下，為什么求完S的值，回帶的時(shí)候結(jié)果會(huì)不一樣呢？

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

8樓2016-12-14 20:07:19

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

pdl9527

專家顧問(wèn) (小有名氣)

專家經(jīng)驗(yàn): +8
應(yīng)助: 100 (初中生)
金幣: 2110.3
紅花: 23
帖子: 282
在線: 112小時(shí)
蟲號(hào): 1227333
注冊(cè): 2011-03-09
性別: GG
專業(yè): 分離過(guò)程
管轄: 計(jì)算模擬

引用回帖:

8樓: Originally posted by YunPengWang at 2016-12-14 20:07:19
你好，我還想問(wèn)一下，為什么求完S的值，回帶的時(shí)候結(jié)果會(huì)不一樣呢？
...

不知道你指哪里不一樣，你將下列代碼加到程序的最末，可以看到將s(1)帶到原方程中求得的F的值是5.8500e+08，和條件中給的F范圍(0.585e9~5.834e9)中的第一個(gè)值完全吻合，同樣你可以將s(20)帶入，求得的結(jié)果和最后一個(gè)F的值是吻合的，說(shuō)明程序沒(méi)有問(wèn)題。

CODE:

yy=@(s) (1+B/f/s)*log(1+f*s)/(f*quadgk( @(x) exp(-x*(a^2+1))./(1+f*s*exp(-a^2*x)),0,Inf));

feval(yy,s(1))

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

YunPengWang

9樓2016-12-14 20:26:10

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

YunPengWang

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 192.6
帖子: 25
在線: 3.8小時(shí)
蟲號(hào): 5343104
注冊(cè): 2016-12-11
性別: GG
專業(yè): 激光

送紅花一朵

引用回帖:

9樓: Originally posted by pdl9527 at 2016-12-14 20:26:10
不知道你指哪里不一樣，你將下列代碼加到程序的最末，可以看到將s(1)帶到原方程中求得的F的值是5.8500e+08，和條件中給的F范圍(0.585e9~5.834e9)中的第一個(gè)值完全吻合，同樣你可以將s(20)帶入，求得的結(jié)果和最后一 ...

錯(cuò)誤找到了，謝謝大神�。�！

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

10樓2016-12-14 21:53:26

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

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

[考研] 301求調(diào)劑 +12	A_JiXing 2026-04-01	12/600	2026-04-01 22:47 by peike
[考研] 求調(diào)劑推薦材料 304 +21	荷包蛋hyj 2026-03-26	21/1050	2026-04-01 21:09 by lijunpoly
[考研] 食品學(xué)碩362求調(diào)劑 +3	xuanxianxian 2026-04-01	3/150	2026-04-01 21:05 by 啊李999
[考研] 085601材料工程找調(diào)劑 +20	oatmealR 2026-03-29	21/1050	2026-04-01 21:00 by lijunpoly
[考研] 0710生物學(xué)考研調(diào)劑 +3	李多米lee. 2026-03-27	4/200	2026-04-01 16:21 by zzchen2000
[考研] 調(diào)劑 +3	好好讀書。 2026-04-01	5/250	2026-04-01 15:39 by 好好讀書。
[考研] 324求調(diào)劑 +4	想上學(xué)求調(diào) 2026-04-01	5/250	2026-04-01 14:19 by ZXlzxl0425
[考研] 一志愿武理材料工程302調(diào)劑環(huán)化或化工 +15	Doleres 2026-03-31	16/800	2026-04-01 09:49 by lfj11
[考研] 289求調(diào)劑 +3	Acesczlo 2026-03-29	4/200	2026-03-31 14:48 by 熱情沙漠
[考研] 0703化學(xué) +20	妮妮ninicgb 2026-03-27	20/1000	2026-03-31 13:33 by 無(wú)際的草原
[考研] 08工科，295，接受跨專業(yè)調(diào)劑 +6	lmnlzy 2026-03-30	6/300	2026-03-31 10:04 by cal0306
[考研] 抱歉 +4	田洪有 2026-03-30	4/200	2026-03-30 21:26 by mumin1990
[考研] 085701環(huán)境工程求調(diào)劑 +11	多久上課 2026-03-27	12/600	2026-03-30 21:21 by 研究僧導(dǎo)導(dǎo)
[考研] 材料專碩 085600求調(diào)劑 +7	BBQ233 2026-03-30	7/350	2026-03-30 17:44 by oooqiao
[考研] 317分一志愿南理工材料工程本科湖工大求調(diào)劑 +12	芋泥小鈴鐺 2026-03-28	12/600	2026-03-30 17:06 by wangjy2002
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4	qianbai11 2026-03-29	4/200	2026-03-30 11:34 by 唐沐兒
[考研] 材料與化工304求B區(qū)調(diào)劑 +4	邱gl 2026-03-26	7/350	2026-03-30 08:39 by 探123
[考研] 332求92調(diào)劑 +8	蕉蕉123 2026-03-28	8/400	2026-03-29 10:46 by 周梓丹
[考研] 085405 考的11408求各位老師帶走 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 09:19 by 樂(lè)呵呵的追夢(mèng)人
[考研] 材料求調(diào)劑一志愿哈工大324 +7	閆旭東 2026-03-28	9/450	2026-03-28 08:51 by Xu de nuo