1樓:匿名使用者
% 知道1206
clear
clcclose all
lamda_p=0.532;
th=[0:pi/50:pi/2];
lamda_s=[0.54:0.1:1.064]; %lamda_s是一個變值,從0.54一直迴圈到1.064
no_p=4.7359+0.01878./(lamda_p^2-0.1822);%no_p是定值,因為lamda_p已知
ne_p=4.3753+0.01224./(lamda_p^2-0.1667);%ne_p也是定值,lamda_p已知
for m=1:length(th)
for n=1:length(lamda_s)
lamda_i=lamda_s(n)*lamda_p/(lamda_s(n)-lamda_p);%lamda_p已知,所以假設給lamda_s取一個值就能得到一個lamda_i
no_i=4.7359+0.01817./((lamda_i)^2-0.1822);% no_i只與lamda_i有關,lamda_i變它也變
ne_i=4.3753+0.01224./((lamda_i)^2-0.1667);%ne_i只與lamda_i有關lamda_i變它也變
no_s=4.7359+0.01878./(lamda_s(n)^2-0.1822);%no_s只與lamda_s有關,lamda_s變它也變
%問題來了,我知不道如何將上面得到的no_i ,ne_i,no_p,ne_p,no_s代人下面的一個式子中,從而得到一系列的neth_p 和neth_i.
neth_p=(cos(th(m))^2/(no_p^2)+sin(th(m))^2./(ne_p^2))^(-0.5);
neth_i=(cos(th(m))^2/(no_i^2)+sin(th(m))^2./(ne_i^2))^(-0.5);
% 下面兩行是你原來的式子,我感覺不對,就修改為上面的式子了。
% neth_p=(cos(m)^2/(no_p^2)+sin(m)^2./(ne_p^2))^(-0.5);
% neth_i=(cos(m)^2/(no_i^2)+sin(m)^2./(ne_i^2))^(-0.5);
%最後將一系列的neth_p 和neth_i代入if這個條件裡,滿足就輸出。來得到 th與lamda_i的影象和 th與lamda_s的影象。
error=neth_p/lamda_p-no_s/lamda_s(n)-neth_i/lamda_i;
if error<=1e-4
lamda_i(m)=lamda_i;
lamda_s(m)=lamda_s(n);
break
endend
endplot(th,lamda_i,'r-*',th,lamda_s,'-.','linewidth',2)
legend('th-lamda-i','th-lamda-s')
grid
ylim([-2,40])
set(gca,'fontsize',14)
2樓:匿名使用者
不能代入,除非你用符號
matlab中如何將一個含有多個變數的方程的某一變數放在方程的一邊,而其他變數放在方程的另一邊? 15
3樓:
有積分變數,它肯定是能求導的,所以把那個積分變數看作是x4,它的導數是x4',從而整個微分方程是一個四階的,相當於括了一階,再進行化簡和整理,用ode45函式解微分方程,一定行!!
matlab程式設計求教,含參變數的方程,給出一系列的參變數的值,求方程相應的解 20
4樓:匿名使用者
樓上 @兔子和小強 從理論上進行了分析,下面我從另外的角度進行說明。
1、可以用控制理論裡面的根軌跡方法來求出a=2:50時方程根的分佈情況:
s=tf('s');
rlocus(s^7/(s^2+1)^7,2:50)
由圖可見,對於給定的a值,方程根本就沒有實數根,更不要說正數解了。如果把a的範圍放大到0~inf的範圍,得到完整的根軌跡如下:
rlocus(s^7/(s^2+1)^7)
ylim([-3 3])
axis equal
可見,對於a>0,方程不存在正數根。對於a<0的情況,可以用下面的命令繪圖:
rlocus(-s^7/(s^2+1)^7)
ylim([-3 3])
axis equal
可見,只有a<0的時候,方程才可能出現正數解(進一步可以確定,方程存在正數解的要求大約是a<-128)。
2、上面用到了控制系統工具箱的rlocus函式,如果未安裝該工具箱,或者對根軌跡的概念不熟悉,也可以用roots函式直接求根:
syms a x
f=a*x^7+(x^2+1)^7;
r=arrayfun(@(a),2:50);
plot(cell2mat(r).')
繪圖結果與上面第一張圖一致。
5樓:兔子和小強
該方程的解是 x^2 - a^(1/7)x + 1 = 0 的解。利用二次方程根的公式可直接寫出解。
但是,當 a = 2:50 時,δ<0,方程不存在實數解。
求一些帶京腔京韻的歌曲,求推薦一些帶京腔京韻的歌曲
張傑的 心bai中一把火du 最後有一句,李玉zhi剛的,周杰倫的dao 雨下一整晚 後半專部分 王力巨集 蓋屬世英雄 陶喆 susan說 王佩 夢北京 李宇春 n 1 魏晨 少年遊 信樂團 one night in beijing霍元甲 周杰倫 梨園英雄 華少翌 陳紅 深夜走過長安街 tank 三...
求助一些六上的數學問題
1.某一年,三年定期存款的年利率是,到期還要繳納20 的利息稅。李明存了1000元,定期三年,到期後可的稅後利息 元。稅前利息是。1000 元。稅後利息是。81 1 20 元。2.為名小學今年招收新生160人,今年新生人數的10分之1相當於去年新生人數的9分之1,去年招收新生多少人?160 10分之...
初二上的一些題 求助求助!!幫幫忙!!
1 3x 2 x 1 3x 2 x 2,比較係數得。m 1,n 2 2由 x 1 2 2 x 2 x 1 4,k 1 43 原式 8x 2y 2 x 3y m 34 選d5 x y 2 x 2 xy y 2 4,m y 2 16 6 1 原式 a 2b a 2b 3a 2 2b 2a 2 4 b a...