Jak to zrobić żeby dany sygnał miał taką własność :
- o dominującej gęstości widmowej w zakresie od 0 do 10 Hz.
%Dane:
Fs= 1; % czętotliwośc próbkowania
t=0:1/Fs:100; % czas
f= 100 ; % czętotliwość Hz
%wartość średnia równa 0, odchylenie standardowe równe 10kN
syg=sin(2*pi*f*t);
syg1 =syg+randn(1, length(syg));
syg1=syg1-mean(syg1);
syg1=10*syg1/std(syg1);
F=syg1;
plot(t, F)
title('Przebieg siły losowej F(t)');
xlabel('t[s]');
ylabel('F[kN]');
hold on;
grid on ;
% odczyt średniej i odchylenia
stats = [mean(syg1) std(syg1)]
%gęstość widowa zakres 0-10 Hz
- o dominującej gęstości widmowej w zakresie od 0 do 10 Hz.
%Dane:
Fs= 1; % czętotliwośc próbkowania
t=0:1/Fs:100; % czas
f= 100 ; % czętotliwość Hz
%wartość średnia równa 0, odchylenie standardowe równe 10kN
syg=sin(2*pi*f*t);
syg1 =syg+randn(1, length(syg));
syg1=syg1-mean(syg1);
syg1=10*syg1/std(syg1);
F=syg1;
plot(t, F)
title('Przebieg siły losowej F(t)');
xlabel('t[s]');
ylabel('F[kN]');
hold on;
grid on ;
% odczyt średniej i odchylenia
stats = [mean(syg1) std(syg1)]
%gęstość widowa zakres 0-10 Hz