How can i transform rad/s to hertz for graph?

218 views (last 30 days)
clc
clear
hold on
fid = fopen('2.1.2.txt');
a=fscanf(fid,'%g %g %g %g',[4 inf]);
a = a';
fclose(fid)
%tm=a(:,1);
s1=a(:,2);
s2=a(:,3);
s3=a(:,4);
ts=0:(1/15000):(1/15000)*(length(s1)-1);
makts=3; %%%%%%%%%%%%%DEĞİŞTİR
s1=s1/max(s1);
s2=s2/max(s2);
s3=s3/max(s3);
subplot(3,2,1)
plot(ts,s1,'k')
axis([0 makts -1 1])
ylabel('Normalized Amplitude')
xlabel('Time (s)')
subplot(3,2,3)
plot(ts,s2,'k')
axis([0 makts -1 1])
ylabel('Normalized Amplitude')
xlabel('Time (s)')
subplot(3,2,5)
plot(ts,s3,'k')
axis([0 makts -1 1])
ylabel('Normalized Amplitude')
xlabel('Time (s)')
%%%%%%%%%DC%%%%%%
x=s1;
sp=fft(x);
nn=length(x);
sp(1:2)=zeros(size(1:2));
sp(nn-1:nn)=zeros(size(1:2));
x=real(ifft(sp));
%%%%%%%%%%%%%%%%%%
fss=15000;
T=1/fss;
upfreq=fss/2;
% ************ FFT ****************
p1=length(x);
z=(2/p1)*abs(fft(x));
z=z(1:p1/2-1);
z(1,1)=0.5*z(1,1);
nff=length(z);
nf1=round(2*nff*upfreq/fss);
z=z(1:nf1);
subplot(3,2,2)
s1z=abs(z)/max(abs(z));
plot((0:nf1-1).*(fss/p1),s1z,'k'),
axis([0 2000 0 max(s1z)])
hold on
ylabel('Normalized Amplitude')
xlabel('Frequence (Hz)')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x=s2;
sp=fft(x);
nn=length(x);
sp(1:2)=zeros(size(1:2));
sp(nn-1:nn)=zeros(size(1:2));
x=real(ifft(sp));
% ************ FFT ****************
p1=length(x);
z=(2/p1)*abs(fft(x));
z=z(1:p1/2-1);
z(1,1)=0.5*z(1,1);
nff=length(z);
nf1=round(2*nff*upfreq/fss);
z=z(1:nf1);
subplot(3,2,4)
s2z=abs(z)/max(abs(z));
plot((0:nf1-1).*(fss/p1),s2z,'k'),
axis([0 2000 0 max(s2z)])
ylabel('Normalized Amplitude')
xlabel('Frequence (Hz)')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x=s3;
sp=fft(x);
nn=length(x);
sp(1:2)=zeros(size(1:2));
sp(nn-1:nn)=zeros(size(1:2));
x=real(ifft(sp));
% ************ FFT ****************
p1=length(x);
z=(2/p1)*abs(fft(x));
z=z(1:p1/2-1);
z(1,1)=0.5*z(1,1);
nff=length(z);
nf1=round(2*nff*upfreq/fss);
z=z(1:nf1);
subplot(3,2,6)
s3z=abs(z)/max(abs(z));
plot((0:nf1-1).*(fss/p1),s3z,'k'),
axis([0 2000 0 max(s3z)])
ylabel('Normalized Amplitude')
xlabel('Frequence (Hz)')
  1 Comment
Torsten
Torsten on 25 Feb 2023
1 Hz equals 2*pi rad/s
and
1 rad/s equals 1/(2*pi) Hz.
Can you take it from here ?

Sign in to comment.

Answers (1)

Bruno Luong
Bruno Luong on 25 Feb 2023
omega = 2*pi*f
with unit
omega [rad/s]
f [1/s]

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!