Matlab Восстанавливая freqz, нормализуя ось x и получая половину графика

Question

У меня есть функция, которая в основном воссоздает команду freqz в matlab. Я выяснил, как построить всю трансформацию моего частотного отклика, но мне нужно только половину его, и мне нужно нормализовать его от pi до 1 (где 0: pi представляет мою ось x, и я хочу, чтобы это 0: 1). Код здесь:

function freqrespplot(b, a) 

hb = fft(b,512); 
ha = fft(a,512); 

magh = (abs(hb) ./ abs(ha)); 
angh = angle(hb ./ ha); 

x = linspace(0,2*pi, 512); 

subplot(2,1,1); 
plot(x,magh,'g'); 
subplot(2,1,2); 
plot(x,angh,'r'); 

end 

В примере б = [1 2 2], а = [0 1 0,8], участки, такие как следующие:

freqrespplot ([1 2 2], [0 1 .8]); Магнитуда

и мой соответствующий фазовый участок

Я хочу, чтобы ось х (омега), чтобы перейти от 0 до 1, в обоих случаях, и соответствуют 0 до пи, принимая только половину граф, как это возможно, если:

Answer 1

Вам нужно только некоторые небольшие изменения, отмеченные с комментариями в коде ниже:

function freqrespplot(b, a) 

hb = fft(b,512); 
ha = fft(a,512); 

magh = (abs(hb) ./ abs(ha)); 
angh = angle(hb ./ ha); 

x = linspace(0,2, 513); %// 2, not 2*pi. And 513. Last value will be discarded 
x = x(1:end-1); %// discard last value to avoid having 0 and 2*pi (they are the same) 

subplot(2,1,1); 
plot(x(1:end/2),magh(1:end/2),'g'); %// plot only half 
subplot(2,1,2); 
plot(x(1:end/2),angh(1:end/2),'r'); %// plot only half 

end 

С вашим примером b и a это производит

Вы можете включить один дополнительный образец для достижения правого края графика. Я комментирую только различие с выше кодом:

function freqrespplot(b, a) 

hb = fft(b,512); 
ha = fft(a,512); 

magh = (abs(hb) ./ abs(ha)); 
angh = angle(hb ./ ha); 

x = linspace(0,2, 513); 
x = x(1:end-1); 

subplot(2,1,1); 
plot(x(1:end/2+1),magh(1:end/2+1),'g'); %// plot only lower half plus one value 
subplot(2,1,2); 
plot(x(1:end/2+1),angh(1:end/2+1),'r'); %// plot only lower half plus one value 

end 

Сравните полученный график (правая часть):