wh1125

木蟲 (正式寫手)

應(yīng)助: 10 (幼兒園)
金幣: 3078.9
紅花: 5
帖子: 498
在線: 431.2小時(shí)
蟲號(hào): 417863
注冊(cè): 2007-07-02
性別: GG
專業(yè): 集成電路設(shè)計(jì)與測(cè)試

[求助] 求助：matlab求解一個(gè)方程組，其中一個(gè)方程有三角函數(shù)

方程組有5個(gè)方程，前4個(gè)方程比較簡(jiǎn)單，就是最后一個(gè)方程有三角函數(shù)，不知道該如何求解，請(qǐng)大神們幫幫忙�。「屑げ槐M�。�！

我寫的代碼：
syms K0 K1 K2 K3 beta0 k0 n0 n1 n2 n3

k0=2*pi/0.88;
n0=3.4;
n1=3.6;
n2=3.4;
n3=3.55;
[K0 K1 K2 K3 beta0]=solve('K0^2=k0^2*n0^2-beta0^2','K1^2=k0^2*n1^2-beta0^2','K2^2=k0^2*n2^2-beta0^2','K3^2=k0^2*n3^2-beta0^2',...
'(K1^2*K2*K3^2+K0*K3^4+K0*K1^2*K2^2+K0^2*K2*K3^2)*(tan(K3*0.2))^2=K3^2*(K0^2*K2+K0*K2^2+K0*K1^2+K1^2*K2)','K0','K1','K2','K3','beta0')

方程組.png

[ Last edited by wh1125 on 2013-9-26 at 09:22 ]

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有151人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(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)問(wèn)還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有26人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有34人回復(fù)

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

solve解超越方程少了一個(gè)解已經(jīng)有3人回復(fù)
請(qǐng)教一個(gè)matlab求解非線性方程組的問(wèn)題已經(jīng)有9人回復(fù)
matlab如何求解一個(gè)非線性微分方程組已經(jīng)有8人回復(fù)
求助，用matlab解微分方程組，希望高手能給予指點(diǎn) 已經(jīng)有17人回復(fù)
這個(gè)方程組怎么解啊已經(jīng)有14人回復(fù)
MATLAB求解非線性方程組已經(jīng)有5人回復(fù)
求Matlab高手解決線性方程組的迭代求解問(wèn)題已經(jīng)有17人回復(fù)
如何求解下列方程：xlnx=a(x-b) 已經(jīng)有14人回復(fù)
各位朋友，如何用matlab編寫這個(gè)方程，并求解？已經(jīng)有15人回復(fù)
求助 MATLAB解方程組-fslove 已經(jīng)有7人回復(fù)
求助matlab---fsolve解非線性方程組已經(jīng)有6人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復(fù)
高金求助matlab解微分方程組已經(jīng)有12人回復(fù)
求matlab求解一個(gè)方程組的問(wèn)題已經(jīng)有6人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
MATLAB解方程出現(xiàn)??? Subscripted assignment dimension mismatch. 已經(jīng)有3人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題已經(jīng)有21人回復(fù)
【求助】用mathematica 5.0求解一個(gè)非線性方程組失敗，特發(fā)帖求助！已經(jīng)有5人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
求助一個(gè)差分方程的MATLAB實(shí)現(xiàn) 已經(jīng)有3人回復(fù)
matlab求解超越方程已經(jīng)有3人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)

1樓 2013-09-25 22:51:38

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

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

月只藍(lán)

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

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
wh1125: 金幣+28, ★★★★★最佳答案, 太感謝了！我對(duì)你的景仰猶如滔滔江水連綿不絕�。� 2013-09-26 19:12:20

function solve_nonl
clear all;clc
x0=[-7 7 7 6.5 24];
options=optimset('MaxIter',1000,'MaxFunEvals',5000);
x=fsolve(@nonl,x0,options)

