how to plot equation containing bessel functions
3 views (last 30 days)
Show older comments
Hello, I am struggling on this equation since some days. because I followed the theoretical equation of plate response but the plot still not correct. Please could you have a look on my script.
clc clear
%propriétés de la plaque
% épaisseur
h=0.001;
% fréquence
f = 100; %Hz
% nombre d'onde
k = 1000;
%%%%%%%%%%%%%
v0 =4 ;
n = 1000;
w = 2*pi*f;
% vecteur position
% longueur
lx=10;
% nb d'échantillons
Nx=512;
dx=lx/Nx;
rx=dx:dx:lx-dx;
% champ vx
J0 = besselj(0,k*rx); %fonction de bessel d ordre 0
J0i = besselj(0,1i*k*rx);
N0 = besseli(0,k*rx); %fonction de bessel d ordre 0 modifiée
N0i= besseli(0,1i*k*rx);
t1=0;
t2=0.023;
t3=2*t2;
v1t= v0*(J0+1i*N0-1i*J0i-1i*N0i)*exp(-1i*2*pi*f*t1);
v2t= v0*(J0+1i*N0-1i*J0i-1i*N0i)*exp(-1i*2*pi*f*t2);
v3t= v0*(J0+1i*N0-1i*J0i-1i*N0i)*exp(-1i*2*pi*f*t3);
figure(1)
plot(rx,real(v1t),'r',rx,real(v2t),'b',rx,real(v3t),'g',rx,real(v3t)*0,'k');
break;
Thank you in advance.
0 Comments
Accepted Answer
More Answers (2)
See Also
Categories
Find more on Bessel functions in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!