Signal super gaussienne

[invitec8026e9d]

Bonjour tout le monde ;

Voila mon but dans ce programme est d'introduire a l'entrée une gaussienne et je veux avoir en sortie le méme signal ( j'utilise pour cela Matlab ) quand je simule ya une erreure si une personne peut m'aidée ça serrait sympa voici mon algorithme :

 Po= 4; % [mW] initial peak power of signal source 
C = 0; % Chirping Parameter 
m = 1; % Super Gaussian parameter (m=1 ==> Gaussian) 
t0 = 100e-12; %[sec] initial pulse width 
vgroup= 2e8; 
alpham=0.24; 
attn=0.24; 
taum = 10000; 
dtau = 2*taum/2^12; %1.953 
tunit= 1e-12; % make time unit in psec 
tau = (-1000:dtau:(taum-dtau))*tunit; % 
% fs = 1/(dtau*tunit);

at = sqrt(Po)*exp(-0.5*(1+i*C)*(tau/t0).^(2*m)); % Input field in the time domain 
subplot(2,2,1); 
plot(tau,at); 
title('input wave');

for x=0:16;
    N=length(x);

bt = sqrt(Po)*exp(-0.5*((tau-((2*x)/vgroup))/t0).^2)*exp(3.84);

subplot(2,2,2) 
plot(N,bt,'-rx'); 
end 
l=trapz(x,bt); 
subplot(2,2,3); 
plot(t,l,'-rx');

Merci .
-----

[Antoane]

Bonjour et bienvenue sur Futura!

Les messages d'erreur donnés par Matlab sont généralmenet assez clairs, il faut les lire, et essayer de comprendre. Ici, Matlab te dit:

Error using trapz (line 58)
LENGTH(X) must equal the length of Y in dim 2.

Y, c'est ton vcteur ligne bt.
Une solution (à voir si c'est effectivement ce que tu veux faire physiquement), consiste à écrire :
l=trapz(x*ones(size(bt)), bt) ;

Une fois que tu auras réglé ce problème, il faudra t'intéresser au fait que t n'est pas défini.

Et, accessoirement, remplace i par 1i dans :
at = sqrt(Po)*exp(-0.5*(1+i*C)*(tau/t0).^(2*m)); % Input field in the time domain

et il serait peut-être nécessaire d'ajouter un "hold on" après subplot(222) (qui pourrait d'ailleurs sortir de la boucle), sans quoi seul le dernier plot (i.e. celui pour x=16) sera conservé sur le subplot.

Enfin, pense à bien indenter ton programme.