255
% Desenha PDF onde nenhum passou no teste
plot(xout_trecho(row_nenhum,:),n_pdf_trecho(row_nenhum,:),'dr','MarkerSize',8)
% ------------------------------------------------------------------------
% Desenha PDF onde somente RICE passou no teste
% subplot(2,2,1:2)
plot(xout_trecho(row_so_rice,:),n_pdf_trecho(row_so_rice,:),'dr','MarkerSize',8)
% ------------------------------------------------------------------------
% Desenha PDF onde RICE passou com menor Qui
subplot(2,2,1:2)
plot(xout_trecho(ind_min_rice,:),n_pdf_trecho(ind_min_rice,:),'dr','MarkerSize',8)
hold on
grid on
plot(r_trecho(:,ind_min_rice),p3_rice_trecho(:,ind_min_rice),'-b','LineWidth',3)
plot(r_trecho(:,ind_min_rice),p1_trecho(:,ind_min_rice),'-.g','LineWidth',3)
plot(r_trecho(:,ind_min_rice),p2_trecho(:,ind_min_rice),'-.m','LineWidth',3)
xlabel('valor de r=V/<V> (a)', 'FontSize', 12,'FontWeight', 'bold')
ylabel('PDF (num. observações)', 'FontSize', 12,'FontWeight', 'bold')
titulo=({'Comparação PDF: RICE, Gauss e Rayleigh';...
['Menor Qui - Setor: ' num2str(ind_min_rice) ...
' e Fator K=' num2str(k_rice_trecho(ind_min_rice))]});
title(titulo, 'FontSize', 10,'FontWeight', 'bold');
legend('Empírico','Rice', 'Gauss','Rayleigh');
% figure
% ------------------------------------------------------------------------
% Desenha PDF onde RICE passou com maior Qui
subplot(2,2,3:4)
plot(xout_trecho(ind_max_rice,:),n_pdf_trecho(ind_max_rice,:),'dr','MarkerSize',8)
hold on
grid on
plot(r_trecho(:,ind_max_rice),p3_rice_trecho(:,ind_max_rice),'-b','LineWidth',3)
plot(r_trecho(:,ind_max_rice),p1_trecho(:,ind_max_rice),'-.g','LineWidth',3)
plot(r_trecho(:,ind_max_rice),p2_trecho(:,ind_max_rice),'-.m','LineWidth',3)
xlabel('valor de r=V/<V> (b)', 'FontSize', 12,'FontWeight', 'bold')
ylabel('PDF (num. observações)', 'FontSize', 12,'FontWeight', 'bold')
titulo=({'Comparação PDF: RICE, Gauss e Rayleigh';...
['Maior Qui - Setor: ' num2str(ind_max_rice) ...
' e Fator K=' num2str(k_rice_trecho(ind_max_rice))]});
title(titulo, 'FontSize', 10,'FontWeight', 'bold');
legend('Empírico','Rice', 'Gauss','Rayleigh');
% ************************************************************************
% Fim
% ************************************************************************
% Nome do programa: m_taxa_de_cruzamento_de_nivel_rev_d
% ------------------------------------------------------------------------
% Este programa calcula a taxa de cruzamento de nível e a duração média do
% desvanecimento em cada setor de um trecho
% ************************************************************************
% Arquivo de entrada: dados_pdf_teste_qui_trecho_1.mat
% Variável(eis) de entrada:
% ------------------------------------------------------------------------
% Arquivo de saída1: .mat
% Variável(eis) de saída:
% Arquivo de saída2: .mat
% Variável(eis) de saída:
% ************************************************************************
clear all; % Apaga todas as variáveis do Espaço de Trabalho
clc; % Apaga área de comando
% clf; % Apaga gráficos
% ------------------------------------------------------------------------
% Variáveis de leitura do arquivo
% -----------------------------------------------------------------------
% caminho_entrada='C:\Users\Pedro\1-Projetos\Mestrado\Pesquisa de
caminho_entrada='C:\Users\Pedro\1-Projetos\Mestrado\Pesquisa de campo\Banda
Estreita\Medidas em Campo\Petrobras_3_5\Maritimo_Macae\PDF_desv_rapido_2\';
arquivo_entrada='dados_pdf_teste_qui_trecho_'; % Nome do arquivo de leitura