how to plot equation containing bessel functions

4 visualizzazioni (ultimi 30 giorni)
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.

Risposta accettata

Youssef  Khmou
Youssef Khmou il 8 Apr 2013
hi,
change lx to :
lx=0.7;
fine?

Più risposte (2)

philippe
philippe il 9 Apr 2013
hi,
Thank you for the answer.
But the problem is that the three functions shouldn't cross in zero at the same time because they are shifted between the three times. this is how the plot should be.
Thanks

philippe
philippe il 10 Apr 2013
any one has tried it before ?

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by