time_fun_s.m

我做了很久的

function y=time_fun_s(t)
syms a a1 
M=input('周期与脉冲宽度之比M=');
T=input('please input信号的周期T=')
A=1;tao=T/M;a=tao/2;
y1=sym('Heaviside(t+a1)')*A;
y=y1-sym('Heaviside(t-a1)')*A;
y=subs(y,a1,a);
y=simple(y);