function f=nonl(x)
K0=x(1);
K1=x(2);
K2=x(3);
K3=x(4);
beta0=x(5);
k0=2*pi/0.88;
n0=3.4;
n1=3.6;
n2=3.4;
n3=3.55;
f(1)=K0^2-K0^2*n0^2+beta0^2;
f(2)=K1^2-k0^2*n1^2+beta0^2;
f(3)=K2^2-k0^2*n2^2+beta0^2;
f(4)=K3^2-k0^2*n3^2+beta0^2;
f(5)=(K1^2*K2*K3^2+K0*K3^4+K0*K1^2*K2^2+K0^2*K2*K3^2)*(tan(K3*0.2))^2-(K3^2*(K0^2*K2+K0*K2^2+K0*K1^2+K1^2*K2));

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

6樓2013-09-26 18:17:37

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

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

試試1stOpt，不用初值（1stOpt不區(qū)分大小寫，小寫k0用k01代替）：

CODE:

Constant k01=2*pi/0.88,n0=3.4,n1=3.6,n2=3.4,n3=3.55;

Function

K0^2=k01^2*n0^2-beta0^2;

K1^2=k01^2*n1^2-beta0^2;

K2^2=k01^2*n2^2-beta0^2;

K3^2=k01^2*n3^2-beta0^2;

(K1^2*K2*K3^2+K0*K3^4+K0*K1^2*K2^2+K0^2*K2*K3^2)*(tan(K3*0.2))^2=K3^2*(K0^2*K2+K0*K2^2+K0*K1^2+K1^2*K2);

k0: -9.57452012855549E-12
beta0: -24.2759432314211
k1: -8.44814217715797
k2: 7.09783257508414E-12
k3: 7.29012925007402

贊一下(1人)

回復(fù)此樓

8樓2013-09-26 20:39:37

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

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

1stOpt是個(gè)優(yōu)化計(jì)算軟件，與Matlab無(wú)關(guān)，谷歌一下。

贊一下(1人)

回復(fù)此樓

10樓2013-09-26 22:07:42

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

普通回帖

lihongyu0526

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1252.3
帖子: 176
在線: 52小時(shí)
蟲號(hào): 191329
注冊(cè): 2006-02-20
性別: GG
專業(yè): 無(wú)機(jī)材料化學(xué)

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

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

你有6個(gè)未知數(shù)，卻只有5個(gè)方程。

贊一下

回復(fù)此樓

2樓2013-09-26 14:46:02

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

wh1125

木蟲 (正式寫手)

應(yīng)助: 10 (幼兒園)
金幣: 3078.9
紅花: 5
帖子: 498
在線: 431.2小時(shí)
蟲號(hào): 417863
注冊(cè): 2007-07-02
性別: GG
專業(yè): 集成電路設(shè)計(jì)與測(cè)試

引用回帖:

2樓: Originally posted by lihongyu0526 at 2013-09-26 14:46:02
你有6個(gè)未知數(shù)，卻只有5個(gè)方程。

沒有啊，就5個(gè)未知數(shù)，那個(gè)小寫k0上面已經(jīng)賦值了。

贊一下

回復(fù)此樓

3樓2013-09-26 16:03:49

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

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

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

★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
wh1125: 金幣+2, ★有幫助, 沒有解決我的問(wèn)題 2013-09-26 17:41:38
ben_ladeng: 金幣+1, 感謝參與 2013-09-26 19:13:25

1：Solve函數(shù)是求符號(hào)解，很多情況下符號(hào)解是無(wú)法求出的，此時(shí)可改用數(shù)值解命令fsolve；
2：5個(gè)方程卻只有4個(gè)未知數(shù)（k1,k2,k3,beta0），屬超越方程，無(wú)精確解；
3：如果去掉最后一個(gè)方程，可得唯一解：

beta0: 23.2021994252868
k1: -11.0612145738144
k2: -7.13998330841327
k3: -10.2041827740037

贊一下

回復(fù)此樓

4樓2013-09-26 16:12:27

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

wh1125

木蟲 (正式寫手)

應(yīng)助: 10 (幼兒園)
金幣: 3078.9
紅花: 5
帖子: 498
在線: 431.2小時(shí)
蟲號(hào): 417863
注冊(cè): 2007-07-02
性別: GG
專業(yè): 集成電路設(shè)計(jì)與測(cè)試

