matlab求助,將一些帶變數的方程代入到另方程裡,大神都在睡覺嗎?跪求了

2021-04-17 22:04:29 字數 2638 閱讀 6509

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