shuitong110

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 209.3
散金: 1042
紅花: 6
帖子: 263
在線: 85.7小時(shí)
蟲號(hào): 2772293
注冊(cè): 2013-11-02
性別: MM
專業(yè): 半導(dǎo)體微納機(jī)電器件與系統(tǒng)

[求助] 程序有誤！！MTLAB 對(duì)函數(shù)積分（變量L），積分后的函數(shù)M再進(jìn)行sin函數(shù)計(jì)算得到函數(shù)MM 已有2人參與

程序有誤��！MTLAB 對(duì)函數(shù)積分（變量L），積分后的函數(shù)M再進(jìn)行sin函數(shù)計(jì)算得到函數(shù)MM，最后繪制關(guān)于L，MM的圖。繪制的代碼如下，能出圖，但是和書上的結(jié)果差別很大，出來的圖形應(yīng)該是近似于正弦函數(shù)的圖。我運(yùn)算出來的怎么是一條直線，求大俠解答！�。�！
clear;
clear;
L=0:0.02*10^-6:70*10^-6;

d=0.1*10^-6;
R=100*10^-6;
P=0.3;
k=2*pi./(1.55*10^-6);
B=11*10^6;
n1=3.45;
n2=1;
w=0.5*10^-6;
y1=sqrt(abs(k.^2*n1.^2-B.^2));
y2=sqrt(abs(B.^2-k.^2*n2.^2));

A=2*P*y1.^2*y2.^2;
C=B*k*k*(n1.^2-n2.^2)*(2+y2*w);
D=exp(-y2*d-y2*R+y2.*sqrt(R*R-L.*L));
M=A.*D./C;

S=trapz(L,M);
MM=sin(S);
plot(L,MM);

1111.gif

回復(fù)此樓

» 猜你喜歡

2026年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有6人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有280人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有65人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)

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

【求助】函數(shù)sin(ax)*exp(-ibx)/ax在積分區(qū)間m1到m2的積分【已解決】已經(jīng)有12人回復(fù)

1樓 2014-04-25 11:01:49

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

回帖支持 ( 顯示支持度最高的前 50 名 )

alorsmail

鐵桿木蟲 (著名寫手)

應(yīng)助: 52 (初中生)
金幣: 7398.3
紅花: 41
帖子: 2104
在線: 237小時(shí)
蟲號(hào): 2537936
注冊(cè): 2013-07-09
專業(yè): 系統(tǒng)生物學(xué)

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

★
fegg7502: 金幣+1, 應(yīng)助指數(shù)+1, 鼓勵(lì)交流 2014-05-09 08:10:37

如果要得到一個(gè)曲線的話，你需要加一個(gè)循環(huán)，
示例如下，結(jié)果見圖片
==========================
L=0:0.02*10^-6:140*10^-6 ;

d=0.1*10^-6 ;
R=100*10^-6 ;
P=0.3 ;
k=2*pi./(1.55*10^-6) ;
B=11*10^6 ;
n1=3.45 ;
n2=1 ;
w=0.5*10^-6 ;
y1=sqrt(abs(k.^2*n1.^2-B.^2)) ;
y2=sqrt(abs(B.^2-k.^2*n2.^2)) ;

MM = zeros(size(L)) ;
for i = 10:length(L)
Lp = L(1:i) ;
A=2*P*y1.^2*y2.^2 ;
C=B*k*k*(n1.^2-n2.^2)*(2+y2*w) ;
D=exp(-y2*d-y2*R+y2.*sqrt(R*R-Lp.*Lp)) ;
M=A.*D./C ;
S=trapz(Lp,M) ;
MM(i)=sin(S) ;
end

plot(L,MM) ;
==========================

程序有誤��！MTLAB 對(duì)函數(shù)積分（變量L），積分后的函數(shù)M再進(jìn)行sin函數(shù)計(jì)算得到函數(shù)MM-1

figure.png

贊一下(2人)

回復(fù)此樓

6樓2014-04-25 16:12:52

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

普通回帖

alorsmail

鐵桿木蟲 (著名寫手)

應(yīng)助: 52 (初中生)
金幣: 7398.3
紅花: 41
帖子: 2104
在線: 237小時(shí)
蟲號(hào): 2537936
注冊(cè): 2013-07-09
專業(yè): 系統(tǒng)生物學(xué)

最好能提供積分的原方程，這樣方便大家對(duì)比

贊一下

回復(fù)此樓

2樓2014-04-25 12:30:48

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

shuitong110

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 209.3
散金: 1042
紅花: 6
帖子: 263
在線: 85.7小時(shí)
蟲號(hào): 2772293
注冊(cè): 2013-11-02
性別: MM
專業(yè): 半導(dǎo)體微納機(jī)電器件與系統(tǒng)

引用回帖:

2樓: Originally posted by alorsmail at 2014-04-25 12:30:48
最好能提供積分的原方程，這樣方便大家對(duì)比

已經(jīng)將函數(shù)上傳至圖片格式，希望給與指點(diǎn)。謝謝！

贊一下

回復(fù)此樓

3樓2014-04-25 14:50:45

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

shuitong110

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 209.3
散金: 1042
紅花: 6
帖子: 263
在線: 85.7小時(shí)
蟲號(hào): 2772293
注冊(cè): 2013-11-02
性別: MM
專業(yè): 半導(dǎo)體微納機(jī)電器件與系統(tǒng)

