MATLAB Answers

how to plot equation containing bessel functions

26 views (last 30 days)
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

Sign in to comment.

Accepted Answer

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

  1 Comment

philippe
philippe on 9 Apr 2013
hi, no it is just the factor of the samples number.

Sign in to comment.

More Answers (2)

philippe
philippe on 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

  0 Comments

Sign in to comment.


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

  0 Comments

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!