notdingy

鐵蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 3963.2
散金: 27
帖子: 680
在線: 31.8小時
蟲號: 3429502
注冊: 2014-09-20
性別: MM
專業(yè): 光學

[求助] 利用MATLAB計算如下方程組，求大神算出結(jié)果。解析法顯示為空，可能要利用數(shù)值法已有3人參與

L1='0.8334=(x^2*(sinh(sqrt(x^2-(2*pi*y*(1/1544.14-1/z)^2)))*u)^2)/(x^2*(cosh(sqrt(x^2-(2*pi*y*(1/1544.14-1/z)^2)))*u)^2-(2*pi*y*(1/1544.14-1/z)^2))';
L2='0.8799=(x^2*(sinh(sqrt(x^2-(2*pi*y*(1/1544.145-1/z)^2)))*u)^2)/(x^2*(cosh(sqrt(x^2-(2*pi*y*(1/1544.145-1/z)^2)))*u)^2-(2*pi*y*(1/1544.145-1/z)^2))';
L3='0.91092=(x^2*(sinh(sqrt(x^2-(2*pi*y*(1/1544.15-1/z)^2)))*u)^2)/(x^2*(cosh(sqrt(x^2-(2*pi*y*(1/1544.15-1/z)^2)))*u)^2-(2*pi*y*(1/1544.15-1/z)^2))';
L4='0.92915=(x^2*(sinh(sqrt(x^2-(2*pi*y*(1/1544.155-1/z)^2)))*u)^2)/(x^2*(cosh(sqrt(x^2-(2*pi*y*(1/1544.155-1/z)^2)))*u)^2-(2*pi*y*(1/1544.155-1/z)^2))';
solve(L1,L2,L3,L4,x,y,z,u)

x，y，z，u為未知數(shù)，四個方程形式都一樣。利用fsolve？

回復此樓

» 猜你喜歡

哈爾濱理工大學物理系招收物理學考研調(diào)劑已經(jīng)有7人回復
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有111人回復
求調(diào)劑已經(jīng)有0人回復
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復
光學工程學碩調(diào)劑信息已經(jīng)有31人回復
歡迎加入課題組已經(jīng)有0人回復
散金幣，求好運，祝面上順利！已經(jīng)有68人回復
華東師范大學芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復
調(diào)劑已經(jīng)有0人回復

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

如何利用MATLAB模擬建筑物，用射線跟蹤法計算接收機狀態(tài) 已經(jīng)有7人回復
怎么利用matlab計算如下函數(shù)的積分已經(jīng)有8人回復
matlab計算數(shù)據(jù)的鄰接矩陣已經(jīng)有0人回復
MATLAB計算問題已經(jīng)有0人回復
matlab計算晶相含量已經(jīng)有3人回復
求一份利用蒙特卡羅方法，并結(jié)合MATLAB計算光在膠體中散射系數(shù)的代碼已經(jīng)有1人回復
求教如何利用matlab 模擬計算兩塊異性磁體之間磁感應強度分布？已經(jīng)有0人回復
如何利用matlab 模擬計算兩塊異性磁體之間磁感應強度分布？已經(jīng)有0人回復
利用MATLAB計算熱管軸向參數(shù)的問題已經(jīng)有1人回復
怎么利用matlab計算出非線性函數(shù)模型的系數(shù)？急求�。。。�！已經(jīng)有4人回復
利用Matlab模擬空間濾波實驗已經(jīng)有43人回復
計算數(shù)據(jù)，利用matlab畫三維曲面圖已經(jīng)有2人回復
如何利用Matlab編程計算熱加工圖中功率耗散和失穩(wěn)值（利用Origin畫圖）已經(jīng)有0人回復
【求助】如何利用MATLAB計算N元N次方程已經(jīng)有9人回復
【求助】如何利用matlab計算這個無窮級數(shù)的和已經(jīng)有12人回復

1樓 2015-11-17 19:55:08

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

dingd

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

計算強帖: 4
應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應助回帖

★ ★
感謝參與，應助指數(shù) +1
notdingy: 金幣+2, ★★★很有幫助 2015-11-18 17:39:05

數(shù)值解的話用1stOpt試試，很簡單的。

