NumericalEvaluationEvalution.zip

% 均值和方差的对比 clear close all HRTF_LPB = zeros(200, 3); HRTF_multi = zeros(200, 3); load('HRTF_total.mat'); % 储存了数据 HRTF_LPB(:, 1:3) = HRTFeq_0(:, 2:4); % HRTFeq_0的2-4列为同一个人的3次测量 HRTF_LPB_log = 20*log10(HRTF_LPB)-20*log10(HRTF_LPB(1,:)); HRTF_LPB_log_ave = mean(HRTF_LPB_log, 2); HRTF_LPB_log_var = var(HRTF_LPB_log,0,2); LPB_var = sum(HRTF_LPB_log_var(1:180))/length(HRTF_LPB_log_ave(1:180)); HRTF_multi(:, 1) = HRTFeq_0(:, 5); load('ZYJ_YJR_0_L.mat'); HRTF_multi(:, 2) = P_tf_fre_sm(:, 1); %两列分别为两个人的测量 HRTF_multi(:, 3) = P_tf_fre_sm(:, 2); HRTF_multi_log = 20*log10(HRTF_multi)-20*log10(HRTF_multi(1,:)); HRTF_multi_log_ave = mean(HRTF_multi_log, 2); HRTF_multi_log_var = var(HRTF_multi_log,0,2); multi_var = sum(HRTF_multi_log_var(1:180))/length(HRTF_multi_log_ave(1:180)); %频率索引 N_fft = 400; %fft的点数，大于裁剪后信号的点数 f = linspace(0, 48000/2, N_fft/2+1); f = f(2:end); [~, freq_idx_3k] = ismember(3000, f); [~, freq_idx_8k] = ismember(8040, f); [~, freq_idx_18k] = ismember(18000, f); var_multi_0_3k = sum(HRTF_multi_log_var(1:freq_idx_3k))/length(HRTF_multi_log_var(1:freq_idx_3k)); var_LPB_0_3k = sum(HRTF_LPB_log_var(1:freq_idx_3k))/length(HRTF_LPB_log_var(1:freq_idx_3k)); var_multi_3k_8k = sum(HRTF_multi_log_var(freq_idx_3k:freq_idx_8k))/length(HRTF_multi_log_var(freq_idx_3k:freq_idx_8k)); var_LPB_3k_8k = sum(HRTF_LPB_log_var(freq_idx_3k:freq_idx_8k))/length(HRTF_LPB_log_var(freq_idx_3k:freq_idx_8k)); var_multi_8k_18k = sum(HRTF_multi_log_var(freq_idx_8k:freq_idx_18k))/length(HRTF_multi_log_var(freq_idx_8k:freq_idx_18k)); var_LPB_8k_18k = sum(HRTF_LPB_log_var(freq_idx_8k:freq_idx_18k))/length(HRTF_LPB_log_var(freq_idx_8k:freq_idx_18k)); fprintf('0-3k Hz : multi:%f; single:%f \n', var_multi_0_3k, var_LPB_0_3k); fprintf('3k-8k Hz : multi:%f; single:%f \n', var_multi_3k_8k, var_LPB_3k_8k); fprintf('8k-18k Hz : multi:%f; single:%f \n', var_multi_8k_18k, var_LPB_8k_18k); fprintf('0-18k Hz : multi:%f; single:%f \n', multi_var, LPB_var ); %% figure semilogx(f, HRTF_LPB_log); xlim([100,18000]); xlabel('Frequency/Hz') ylabel('Magnitude/dB') legend ('LPB_1', 'LPB_2', 'LPB_3') figure semilogx(f, HRTF_multi_log); xlim([100,18000]); xlabel('Frequency/Hz') ylabel('Magnitude/dB') legend ('LPB', 'ZYJ', 'YJR') figure; semilogx(f, HRTF_LPB_log_var); hold on semilogx(f, HRTF_multi_log_var); hold off xlabel('Frequency/Hz') ylabel('Magnitude/dB') legend ('single', 'multi') title('Var') xlim([100,18000]) figure; semilogx(f, HRTF_LPB_log_ave); hold on semilogx(f, HRTF_multi_log_ave); hold off xlabel('Frequency/Hz') ylabel('Magnitude/dB') legend ('single', 'multi') xlim([100,18000]) title('Ave')