 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 3508 | 回復(fù): 2 | ||
[求助]
matlab解一元一次方程的問題 已有1人參與
|
|
matlab解一元一次方程的時候，用到fzero這個命令，但是一元一次方程中我有一個矩陣，matlab總是提示錯誤，但是當(dāng)我把矩陣中的每一個值依次手動帶入進去算就會出結(jié)果。。。這怎么解決呢。。。如果矩陣就幾個值還可以手動挨個帶入算，這要是幾百個我會累死噠。。。。。 比方已知矩陣c=[0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 2 3 4 5 6 7 8 9 10 20 30 40 50 60 70 80 90 100]; 用fzero解下面這個方程，零點大概位置都在1000附近, Q=fzero('x./(sqrt(x.^2+c(1,1))).*besselj(-1,x*1.75*10^-3).*besseli(0,1.75*10^-3.*sqrt(c(1,1)+x.^2))-besselj(0,x*1.75*10^-3).*besseli(-1,sqrt(c(1,1)+x.^2)*1.75*10^-3)',1000)； 這樣解matlab就報錯。。。： ??? Error using ==> fzero at 307 FZERO cannot continue because user supplied expression ==> x./(sqrt(x.^2+c(1,1))).*besselj(-1,x*1.75*10^-3).*besseli(0,1.75*10^-3.*sqrt(c(1,1)+x.^2))-besselj(0,x*1.75*10^-3).*besseli(-1,sqrt(c(1,1)+x.^2)*1.75*10^-3) failed with the error below. Error in inline expression ==> x./(sqrt(x.^2+c(1,1))).*besselj(-1,x*1.75*10^-3).*besseli(0,1.75*10^-3.*sqrt(c(1,1)+x.^2))-besselj(0,x*1.75*10^-3).*besseli(-1,sqrt(c(1,1)+x.^2)*1.75*10^-3) Undefined function or method 'c' for input arguments of type 'double'. 如果我把c（1,1）帶入就沒問題，就出結(jié)果了： Q=fzero('x./(sqrt(x.^2+0.1)).*besselj(-1,x*1.75*10^-3).*besseli(0,1.75*10^-3.*sqrt(0.1+x.^2))-besselj(0,x*1.75*10^-3).*besseli(-1,sqrt(0.1+x.^2)*1.75*10^-3)',1000) 這樣就木有問題。。 我想問下怎樣才能不挨個代數(shù)進去。。。 |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有285人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有68人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
whyjackeyson
金蟲 (小有名氣)
- 應(yīng)助: 20 (小學(xué)生)
- 金幣: 574.1
- 散金: 2
- 紅花: 3
- 帖子: 83
- 在線: 32.3小時
- 蟲號: 3529482
- 注冊: 2014-11-10
- 性別: GG
- 專業(yè): 傳遞過程
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
會飛的泡123: 金幣+20, ★★★★★最佳答案 2015-11-09 21:20:09
感謝參與，應(yīng)助指數(shù) +1
會飛的泡123: 金幣+20, ★★★★★最佳答案 2015-11-09 21:20:09
|
fzero函數(shù)使用錯誤，對于存在符號變量的求解使用fzero可以采用句柄方式代入，或者外部定義函數(shù)。 求解可使用循環(huán)語句。 c = [0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 2 3 4 5 6 7 8 9 10 20 30 40 50 60 70 80 90 100]; for i = 1:28 Q(i) = fzero(@(x)x./(sqrt(x.^2+c(i))).*besselj(-1,x*1.75*10^-3).*besseli(0,1.75*10^-3.*sqrt(c(i)+x.^2))-besselj(0,x*1.75*10^-3).*besseli(-1,sqrt(c(i)+x.^2)*1.75*10^-3),1000); end |
2樓2015-11-09 18:31:07
3樓2015-11-09 21:19:50
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 會飛的泡123 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷 +14 | 曦熙兮 2026-03-15 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 生物學(xué)調(diào)劑 +3 | Surekei 2026-03-21 | 3/150 |
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4 | lemonzzn 2026-03-17 | 7/350 |
|
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 | ￣^￣゜汗 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|