| 5 | 1/1 | 返回列表 |
| 查看: 1711 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
xiaojixiong金蟲 (小有名氣)
|
[求助]
求助:如何求含貝塞爾函數(shù)的超越方程的復(fù)數(shù)解
|
|
|
如題,小弟最近碰到一含貝塞爾函數(shù)的超越方程。程序如下: i=1; for x=0:0.001:5 % x=kc/Wp for y=0:0.001:5 % y=w/Wp if (y-x)<=4.2&&0<(y-x) t=0.6*1000*sqrt(y.*y-x.*x); ta=0.6*sqrt(y.*y-x.*x); tc=0.6*2* sqrt(y.*y-x.*x); t1=(y.*y-1); t2=x.*x*(3.4596-y.*y); t3=y.*y*(y*y-4.4596); T=0.6*1000*sqrt((t1.*(t2+t3))/t3); Ta=0.6*sqrt((t1.*(t2+t3))/t3); pTa=(Ta/4).*(Ta/4); pta=(ta/4).*(ta/4); ptc=(tc/4).*(tc/4); % ta有關(guān)的貝塞爾函數(shù)(4個); cta0=besselj(0,ta(:)); cta1=besselj(1,ta(:)); dta0=bessely(0,ta(:)); dta1=bessely(1,ta(:)); % tc 有關(guān)的貝塞爾函數(shù)(2個) ctc0=besselj(0,tc(:)); dtc0=bessely(0,tc(:)); A=ctc0.*dta0-cta0.*dtc0; B=-ctc0.*dta1+cta1.*dtc0; cTa0=besselj(0,Ta(:)); cTa1=besselj(1,Ta(:)); f=((y.*y-1)/T).*cTa1.*A+(y.*y/t)*cTa0 *B; if abs(f)<=0.00000001 a0(i)=x;b0(i)=y; i=i+1; disp([num2str(x) ' ' num2str(y)]); end end end end 我通過限制f絕對值的精度,去求出了滿足f=0的實數(shù)解,但是我若想求出滿足f=0,x和對應(yīng)的y的復(fù)數(shù)解。這個程序應(yīng)該怎么改啊,謝謝大家了。 |
金蟲 (小有名氣)
銅蟲 (著名寫手)
金蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 348求調(diào)劑 +3 | 小懶蟲不懶了 2026-03-28 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 0856求調(diào)劑 +6 | 楒桉 2026-03-28 | 6/300 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
[考研]
|
y7czhao 2026-03-26 | 10/500 |
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機用戶 2026-03-26 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我!。 +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 307求調(diào)劑 +8 | 超級伊昂大王 2026-03-24 | 9/450 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 329求調(diào)劑 +7 | 鈕恩雪 2026-03-25 | 7/350 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 10/500 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
[考研] 材料/農(nóng)業(yè)專業(yè),07/08開頭均可,過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|