%% Historico aula 9
cd aula_9\
what
addpath c:\Users\zwinglio\Documents\MATLAB\TCABR\ -end
which coisas_em_janelas.m
which coisas_em_janelas.m -all
which leia_TCABR.m -all
edit Visualiza_TCABR_com_fit.m
help uigetfile
ind_ultimo_under = find( FileName=='_', 1, 'last' )
ind_ultimo_dot = find( FileName=='.', 1, 'last' )
FileName
DISP = str2num( FileName( (ind_ultimo_under+1):(ind_ultimo_dot-1) ) )
FilterSpec = 'IPlasma_*.dat';
DialogTitle = 'Indique o disparo a ser lido';
DefaultName = 'c:\Users\zwinglio\Documents\MATLAB\TCABR\IPlasma_25285.dat';
[FileName,PathName,FilterIndex] = ...
uigetfile(FilterSpec,DialogTitle,DefaultName)
visualiza_MHD
visualiza_MHD(25285, 2 )
5
help inputdlg
help fminsearch
edit coisas_em_janelas.m
[tM,yM] = leia_TCABR( 25285, 'MirR' );
[tM,yM] = leia_TCABR( 25285, 'MRef' );
clear
[tM,yM] = leia_TCABR( 25285, 'MRef' );
figure, plot( tM, yM )
help coisas_em_janelas
[Yj,ind] = coisas_em_janelas(yM,[250], @(x)std(x) );
figure, plot( Yj )
figure, plot( mean(tM(ind),2), Yj, '.' )
tok = tM( tM>80 & tM<90 );
yok = yM( tM>80 & tM<90 );
[Yj,ind] = coisas_em_janelas(yok,[250], @(x)ajusta_MHD((1:length(x))*4e-3,x) );
[Yj,ind] = coisas_em_janelas(yok,[250], @(x)ajusta_MHD((1:length(x))'*4e-3,x) );
Yj
tok = tM( tM>80 & tM<100 );
yok = yM( tM>80 & tM<100 );
[Yj,ind] = coisas_em_janelas(yok,[250], @(x)ajusta_MHD((1:length(x))'*4e-3,x) );
yok
Yj
figure, plot( 1./Yj )
visualiza_MHD(25285, 1 )
help fft
clear
home
[tM,yM] = leia_TCABR( 25285, 'MRef' );
tok = tM( tM>80 & tM<90 );
yok = yM( tM>80 & tM<90 );
help my_fft
[Y,f] = my_fft( yok, tok );
figure, plot( f, Y )
tok = tM( tM>80 & tM<81 );
yok = yM( tM>80 & tM<81 );
[Y,f] = my_fft( yok, tok );
figure, plot( f, Y )
[Yj,ind] = coisas_em_janelas(yok,[250], @(x)my_fft(x,(1:length(x))*4e-4) );
Yj
tok = tM( tM>80 & tM<100 );
yok = yM( tM>80 & tM<100 );
[Yj,ind] = coisas_em_janelas(yok,[250], @(x)my_fft(x,(1:length(x))*4e-4) );
whos
figure, imagesc( Yj )
figure, imagesc( Yj' )
set( gca, 'Ydir', 'Normal' )
help spectrograma
[ Yj, f, tj ] = spectrograma( yM, tM, [250 125] );
figure, imagesc( tj, f, Yj' ), set( gca, 'Ydir', 'Normal' )
tok = tM( tM>50 & tM<120 );
yok = yM( tM>50 & tM<120 );
[Tj,ind] = coisas_em_janelas(yok,[250], @(x)ajusta_MHD((1:length(x))'*4e-3,x) );
hold on
plot( mean(tok(ind),2), 1./Tj, '.-w' )
home
clear
home
what
edit mmq_lin
edit mmq_naolin.m
