function y=squarewav(Vm,f,K,n)      %参数为幅度，频率,占空比，以及波形个数

% 正弦信号
% f = 1000;% 正弦波频率
w = 2*pi*f;% 正弦波角频率
% Vm = 2;% 正弦幅值
u = (0.5-K)*pi;% 相位
N = 4000;% 调节采样率，采样率越大，波形误差越小
M = n*N;
Fs = N*f;% 采样率：1/f_sameple为采样时间间隔，通常以N倍的信号源来表示
t = [0:M]/Fs;% 采样时间向量(始终采N个点，即刚好采信号源的一个周期)
y_sin = Vm*sin(w*t+u);

for i=1:M+1
    if(y_sin(i)>=y_sin(1))
        y_plus(i) = Vm;
    else
        y_plus(i) = -Vm;
    end
end

figure
subplot(2,1,1)
plot(t,y_sin,t,y_plus,'r',t,0,'-')
subplot(2,1,2)
ylim([-2.5 2.5])
plot(t,y_plus,'r',t,0,'-')
ylim([-2.5 2.5])