Сообщество Engee

Моделирование LC-осциллятора на основе схемы Колпитса

Автор
avatar-mikhailpetrovmikhailpetrov
Notebook

Моделирование LC-осциллятора на основе схемы Колпитса

В этом примере представлена реализация схемы осциллятора Колпитса с номинальной частотой 9 МГц. Частота колебаний задается формулой:

LC-осцилляторы обладают хорошей частотной избирательностью благодаря более высоким уровням добротности, по сравнению с RC-осцилляторами.

Схема модели:

lc_oscillator_1734963351199.png

Определение функции для загрузки и запуска модели:

function start_model_engee()
    try
        engee.close("lc_oscillator", force=true) # закрытие модели 
        catch err # в случае, если нет модели, которую нужно закрыть и engee.close() не выполняется, то будет выполнена её загрузка после catch
            m = engee.load("$(@__DIR__)/lc_oscillator.engee") # загрузка модели
        end;

    try
        engee.run(m) # запуск модели
        catch err # в случае, если модель не загружена и engee.run() не выполняется, то будут выполнены две нижние строки после catch
            m = engee.load("$(@__DIR__)/lc_oscillator.engee") # загрузка модели
            engee.run(m) # запуск модели
        end
end
start_model_engee (generic function with 1 method)

Запуск симуляции

try
    start_model_engee() # запуск симуляции с помощью специальной функции, реализованной выше
    catch err
    end;

Запись из simout в переменные сигнала напряжения:

res = collect(simout)
V = collect(res[1])
7001×2 DataFrame
6976 rows omitted
Rowtimevalue
Float64Float64
10.04.08249
21.0e-811.9085
32.0e-89.00624
43.0e-85.58947
54.0e-86.79576
65.0e-87.37572
76.0e-87.26239
87.0e-86.60988
98.0e-85.7106
109.0e-85.90388
111.0e-78.27166
121.1e-79.73887
131.2e-79.79918
69906.989e-59.92202
69916.99e-59.48484
69926.991e-59.03062
69936.992e-58.69916
69946.993e-58.59152
69956.994e-58.73975
69966.995e-59.09795
69976.996e-59.55646
69986.997e-59.97573
69996.998e-510.2294
70006.999e-510.2408
70017.0e-510.0048

Визуализация результатов моделирования

using Plots
plot(V[:,1], V[:,2], label="Напряжение, В", linewidth=2, title="Напряжение LC-осциллятора")
plot(V[:,1], V[:,2], label="Напряжение, В", linewidth=2, xlim=(4e-5, 4.1e-5), title="Напряжение LC-осциллятора в течение 1 микросекунды")