[ 發(fā)自手機版 http://www.gaoyang168.com/3g ]

贊一下(1人)

回復此樓

2樓2015-11-17 22:51:58

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

凌風7388

鐵蟲 (正式寫手)

應助: 4 (幼兒園)
金幣: 1153.2
散金: 72
紅花: 5
帖子: 573
在線: 160.6小時
蟲號: 759985
注冊: 2009-04-29
專業(yè): 化學反應工程

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
notdingy: 金幣+8, ★★★很有幫助 2015-11-18 17:38:48

用fsolve可以解，但是解之間需要給定初值，初值不同得到的結(jié)果就不同。
比如假設(shè)x=y=z=u=1 1 1 1，得到的結(jié)果依次為：1.8634 0.3599 1.9406 1.1542
比如假設(shè)x=z=1，y=u=2，得到的結(jié)果依次為：1.6796 2.0456 0.8471 2.4493

贊一下(1人)

回復此樓

3樓2015-11-18 09:36:02

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

notdingy

鐵蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 3963.2
散金: 27
帖子: 680
在線: 31.8小時
蟲號: 3429502
注冊: 2014-09-20
性別: MM
專業(yè): 光學

引用回帖:

3樓: Originally posted by 凌風7388 at 2015-11-18 09:36:02
用fsolve可以解，但是解之間需要給定初值，初值不同得到的結(jié)果就不同。
比如假設(shè)x=y=z=u=1 1 1 1，得到的結(jié)果依次為：1.8634 0.3599 1.9406 1.1542
比如假設(shè)x=z=1，y=u=2，得到的結(jié)果依次為：1.6796 ...

具體的程序要怎么寫呢，沒用過fsolve

贊一下

回復此樓

4樓2015-11-18 10:39:19

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

notdingy

鐵蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 3963.2
散金: 27
帖子: 680
在線: 31.8小時
蟲號: 3429502
注冊: 2014-09-20
性別: MM
專業(yè): 光學

引用回帖:

2樓: Originally posted by dingd at 2015-11-17 22:51:58
數(shù)值解的話用1stOpt試試，很簡單的。

網(wǎng)上找了很多1stopt，都安裝不成功

贊一下

回復此樓

5樓2015-11-18 10:40:19

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

dingd

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

計算強帖: 4
應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

引用回帖:

5樓: Originally posted by notdingy at 2015-11-18 10:40:19
網(wǎng)上找了很多1stopt，都安裝不成功...

壇子里就有啊

回復此樓

6樓2015-11-18 11:44:10

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

1314168apple

金蟲 (知名作家)

應助: 68 (初中生)
金幣: 677
紅花: 12
帖子: 6872
在線: 1462.3小時
蟲號: 287760
注冊: 2006-10-21
專業(yè): 色譜分析

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
notdingy: 金幣+10, ★★★很有幫助 2015-11-18 17:38:30

fun=@(x)[0.8334-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2));0.8799-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.145-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.145-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.145-1/x(3))^2));0.91092-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.15-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.15-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.15-1/x(3))^2));0.92915-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.155-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.155-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.155-1/x(3))^2))];

TT=fsolve(fun,randn(4,1));

都沒有精確的解，可能不存在？檢查檢查。
誤差(好些解得到的差不多)：

fun(TT)

-0.0549
-0.0084
0.0226
0.0408

贊一下

回復此樓

為了最終理解你所不理解的，你必須經(jīng)歷一條愚昧無知的道路。為了占有你從未占有的東西，你必須經(jīng)歷被剝奪的道路。為了達到你現(xiàn)在所不在的名位，你必須經(jīng)歷那...

7樓2015-11-18 12:19:55

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

凌風7388

鐵蟲 (正式寫手)

應助: 4 (幼兒園)
金幣: 1153.2
散金: 72
紅花: 5
帖子: 573
在線: 160.6小時
蟲號: 759985
注冊: 2009-04-29
專業(yè): 化學反應工程

引用回帖:

4樓: Originally posted by notdingy at 2015-11-18 10:39:19
具體的程序要怎么寫呢，沒用過fsolve...

function xFsolve
clear all;clc
x0=[1 1 1 1]';x1=fsolve(@nonlineqs,x0)
x0=[1 2 1 2]';x2=fsolve(@nonlineqs,x0)

function f=nonlineqs(x)
f(1)=0.8334-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2));
f(2)=0.8799-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2));
f(3)=0.91092-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2));
f(4)=0.92915-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2));

贊一下

回復此樓

8樓2015-11-18 13:28:44

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

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

[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 0805 316求調(diào)劑 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 學員8dgXkO
[考研] 材料學碩333求調(diào)劑 +3	北道巷 2026-03-18	3/150	2026-03-21 18:17 by 學員8dgXkO
[考研] 278求調(diào)劑 +9	煙火先于春 2026-03-17	9/450	2026-03-21 17:47 by 學員8dgXkO
[考研] 0703化學297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學員8dgXkO
[考研] 346求調(diào)劑[0856] +4	WayneLim327 2026-03-16	7/350	2026-03-21 04:02 by JourneyLucky
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 材料專碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 295復試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染