%MATLAB實現壓控振蕩器

clear all;
clc;
t0 = 0.15;                          %定義壓控信號持續時間
ts = 0.0001;                        %定義信號采樣率
fc = 50;                            %定義振蕩頻率
t = [0:ts:t0];                      %時間矢量
u0 = 20*t(1:length(t)/3);           %定義壓控信號（單周期）
u = [u0,u0,u0,0];                   %定義壓控信號（3個周期）
Ac = 1;                             %定義振幅
kc = 0.1;                           %定義輸入信號靈敏度
fi = pi/4;                             %定義初始相位
%對壓控信號進行積分
u_int(1) = 0;                       %定義壓控信號積分初值
for i = 1:length(u)-1               %進行離散積分
    u_int(i+1) = u(i)+u_int(i);
end
% nn=u_int(540)
%壓控振蕩信號
y = Ac*cos(2*pi*fc*t+2*pi*kc*u_int+fi);
% y = Ac*cos(2*pi*fc*t+2*pi*kc*u_int(i+1)+fi);
% nn= Ac*cos(2*pi*fc*t+2*pi*kc*u_int(10)+fi);
%作圖
subplot(3,1,1)                      %顯示壓控信號
plot(t,u)
axis([0 t0 -1.5 1.5])
title('控制信號')
subplot(3,1,2)                      %顯示振蕩信號
plot(t,y)
axis([0 t0 -1.5 1.5])
title('壓控振蕩信號')
% subplot(3,1,3)  
% plot(t,nn)
% axis([0 t0 -1.5 1.5])