Создайте передаточную функцию

Считайте в 4-портовой объединительной плате данные S-параметра от 'default.s4p'.

S = sparameters('default.s4p');
fourportdata = S.Parameters;
freq = S.Frequencies;
fourportZ0 = S.Impedance;

Преобразуйте 4-портовые одноконцевые S-параметры в 2-портовые дифференциальные S-параметры

diffdata = s2sdd(fourportdata);
diffZ0 = 2*fourportZ0;

Создайте передаточную функцию из дифференциальных 2-портовых данных

tfdata = s2tf(diffdata,diffZ0,diffZ0,diffZ0);

Анализ выхода рационализации при использовании значения по умолчанию для 'DelayFactor'

Используйте freqresp функция для вычисления отклика выхода rationalfit.

defaultfit = rationalfit(freq,tfdata)

Warning: Achieved only -10.2 dB accuracy with 48 poles, not -40.0 dB.  Consider specifying a larger number of poles using the 'NPoles' parameter.

defaultfit = 
   rfmodel.rational with properties:

        A: [48x1 double]
        C: [48x1 double]
        D: 0
    Delay: 0
     Name: 'Rational Function'

respfreq = 0:4e6:20e9;
defaultresp = freqresp(defaultfit,respfreq);

Обратите внимание, что параметр 'Delay' равен нулю (без задержки, удаленной из данных).

Постройте график исходных данных по сравнению с выходами по умолчанию rationalfit.

figure
subplot(2,1,1)
tfdataDB = 20*log10(abs(tfdata));
plot(freq,tfdataDB,'.-')
hold on
plot(respfreq,20*log10(abs(defaultresp)))
hold off
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')
defaultnpoles = numel(defaultfit.A);
defstr = ['Default DelayFactor (Uses ',num2str(defaultnpoles),' poles)'];
title(defstr)
legend('Original Data','Default rationalfit','Location','best')
subplot(2,1,2)
tfdataphase = 180*unwrap(angle(tfdata))/pi;
plot(freq,tfdataphase,'.-')
hold on
plot(respfreq,180*unwrap(angle(defaultresp))/pi)
hold off
xlabel('Frequency (Hz)')
ylabel('Angle (degrees)')
legend('Original Data','Default rationalfit','Location','best')

Обратите внимание, что результаты при использовании настроек по умолчанию rationalfit бедны. Поскольку фаза исходных данных имеет очень большой отрицательный наклон, может оказаться возможным улучшить точность рациональной функции при помощи параметра 'DelayFactor'.

Анализ выхода рационализации при использовании пользовательского значения для 'DelayFactor'

Значение 'DelayFactor' должно быть задано в диапазоне от 0 до 1. Выбор значения является упражнением в пробе и ошибке. Для некоторых наборов данных (тех, чья фаза имеет общий уклон вверх) изменение значения 'DelayFactor' не будет иметь никакого эффекта на результат.

Хранение всех других возможных параметров rationalfit константа, 0,98 найдено, чтобы создать хорошую подгонку.

customfit = rationalfit(freq,tfdata,'DelayFactor',0.98)

customfit = 
   rfmodel.rational with properties:

        A: [31x1 double]
        C: [31x1 double]
        D: 0
    Delay: 6.5521e-09
     Name: 'Rational Function'

customresp = freqresp(customfit,respfreq);

Обратите внимание, что параметр 'Delay' не равен нулю (rationalfit удалил некоторую задержку из данных).

Постройте график зависимости исходных данных от пользовательского выхода rationalfit.

subplot(2,1,1)
plot(freq,tfdataDB,'.-')
hold on
plot(respfreq,20*log10(abs(customresp)))
hold off
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')
customnpoles = numel(customfit.A);
customstr = ['DelayFactor = 0.98 (Uses ',num2str(customnpoles),' poles)'];
title(customstr)
legend('Original Data','Custom rationalfit','Location','best')
subplot(2,1,2)
plot(freq,tfdataphase,'.-')
hold on
plot(respfreq,180*unwrap(angle(customresp))/pi)
hold off
xlabel('Frequency (Hz)')
ylabel('Angle (degrees)')
legend('Original Data','Custom rationalfit','Location','best')

Рациональная функция, созданная при помощи пользовательского значения для 'DelayFactor', намного точнее и использует меньше полюсов.