| 6 | 1/1 | 返回列表 |
| 查看: 1749 | 回復(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 |
金蟲 (正式寫手)
新蟲 (初入文壇)
金蟲 (正式寫手)
新蟲 (初入文壇)
鐵桿木蟲 (小有名氣)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷 +9 | 材料學(xué)257求調(diào)劑 2026-03-23 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +10 | Zhao anqi 2026-03-22 | 10/500 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
[考研] 考研調(diào)劑 +4 | 來好運來來來 2026-03-21 | 4/200 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 一志愿南大,0703化學(xué),分?jǐn)?shù)336,求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 初始318分求調(diào)劑(有工作經(jīng)驗) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 生物學(xué)調(diào)劑招人!。 +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|