YunPengWang

新蟲 (初入文壇)

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

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

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

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

回復此樓

» 猜你喜歡

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

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

月只藍

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

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

【答案】應助回帖

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

m文件如下：
簡便起見，我用的窮舉法（肯定還有其他效率更高的方法）
參數(shù)f_number表示在區(qū)間內計算的次數(shù)，取整數(shù)，數(shù)值越大，計算精度越高，以下代碼中取為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)

計算結果：
方程殘差 residual_errory =5.464069545269013e-04

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

贊一下

回復此樓

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

YunPengWang

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

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

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

YunPengWang

新蟲 (初入文壇)

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

送紅花一朵

引用回帖:

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

function YunPengWang
clear;clc
forma ...

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

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

贊一下

回復此樓

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

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

tuburfransth

金蟲 (正式寫手)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

pdl9527

專家顧問 (小有名氣)

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

引用回帖:

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

首先應該是個較大的范圍，然后根據畫圖的結果人為進行縮小。

贊一下

回復此樓

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

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

pdl9527

專家顧問 (小有名氣)

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

【答案】應助回帖

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

引用回帖:

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

CODE:

function question_12 

%2016-12-13 

clear;clc

B=2; 

a=0.5; 

f=0.5584; 

%在值域內取20個F的值來求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的關系圖像

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

xlabel('F');

ylabel('s')

%最大誤差

max(abs(fval))

贊一下

回復此樓

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

YunPengWang

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

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

YunPengWang

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 192.6
帖子: 25
在線: 3.8小時
蟲號: 5343104
注冊: 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;
%在值域內取20個F的值來求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客戶端

贊一下

回復此樓

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

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

YunPengWang

新蟲 (初入文壇)

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

引用回帖:

你好，我還想問一下，為什么求完S的值，回帶的時候結果會不一樣呢？

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

贊一下

回復此樓

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

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

pdl9527

專家顧問 (小有名氣)

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

引用回帖:

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

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

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))

贊一下

回復此樓

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

YunPengWang

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

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

YunPengWang

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 192.6
帖子: 25
在線: 3.8小時
蟲號: 5343104
注冊: 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)中的第一個值完全吻合，同樣你可以將s(20)帶入，求得的結果和最后一 ...

錯誤找到了，謝謝大神�。。�

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

贊一下

回復此樓

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

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

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

[考研] 336求調劑 +5	kiyy 2026-04-01	5/250	2026-04-02 07:14 by chixmc
[考研] 一志愿上海海洋大學083200食品學碩，求調劑，接受其他專業(yè) +5	what張 2026-04-01	5/250	2026-04-01 22:35 by guanxin1001
[考研] 085602化學工程268分蹲調劑 +8	月照花林。 2026-04-01	8/400	2026-04-01 22:08 by 無際的草原
[考研] 求調劑推薦材料 304 +21	荷包蛋hyj 2026-03-26	21/1050	2026-04-01 21:09 by lijunpoly
[考研] 310分求調劑 +4	成功上岸wang 2026-04-01	4/200	2026-04-01 20:35 by liu823948201
[考研] 332求調劑 +8	Lyy930824@ 2026-03-29	8/400	2026-04-01 18:40 by 千鈞澄玉宇啊
[考研] 353求調劑 +4	拉鉤不許變 2026-04-01	4/200	2026-04-01 18:10 by 記事本2026
[考研] 283求調劑 +9	A child 2026-03-28	9/450	2026-04-01 14:20 by Jaylen.
[考研] 08工科，295，接受跨專業(yè)調劑 +6	lmnlzy 2026-03-31	6/300	2026-04-01 11:02 by 逆水乘風
[考研] 318求調劑 +3	篤行致遠. 2026-03-31	3/150	2026-03-31 20:27 by 求調劑zz
[考研] 一志愿北京科技大學材料學碩328分求調劑 +4	1段時間 2026-03-31	5/250	2026-03-31 19:56 by wxiongid
[考研] 318求調劑 +10	陳晨79 2026-03-30	10/500	2026-03-31 17:37 by 544594351
[考研] 269求調劑 +4	我想讀研11 2026-03-31	4/200	2026-03-31 10:04 by cal0306
[考研] 抱歉 +4	田洪有 2026-03-30	4/200	2026-03-30 21:26 by mumin1990
[考研] 310求調劑 +10	爭取九點睡 2026-03-30	10/500	2026-03-30 16:45 by ztnimte
[考研] 086000生物與醫(yī)藥調劑 +5	Feisty。 2026-03-28	9/450	2026-03-29 12:02 by longlotian
[考研] 356求調劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 材料求調劑一志愿哈工大324 +7	閆旭東 2026-03-28	9/450	2026-03-28 08:51 by Xu de nuo
[考研] 315調劑 +4	0860求調劑 2026-03-26	5/250	2026-03-27 11:23 by wangjy2002
[考研] 341求調劑 +7	青檸檬1 2026-03-26	7/350	2026-03-27 00:19 by wxiongid