Что можно сделать чтоб ФЧХ совпали?(

clc
close all
clear
syms w
w1 = logspace(-5,1,1000);
den = [1 6 5 8];
num = [1];
sys = tf(num, den)
 
[comp2, zcomp2, Kcomp2, prost2, Kprost2] = pa3del(den);
[fchx2 ] = func6(comp2, zcomp2, prost2);
 
 
figure
grid on 
hold on
phase = unwrap(angle(h));
plot(w1,phase*180/pi,'black','linewidth',3);
title('Фазо-частотная характеристика')
fchx22=subs(fchx2, w, w1, 0);
plot(w1,(fchx22)*180/pi,'red')

function [ X1  ] = func6(comp1, comp2, prost1)
comp1
comp2
prost1
x1=0;
x2=0;
syms w 
for i = 1:length(prost1)
    F = -atan(w*prost1(i));
    x1 = F+x1;
end
 
for i = 1:length(comp1)
     F = -atan((2*comp1(i)*comp2(i)*w)/(1-comp1(i)^2*w^2));
     x2 = x2+ F;
end
X1 = x1+x2;
end

function [comp1, comp2, Kcomp, prost1, Kprost] = pa3del(X)
syms p
comp = [];
prost = [];
comp1 = [];
comp2 = [];
prost1 = [];
Kprost = 1;
Kcomp = 1;
korni = round(roots(X)*1e6)/1e6
for i = 1:length(korni)
    if ~isreal(korni(i)) % если корень комплексный
        comp = [comp korni(i)]; % добавляем в комплексные корни
    else % остались простые корни
        prost = [prost korni(i)];  % добавляем в простые корни
    end
end
for i = 1:2:length(comp)
    Re = real(comp(i));
    Im = imag(comp(i));
    ReIm = Re*Re+Im*Im;
    Kcomp = Kcomp*ReIm; 
    sdd1 = sqrt(1/ReIm);
    sdd2 = -sdd1*Re;
    comp1 = [comp1 sdd1];
    comp2 = [comp2 sdd2];
end
for i = 1:length(prost)
    Kprost = Kprost*abs(prost(i));
    prost1 =[prost1 1/abs(prost(i))];
end
end

        for j=1:length(Fiw3)
            if ((Fiw3(j))>0);
            Fiw4(j)=Fiw3(j)-pi;
            else
            Fiw4(j)=Fiw3(j);
            end
        end