 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1753 | 回復(fù): 5 | ||
cologate007新蟲 (初入文壇)
|
[求助]
matlab對含貝塞爾函數(shù)的超越方程求根 已有1人參與
|
|
有一個含有貝塞爾函數(shù)的超越方程，用solve函數(shù)求解，輸入matlab直接給出了一個復(fù)數(shù)結(jié)果，但是實際上超越方程是有很多個根的，不知道給出的這個解是第幾階的？怎么樣才能得到它的第n（任意的）階根呢？采用牛頓迭代法給定不同的初值可以計算得到不同的根，但是也不能確定是第幾階。不知道論壇里有沒有朋友有更好的計算方法？函數(shù)就不貼了，下面是牛頓迭代法的matlab程序，函數(shù)見里面的f。請了解的朋友給點意見，感謝！ clc;clear;close all; r0=0.0125; r=r0; fre=1000; z=0; l=0.122; r0=0.0125; r=r0; rou0=1000; a0=1500; pi=3.1415926; k0=2*pi*f/a0; Zw=1800000*(1+i); %牛頓迭代法 syms f kr; f=subs('besselj(0,kr*r0)-kr*r0*i*Zw/(rou0*a0*k0*r0)*besselj(1,kr*r0)'); df=diff(f); kr0=0.1-0.1*i; e1=0.1e-10; e2=0.000005; N=100; k=1; while (k<N) kr=kr0; if abs(eval(f))<e1 fprintf('奇異!\nkr=%.6f\n迭代次數(shù)為:%d\n',kr0,k) break else kr1=kr0-eval(f)/eval(df); if abs(kr1-kr0)<e2 fprintf('kr=%.6f\n迭代次數(shù)為:%d\n',kr1,k) break else kr0=kr1; k=k+1; end end end if k>=N fprintf('失敗\n') end |
回復(fù)此樓» 猜你喜歡
285求調(diào)劑
已經(jīng)有5人回復(fù)
-
289求調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
【求調(diào)劑】085601材料工程專碩 | 總分272 |
已經(jīng)有4人回復(fù)
-
286求調(diào)劑
已經(jīng)有8人回復(fù)
-
083000學(xué)碩274求調(diào)劑
已經(jīng)有3人回復(fù)
-
085701環(huán)境工程求調(diào)劑
已經(jīng)有4人回復(fù)
-
070300求調(diào)劑306分
已經(jīng)有3人回復(fù)
-
266求調(diào)劑
已經(jīng)有12人回復(fù)
-
化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助matlab寫變限積分
已經(jīng)有4人回復(fù)
- 求助：如何求含貝塞爾函數(shù)的超越方程的復(fù)數(shù)解
已經(jīng)有3人回復(fù)
- 有關(guān)于mathematics的問題，用mathematics對一些點擬合出一個橢圓方程
已經(jīng)有24人回復(fù)
- [求助]求解極座標(biāo)上的偏微分方程(Laplace’s Equation)!!!!
已經(jīng)有6人回復(fù)
- 求教關(guān)于origin8.0 自定義函數(shù)擬合方面的問題
已經(jīng)有5人回復(fù)
- 急！在用origin進(jìn)行非線性擬合怎樣創(chuàng)建新的擬合函數(shù)？
已經(jīng)有14人回復(fù)
- 師兄師姐們，幫小弟看看這個 Laplace變換后的二階常微分方程怎么解啊
已經(jīng)有23人回復(fù)
- VS2012如何快速輸入函數(shù)？
已經(jīng)有6人回復(fù)
- origin 里 “三次樣條函數(shù)擬合” 與“三次多項式擬合”是一回事嗎？如圖
已經(jīng)有3人回復(fù)
- 用matlab求解含有bessel函數(shù)的方程
已經(jīng)有5人回復(fù)
- 求解高次復(fù)雜函數(shù)
已經(jīng)有11人回復(fù)
- 微分方程求解（通解中含hypergeom函數(shù)）
已經(jīng)有17人回復(fù)
- 用matlab求bessel函數(shù)的根
已經(jīng)有6人回復(fù)
- 一個高一函數(shù)題求解
已經(jīng)有7人回復(fù)
- 【求助】復(fù)超越方程的解法
已經(jīng)有19人回復(fù)
mathstudy
金蟲 (正式寫手)
- 數(shù)學(xué)EPI: 2
- 應(yīng)助: 161 (高中生)
- 金幣: 4320.1
- 散金: 38
- 紅花: 16
- 帖子: 446
- 在線: 141.1小時
- 蟲號: 2515489
- 注冊: 2013-06-20
- 性別: GG
- 專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
2樓2015-03-23 09:33:06
cologate007
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 407.5
- 帖子: 10
- 在線: 43.1小時
- 蟲號: 1522241
- 注冊: 2011-12-03
- 專業(yè): 聲學(xué)
3樓2015-03-23 14:15:33
mathstudy
金蟲 (正式寫手)
- 數(shù)學(xué)EPI: 2
- 應(yīng)助: 161 (高中生)
- 金幣: 4320.1
- 散金: 38
- 紅花: 16
- 帖子: 446
- 在線: 141.1小時
- 蟲號: 2515489
- 注冊: 2013-06-20
- 性別: GG
- 專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
4樓2015-03-24 09:44:38
cologate007
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 407.5
- 帖子: 10
- 在線: 43.1小時
- 蟲號: 1522241
- 注冊: 2011-12-03
- 專業(yè): 聲學(xué)
5樓2015-03-24 14:20:53
yinmengyuhen
鐵桿木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 7314
- 帖子: 183
- 在線: 97.6小時
- 蟲號: 1722129
- 注冊: 2012-03-28
- 性別: MM
- 專業(yè): 力學(xué)
6樓2019-04-10 10:11:46
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cologate007 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +7 | 新時代材料 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 286求調(diào)劑 +8 | PolarBear11 2026-03-26 | 8/400 |
|
|
|
[考研] 291求調(diào)劑 +12 | hhhhxn.. 2026-03-23 | 18/900 |
|
|
|
[考研] 265求調(diào)劑 +5 | 小木蟲085600 2026-03-27 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +6 | 嗯嗯嗯嗯2 2026-03-27 | 6/300 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +5 | SYA！ 2026-03-23 | 5/250 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 329求調(diào)劑 +7 | 鈕恩雪 2026-03-25 | 7/350 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +6 | 紅衣隱官 2026-03-21 | 6/300 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 10/500 |
|
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +5 | WW.' 2026-03-23 | 7/350 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|