产生一个占空比呈正弦周期变化的方波

%%% 正弦波  对应的方波和三角波
clc
clear
close

%% 参数
fz = 600000;    %%采样频率
count = 30;     %%三角波个数及方波个数

%% 主程序
xL = [1/fz:(1/fz):1];
accuracy = fz/count;
x = 1:1:fz;

y_sin = sin(2*pi*x/fz);
y_sin = (0.8*y_sin)/2+0.5; %使正弦波幅值在0.1~0.9之间
flag = 1;
y_san(1)=0;

for i=1:fz
    y_san(i+1) = y_san(i)+flag * (1/accuracy); 
    if rem(i,accuracy) == 0
       flag = -flag;      
    end
end
y_san = y_san(1:fz);
for i = 1:fz
   if y_san(i)>y_sin(i)
      y_fang(i) = 0;
   else
      y_fang(i) = 1;
   end
end
figure(1)
plot(xL,y_san,'k')
hold on
plot(xL,y_sin,'k')
axis([0 1 -0.1 1.1])
%axis off
figure(2)
plot(xL,y_fang,'k')
hold on
plot(xL,y_sin,'k')
axis([0 1 -0.1 1.1])
%axis off