У меня половина синуса со временем 0:2*T:
Rc = 1e3;
T = 1/Rc;
Fs = 2e3; % sampling frequency
dt = 1/Fs;
over = Fs/Rc; % sampling factor - 2
sps = 10;
time = 0:dt/sps:2*T;
half_Sine = sin(pi*time/(2*T)).^3;
figure(1);
plot(time,half_Sine, 'b--o');
grid on
xlabel('time','FontSize',13);
ylabel('a(t)','FontSize',13);
Но мне нужно время -T/2<= time<= T/2. И представить ось времени как time/T. Когда я делаю
time = -T/2:dt/sps:T/2;
Это дает мне не половину синуса. Итак, мне нужно что-то вроде этого:
1.- куб на функции sin
предотвращает симметрию результирующего графика по оси Y.
Rc = 1e3;
T = 1/Rc;
Fs = 2e3; % sampling frequency
dt = 1/Fs;
over = Fs/Rc; % sampling factor - 2
sps = 10;
t =-2*T :dt/sps:2*T;
y= sin(pi*t/(2*T)).^3;
figure;
plot(t,y, 'b--o');
grid on
xlabel('time','FontSize',13);
ylabel('a(t)','FontSize',13);
2.- Чтобы иметь максимум на t=0
, вам нужно использовать функцию cos
, а не sin
,
и квадрат, а не куб
Rc = 1e3;
T = 1/Rc;
Fs = 2e3; % sampling frequency
dt = 1/Fs;
over = Fs/Rc; % sampling factor - 2
sps = 10;
t =-T :dt/sps:T;
y= cos(pi*t/(2*T)).^2;
figure(1);
plot(t,y, 'b--o');
grid on
xlabel('time','FontSize',13);
ylabel('a(t)','FontSize',13);
Теперь у вас есть [-T T]
сюжет,
3.- нужный вам интервал [-T/2 T/2]
t = -T/2:dt/sps:T/2; y= cos(pit/(2T)).^2;
фигура; сюжет(t,y, 'b--o'); сетка на xlabel('время','Размер шрифта',13); ylabel('a(t)','Размер шрифта',13);
4.- Вы упомянули, что хотите нормализовать ось времени.
Если вы измените t делением на T, результирующий график будет очень узким временным интервалом около t = 0 и почти постоянным y = 1.
Вместо этого просто измените аннотацию оси x следующим образом.
figure;
hp1=plot(t,y, 'b--o');
hp1.XData=hp1.XData/T
grid on
xlabel('time/T','FontSize',13);
ylabel('a(t)','FontSize',13);