Выполните рациональную аппроксимацию, используя представление компонента с шестовым остатком
Используйте 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 | Написание Verilog-A описание 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 Mrad/s,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')
Схема сохранена в текущей папке.