Выполнить рациональную подгонку с использованием полюсно-остаточного представления компонента
Используйте rational класс для представления РЧ-компонентов с использованием объекта рациональной функции вида:
Существует два способа построения объекта рациональной функции:
Объект рациональной функции можно подогнать к данным компонента с помощью rationalfit функция.
Вы можете использовать rfmodel.rational для непосредственного задания представления «полюс-остаток» компонента.
h = rfmodel.rational возвращает объект рациональной функции, для свойств которого заданы значения по умолчанию.
h = rfmodel.rational('Property1',value1,'Property2',value2,...) задает свойства, используя одну или несколько пар имя-значение. Можно указать несколько пар имя-значение. Заключить имя каждого свойства в кавычку
timeresp | Временной отклик для рационального объекта и rationalfit объект функции |
stepresp | Реакция ступенчатого сигнала для рационального объекта и rationalfit объект функции |
freqresp | Частотная характеристика рационального объекта и rationalfit объект функции |
impulse | Импульсная характеристика для объекта рациональной функции |
ispassive | Вернуть значение true, если rationalfit выход пассивен на всех частотах |
makepassive | Обеспечение пассивности rationalfit выход или рациональный объект |
passivity | Пассивность графика N-by-N
rationalfit вывод функции |
pwlresp | Расчет временной характеристики кусочно-линейного входного сигнала |
generateSPICE | Создать файл SPICE из rationalfit S-параметров |
writeva | Write Verilog - описание rational объект |
Соответствие рациональной функции данным из rfdata.data объект.
S = sparameters('defaultbandpass.s2p');
freq = S.Frequencies;
data = rfparam(S,2,1);
fit = rationalfit(freq,data)fit =
rfmodel.rational with properties:
A: [10x1 double]
C: [10x1 double]
D: 0
Delay: 0
Name: 'Rational Function'
Постройте объект рациональной функции, rat, с поляками на уровне-4 мрад/с,-3 выпускниками/с, и-5 выпускниками/с и остатками 600 мрад/с, 2 выпускниками/с и 4 выпускниками/с.
rat=rfmodel.rational('A',[-5e9,-3e9,-4e6],'C',[6e8,2e9,4e9]);
Выполните анализ частотной области от 1,0 МГц до 3,0 ГГц.
f = [1e6:1.0e7:3e9];
Постройте график результирующей частотной характеристики в децибелах на плоскости X-Y.
[resp,freq]=freqresp(rat,f); figure plot(freq/1e9,20*log10(abs(resp))); xlabel('Frequency (GHz)') ylabel('Magnitude (dB)')
Чтение файла с именем passive.s2p и подогнать S-параметры 2 на 2. Создайте файл SPICE этих S-параметров.
S = sparameters('passive.s2p'); fit = rationalfit(S); generateSPICE(fit,'passive.ckt')
Цепь сохраняется в текущей папке.