Mi código solo escribe el nombre del archivo
Publicado por Fredy (8 intervenciones) el 17/06/2013 01:02:50
Hola ya había posteado este código pero como el problema ya no es un error creo que no es repost de lo contrario perdón.
Ejecutando el código pareciera estar todo bien "no genera ningún error" y en teoría por cada valor de P y An toman un valor y se asignan a v1 y v2 pero al final solo arroja el nombre del archivo "epotencia" a que se debe esto?. este es el código.
Ejecutando el código pareciera estar todo bien "no genera ningún error" y en teoría por cada valor de P y An toman un valor y se asignan a v1 y v2 pero al final solo arroja el nombre del archivo "epotencia" a que se debe esto?. este es el código.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
clc;
clear;
[F,Fs]=wavread('u_hombre');
wavplay(F,Fs);
t=size(F,1)/Fs;
n=0:1/Fs:t;
plot(n(2:end),F)
xlim([0 t])
wo=(2*pi)/t
%aproximadamente un ciclo
%T=size(nboch,1)/Fs
for n=0.005283:0.0001:t
ao=0;a(o)=ao +(1/t)*trapz(t,F);
a(n)=0;
u=cos(n*wo*t)
v=sin(n*wo*t)
a(n)= a(n)+(2/t)*trapz(n,F.*cos(n*wo*t));
b(n)=0;
b(n)= b(n)+(2/t)*trapz(n,F.*sin(n*wo*t));
%A[n] toma el valor de a[n] y b[n] en cada ciclo
An=sqrt(a(n)^2+b(n)^2); %calculo la potencia
P=Ao^2+(a(n)^2/sqrt(2))^2
if a(n)>0&&b(n)>0
tetan=atan(-b(n)/a(n));
else if a(n)<0&&b(n)>0
tetan=90+atan(-b(n)/a(n));
else if a(n)<0&&b(n)<0
tetan=180+atan(-b(n)/a(n));
else if a(n)>0&&b(n)<0
tetan=270+atan(-b(n)/a(n));%calculo el angulo
%ingreso el valor de cada A[n] y P en un vector v
v1(n)=P;
v2(n)=An;
T(n)=n
end
plot(T,v1)
plot(T,v2)
Valora esta pregunta
0
