不用混編,直接用_jn。
網(wǎng)站設(shè)計(jì)、成都網(wǎng)站設(shè)計(jì)服務(wù)團(tuán)隊(duì)是一支充滿著熱情的團(tuán)隊(duì),執(zhí)著、敏銳、追求更好,是創(chuàng)新互聯(lián)的標(biāo)準(zhǔn)與要求,同時(shí)竭誠為客戶提供服務(wù)是我們的理念。創(chuàng)新互聯(lián)建站把每個(gè)網(wǎng)站當(dāng)做一個(gè)產(chǎn)品來開發(fā),精雕細(xì)琢,追求一名工匠心中的細(xì)致,我們更用心!
double _jn( int n, double x );
這是Help:
_j0, _j1, _jn
These routines return Bessel functions of the first kind: orders 0, 1, and n, respectively.
Example
#include math.h
#include stdio.h
int main( void )
{
double x = 2.387;
int n = 3, c;
printf( "Bessel functions for x = %f:\n", x );
printf( " Kind Order Function Result\n\n" );
printf( " First 0 _j0( x ) %f\n", _j0( x ) );
printf( " First 1 _j1( x ) %f\n", _j1( x ) );
for( c = 2; c 5; c++ )
printf( " First %d _jn( %d, x ) %f\n", c, c, _jn( c, x ) );
printf( " Second 0 _y0( x ) %f\n", _y0( x ) );
printf( " Second 1 _y1( x ) %f\n", _y1( x ) );
for( c = 2; c 5; c++ )
printf( " Second %d _yn( %d, x ) %f\n", c, c, _yn( c, x ) );
}
語法不對(duì),不能用x作為索引號(hào)
x=0:0.2:10;
y00=(besselj(-1,x(:))-besselj(1,3*x(:)))/2;
y11=(bessely(-1,x(:))-bessely(1,3*x(:)))/2;
y01=(besselj(-1,x(:))-besselj(1,3*x(:)))/2;
y10=(besselj(-1,x(:))-besselj(1,3*x(:)))/2;
y=y00(:).*y11(:)-y01(:).*y10(:);
plot(x,y)
MATLAB提供了計(jì)算貝塞爾函數(shù)的函數(shù),具體包括:
besselj - 第一類貝塞爾函數(shù),或簡稱貝塞爾函數(shù);
bessely - 第二類貝塞爾函數(shù),又稱諾伊曼函數(shù)(Neumann function);
besseli - 第一類修正貝塞爾函數(shù);
besselk - 第二類修正貝塞爾函數(shù);
besselh - 第三類貝塞爾函數(shù),又稱漢克爾函數(shù)(Hankel function).
這幾個(gè)函數(shù)的調(diào)用語法基本相同,例如
J = besselj(nu,Z)
J = besselj(nu,Z,1)
[J,ierr] = besselj(nu,Z)
其中,nu為貝塞爾函數(shù)的階數(shù),Z為函數(shù)自變量.階數(shù)必須為實(shí)數(shù),但Z可以是復(fù)數(shù).
值得一提的是,上述函數(shù)是MATLAB基本模塊(也就是說不需要任何附加的工具箱)提供的特殊函數(shù),采用數(shù)值方法計(jì)算;而符號(hào)數(shù)學(xué)工具箱則提供了第一和第二類的4個(gè)貝塞爾函數(shù),名稱和調(diào)用方式都與MATLAB基本系統(tǒng)的4個(gè)函數(shù)完全一致,但支持微分、積分等符號(hào)運(yùn)算.
當(dāng)前名稱:貝塞爾函數(shù)c語言程序 貝塞爾函數(shù)的應(yīng)用
文章起源:http://jinyejixie.com/article48/dodhihp.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供微信小程序、網(wǎng)站排名、網(wǎng)站設(shè)計(jì)公司、網(wǎng)站策劃、App開發(fā)、用戶體驗(yàn)
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請(qǐng)盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場,如需處理請(qǐng)聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來源: 創(chuàng)新互聯(lián)