clear
clc

Km=1;
Tm=1;
U=0;
Uon=0.1;
Uoff=0.05;
Aout=0;??%取值為0?1?-1

n=1;
t=0;
Dt=0.001;

for?i=1:10000
????Ac=sin（t）;
????Du=-U/Tm+Km/Tm*（Ac-Aout）;
????U=U+Du*Dt;
????if?Aout==-1
????????if?U>-Uoff
????????????Aout=0;
????????else
????????????Aout=-1;
????????end
????elseif?Aout==0
????????if?U>Uon
????????????Aout=1;
????????elseif?U>-Uon
????????????Aout=0;
????????else
????????????Aout=-1;
????????end
????else%if?Aout==1?
????????if?U????????????Aout=0;
????????else
????????????Aout=1;
????????end
????end

????Ac_store（n）=Ac;
????Aout_store（n）=Aout;
????U_store（n）=U;
????n=n+1;
????t=t+Dt;
end

figure（1）
plot（（1:n-1）*DtAc_store（1:n-1）*DtAout_store）
axis（[1?10?-1.2?1.2]）
legend（‘原指令‘‘PWPF調(diào)制后的指令‘）
grid?on

figure（2）
plot（U_storeAout_store）
axis（[-0.2?0.2?-1.2?1.2]）
xlabel（‘U‘）
ylabel（‘Aout‘）
grid?on

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件?????221863??2010-07-07?20:48??PWPF仿真.pdf

?????文件????????893??2010-07-07?19:52??PWPF.m

?????文件??????21754??2010-07-08?15:31??PWPF.mdl

-----------?---------??----------?-----??----

???????????????244510????????????????????3