★
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-05-09 08:10:51

圖示函數(shù)

程序有誤��！MTLAB 對(duì)函數(shù)積分（變量L），積分后的函數(shù)M再進(jìn)行sin函數(shù)計(jì)算得到函數(shù)MM-2

函數(shù).jpg

贊一下

回復(fù)此樓

4樓2014-04-25 14:52:37

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

alorsmail

鐵桿木蟲 (著名寫手)

應(yīng)助: 52 (初中生)
金幣: 7398.3
紅花: 41
帖子: 2104
在線: 237小時(shí)
蟲號(hào): 2537936
注冊(cè): 2013-07-09
專業(yè): 系統(tǒng)生物學(xué)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
shuitong110: 金幣+15, ★★★很有幫助, 謝謝 2014-05-08 21:09:57
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-05-09 08:11:01

引用回帖:

4樓: Originally posted by shuitong110 at 2014-04-25 14:52:37
圖示函數(shù)

函數(shù).jpg

首先在你的程序里面，用 trapz 進(jìn)行積分的話，最后得到的 S=trapz(L,M)， S 只是一個(gè)數(shù)值，這是為什么你 plot(L,sin(M)); 只能得到一條線的原因
回頭我再看看怎么得到你說的正弦曲線

贊一下

回復(fù)此樓

5樓2014-04-25 15:03:43

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

jinjian0120

鐵蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 16.9
帖子: 73
在線: 12.4小時(shí)
蟲號(hào): 1746222
注冊(cè): 2012-04-09
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

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

★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
shuitong110: 金幣+3 2014-05-08 21:10:12
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-05-09 08:11:08

幾個(gè)含有參數(shù)矩陣相乘得到的一個(gè)多項(xiàng)式，怎么能讓參數(shù)連續(xù)變化，輸出多項(xiàng)式的值，例如這個(gè)例子x=1:1:10,A=[x,0;0,0],b=A^2 謝謝了

贊一下

回復(fù)此樓

7樓2014-04-25 17:31:06

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

shuitong110

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 209.3
散金: 1042
紅花: 6
帖子: 263
在線: 85.7小時(shí)
蟲號(hào): 2772293
注冊(cè): 2013-11-02
性別: MM
專業(yè): 半導(dǎo)體微納機(jī)電器件與系統(tǒng)

引用回帖:

5樓: Originally posted by alorsmail at 2014-04-25 15:03:43
首先在你的程序里面，用 trapz 進(jìn)行積分的話，最后得到的 S=trapz(L,M)， S 只是一個(gè)數(shù)值，這是為什么你 plot(L,sin(M)); 只能得到一條線的原因
回頭我再看看怎么得到你說的正弦曲線...

萬分感謝啊，最后結(jié)果就是一個(gè)類似正弦的！謝謝啦等您回復(fù)~

贊一下

回復(fù)此樓

8樓2014-04-25 19:31:48

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

shuitong110

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 209.3
散金: 1042
紅花: 6
帖子: 263
在線: 85.7小時(shí)
蟲號(hào): 2772293
注冊(cè): 2013-11-02
性別: MM
專業(yè): 半導(dǎo)體微納機(jī)電器件與系統(tǒng)

引用回帖:

6樓: Originally posted by alorsmail at 2014-04-25 16:12:52
如果要得到一個(gè)曲線的話，你需要加一個(gè)循環(huán)，
示例如下，結(jié)果見圖片
==========================
L=0:0.02*10^-6:140*10^-6 ;

d=0.1*10^-6 ;
R=100*10^-6 ;
P=0.3 ;
k=2*pi./(1.55*10^-6) ;
B=11*10^6 ;
...

謝謝您的回復(fù)，
MM = zeros(size(L)) ;
for i = 10:length(L)
Lp = L(1:i) ;
這三行是什么意思？

贊一下

回復(fù)此樓

9樓2014-04-25 19:41:16

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

shuitong110

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 209.3
散金: 1042
紅花: 6
帖子: 263
在線: 85.7小時(shí)
蟲號(hào): 2772293
注冊(cè): 2013-11-02
性別: MM
專業(yè): 半導(dǎo)體微納機(jī)電器件與系統(tǒng)

引用回帖:

不知道怎么給你發(fā)圖片，我在帖子回復(fù)中加了一個(gè)圖示，是書上的類似結(jié)果，只是參數(shù)不太一樣，但是基本形式不應(yīng)有太大變化。

贊一下

回復(fù)此樓

10樓2014-04-25 19:43:59

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

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

[考研] 304求調(diào)劑 +5	曼殊2266 2026-03-18	5/250	2026-03-20 09:00 by ZHANG0tao
[考研] 085600材料與化工求調(diào)劑 +17	enenenhui 2026-03-13	18/900	2026-03-20 08:55 by xingguangj
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 梁成偉老師課題組歡迎你的加入 +9	一鴨鴨喲 2026-03-14	11/550	2026-03-19 17:22 by ！本暗一次！
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 299求調(diào)劑 +5	△小透明* 2026-03-17	5/250	2026-03-18 11:49 by 盡舜堯1
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識(shí)可樂 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike