 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 仿真模擬 » MATLAB » 分段函數(shù)的積分問題
| 1 | 1/1 | 返回列表 |
| 查看: 725 | 回復: 0 | |||
narcissus璇銀蟲 (小有名氣)
|
[交流]
分段函數(shù)的積分問題
|
|
求助一下大神，我現(xiàn)在這個代碼能夠畫出一個Ey的分段函數(shù)的圖像，求問一下，我想求這個分段函數(shù)Ey在兩個定值x0到x4的定積分，接下來的代碼應該怎么寫，我之前用了int函數(shù)，總是運行不出來結(jié)果。 dclear%TE close all n1 = 3.476; ns = 1.45; n0 = 1.45; d = 400e-9; a=d/2; lambda = 1.55e-6; k = 2*pi/lambda; e=2.718281828459; r=(ns^2-n0^2)/(n1^2-ns^2); v=sqrt(k^2*a^2*(n1^2-ns^2)); m=0; F1=@(b)(2*v*sqrt(1-b)); ezplot(F1,[0,1,-5,10]); hold on; F2=@(b)(m*pi)+atan(sqrt(b/(1-b)))+atan (sqrt((b+r)/(1-b))); ezplot(F2,[0,1,-5,10]); xlabel('b') Feigin=@(b)(F1(b)-F2(b)); b(1)=fzero(Feigin,0.5); ne=sqrt(b*n1^2+(1-b)*ns^2); beta=ne*k; kappa=sqrt((k^2*n1^2)-beta^2); xi=sqrt(beta^2-k^2*ns^2); sigma=sqrt(beta^2-k^2*n0^2); phi=m*pi/2+atan(sqrt(b/(1-b)))/2-atan(sqrt((b+r)/(1-b)))/2; format short g [ne kappa xi sigma phi]; Np=1001; x2=linspace(-12*a,-a,Np); x1=linspace(-a,a,Np); x3=linspace(a,12*a,Np); Ey1=zeros(1,Np); Ey2=Ey1; Ey3=Ey1; Ey1=cos(kappa*x1-phi); Ey2=cos(kappa*a+phi)*exp(xi*(x2+a)); Ey3=cos(kappa*a-phi)*exp(-sigma*(x3-a)); Ey=[Ey2,Ey1,Ey3]; Ey=Ey/diag(max(abs(Ey))); x=[x2,x1,x3]; hold on %plot(x,Ey,'g+'); plot(x1,Ey1,'g',x2,Ey2,'r',x3,Ey3,'b') %legend('TE_0','TE_1','TE_2') xlabel('x') axis([-12e-7 12e-7 -1.1 1.1]) %plot([-d,-d],[-1.1,1.1],'black--') plot([ 0, 0],[-1.1,1.1],'black--') plot([x(1), x(end)],[0,0],'black') plot([x(1), x(end)],[1/e,1/e],'black') x0=-1/xi*log(e*cos(kappa*a))-a; x4=1/sigma*log(e*cos(kappa*a))+a; |
回復此樓» 猜你喜歡
070303一志愿西北大學學碩310找調(diào)劑
已經(jīng)有9人回復
-
0703化學調(diào)劑，求各位老師收留
已經(jīng)有9人回復
-
0703化學調(diào)劑
已經(jīng)有11人回復
-
0703化學 305求調(diào)劑
已經(jīng)有3人回復
-
一志愿吉林大學材料學碩321求調(diào)劑
已經(jīng)有5人回復
-
274求調(diào)劑
已經(jīng)有3人回復
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有3人回復
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有9人回復
-
【同濟軟件】軟件（085405）考研求調(diào)劑
已經(jīng)有3人回復
-
267一志愿南京工業(yè)大學0817化工求調(diào)劑
已經(jīng)有7人回復
相關版塊跳轉(zhuǎn)
我要訂閱樓主 narcissus璇 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +3 | Ymlll 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 267一志愿南京工業(yè)大學0817化工求調(diào)劑 +7 | SUICHILD 2026-03-12 | 7/350 |
|
|
|
[考研] 化學工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 314求調(diào)劑 +8 | 無懈可擊的巨人 2026-03-12 | 8/400 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 0703化學調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 本人考085602 化學工程 專碩 +16 | 不知道叫什么！ 2026-03-15 | 18/900 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 復試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|