bmarechal / sicav_sim

Blame view

PM/pm_vs_besselj.m 248 Bytes
edit raw normal view history
01f780010   bmarechal   rename ram to pha... 
1
2
3
4
5
6
7
8
9
10
 
  clc
  clear all
  close all
  
  N = 10;
  
  set(0, 'defaultlinelinewidth', 5);
  
  phi = [0:0.01:1]*2*pi;
  m = 1;
73d57d739   bmarechal   start from {J0,J1} 
11
 
  for n = 1:N
01f780010   bmarechal   rename ram to pha... 
12
13
14
15
16
 
  	freqs = -n:n;
  	amp = besselj(freqs, m);
  	plot(phi, abs(sum(amp'.*exp(1i.*phi.*freqs'))))
  	hold on
  end
73d57d739   bmarechal   start from {J0,J1} 
17
 
  legend(strsplit(num2str(1:N)))