| 5 | 1/1 | 返回列表 |
| 查看: 1761 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
cologate007新蟲 (初入文壇)
|
[求助]
matlab對(duì)含貝塞爾函數(shù)的超越方程求根 已有1人參與
|
|
|
有一個(gè)含有貝塞爾函數(shù)的超越方程,用solve函數(shù)求解,輸入matlab直接給出了一個(gè)復(fù)數(shù)結(jié)果,但是實(shí)際上超越方程是有很多個(gè)根的,不知道給出的這個(gè)解是第幾階的?怎么樣才能得到它的第n(任意的)階根呢?采用牛頓迭代法給定不同的初值可以計(jì)算得到不同的根,但是也不能確定是第幾階。不知道論壇里有沒(méi)有朋友有更好的計(jì)算方法?函數(shù)就不貼了,下面是牛頓迭代法的matlab程序,函數(shù)見里面的f。請(qǐng)了解的朋友給點(diǎn)意見,感謝! 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ā)表 | |
|---|---|---|---|---|
|
[考研] 283求調(diào)劑(080500) +7 | A child 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
[考研] 求調(diào)劑 +9 | 張zz111 2026-03-27 | 10/500 |
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
[考研] 求調(diào)劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問(wèn)津 2026-03-23 | 5/250 |
|