引用回帖:

4樓: Originally posted by dingd at 2013-09-26 16:12:27
1：Solve函數(shù)是求符號(hào)解，很多情況下符號(hào)解是無(wú)法求出的，此時(shí)可改用數(shù)值解命令fsolve；
2：5個(gè)方程卻只有4個(gè)未知數(shù)（k1,k2,k3,beta0），屬超越方程，無(wú)精確解；
3：如果去掉最后一個(gè)方程，可得唯一解：

beta0 ...

有5個(gè)未知數(shù)，那個(gè)大寫的K0也是未知數(shù)。

贊一下

回復(fù)此樓

5樓2013-09-26 17:40:24

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

wh1125

木蟲 (正式寫手)

應(yīng)助: 10 (幼兒園)
金幣: 3078.9
紅花: 5
帖子: 498
在線: 431.2小時(shí)
蟲號(hào): 417863
注冊(cè): 2007-07-02
性別: GG
專業(yè): 集成電路設(shè)計(jì)與測(cè)試

引用回帖:

6樓: Originally posted by 月只藍(lán) at 2013-09-26 18:17:37
function solve_nonl
clear all;clc
x0=;
options=optimset('MaxIter',1000,'MaxFunEvals',5000);
x=fsolve(@nonl,x0,options)

function f=nonl(x)
K0=x(1);
K1=x(2);
K2=x(3);
K3=x(4);
beta0=x(5); ...

大神，弱弱的問(wèn)一句：x0=[-7 7 7 6.5 24]，這個(gè)初值你大致是怎么猜的？

贊一下

回復(fù)此樓

7樓2013-09-26 19:20:41

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

wh1125

木蟲 (正式寫手)

應(yīng)助: 10 (幼兒園)
金幣: 3078.9
紅花: 5
帖子: 498
在線: 431.2小時(shí)
蟲號(hào): 417863
注冊(cè): 2007-07-02
性別: GG
專業(yè): 集成電路設(shè)計(jì)與測(cè)試

引用回帖:

8樓: Originally posted by dingd at 2013-09-26 20:39:37
試試1stOpt，不用初值（1stOpt不區(qū)分大小寫，小寫k0用k01代替）：

Constant k01=2*pi/0.88,n0=3.4,n1=3.6,n2=3.4,n3=3.55;
Function
K0^2=k01^2*n0^2-beta0^2;
K1^2=k01^2*n1^2-beta0^2;
K2^2=k01^2*n2^2-be ...

請(qǐng)問(wèn)在我這里1stOpt怎么用��？help 1stOpt都沒有這個(gè)關(guān)鍵字

贊一下

回復(fù)此樓

9樓2013-09-26 21:42:31

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

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

[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +7	yangfz 2026-03-17	7/350	2026-03-17 23:57 by 星空星月
[考研] 299求調(diào)劑 +4	△小透明* 2026-03-17	4/200	2026-03-17 20:09 by peike
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 085600材料與化工求調(diào)劑 +5	緒幸與子 2026-03-17	5/250	2026-03-17 16:40 by laoshidan
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 271求調(diào)劑 +12	生如夏花… 2026-03-11	14/700	2026-03-17 10:56 by lovewei0727
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 286求調(diào)劑 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 311求調(diào)劑 +5	26研0 2026-03-15	5/250	2026-03-16 16:21 by a不易
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[基金申請(qǐng)] 有必要更換申報(bào)口嗎 20+3	fannyamoy 2026-03-11	3/150	2026-03-14 00:52 by zhanghaozhu
[考研] 332求調(diào)劑 +3	zjy101327 2026-03-11	6/300	2026-03-13 22:48 by JourneyLucky
[考研] 0703化學(xué)調(diào)劑 +4	快樂(lè)的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 求調(diào)劑 +7	18880831720 2026-03-11	7/350	2026-03-13 16:10 by JourneyLucky
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[考研] 0856化學(xué)工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant