Wyznaczanie wartość średnia, skuteczna, międzyszczytowa.

\

Wrzucam TSS Czas jak sie otworzy w Matlabie to sa po kolei pliczki z komendami do roznych tam obliczeń czy cuś i nalezy powybierac poszczegolne i JAKOŚ polaczyc w calosc z wytycznymi zadania zeby pasilo tylko jak ? hehe a importuje sie to przez okienko co jest w matlabie u gory jak sie go odpali

function [xave,xrms,xpp]=cechy(x,T,d)
%disp('wartość średnia');
xave=sum(abs(x))/((T/d)+1);
%disp('wartość skuteczna');
xrms=sqrt(sum(x.^2)/((T/d)+1));
%disp('wartość międzyszczytowa');
xpp=max(x)-min(x);

cd ścieżka
load sygnal
[xave,xrms,xpp] = cechy(x, T, d)

Simulink

pkt 1
to się wpisuje w okienku command window (z prawej) w programie matlab
pkt 2
jak się ma otwarte okienko w którym się stawia bloczki to jest to opcja simulation -> simulaion parameters
pkt 3
zmienne stanu: q - położenie, q z kropką - prędkość, zmienna wyjścia- q- położenie, zmienna wejścia F(t)
macierz dx(t)/dt (2x1)=q z kropką; q z 2 kropkami. macierz A(2x2) = 0 1;-k/m -c/m. macierz x(t) (2x1) = q; q z kropką. macierz B (2x1) =0;1/m. macierz u(t) (1x1)=F(t)
q= macierz C(1x2) 1 0 . macierz x(t) (2x1) = q; q z kropką. macierz D(1x1) 0. macierz u(t) (1x1)=F(t)
pkt 4
impuls jednostkowy nie wiem, skok jednostkowy - wpisać te dane do bloku wejscie1 (step), przebieg harmoniczny - wpisać dane do bloku sine wave

Trajektoria fazowa to wykres x' od x i fi' od fi (prędkości od wychylenia) i robi się to za pomocą bloczku XY Graph dostępnego w sekcji sinks. Ma on dwa wejścia, do pierwszego doprowadzamy to co na osi x - czyli rozdzielamy "przewód" za drugim blokiem całkującym, a do drugiego doprowadzamy odgałęzienie przewodu za pierwszym blokiem całkującym.

load sygnal2;
x=sygnal2;
cumA=cumsum(x);
L=length(cumA(cumA<600000000));
nfft=2^nextpow2(L);
Y=fft(x,nfft)/L;
Fs=0.01;

f=Fs/2*linspace(0,1,nfft/2);
plot(f,2*abs(Y(1:nfft/2)));
title('widmo sygnalu');
xlabel('częstotliwość');
ylabel('|Y(f)|'),grid on;

T=1;

xave=sum(abs(x))/((T/Fs)+1);
xrms=sqrt(sum(x.^2)/((T/Fs)+1));
xpp=max(x)-min(x);
wk=xrms/xave;
wsz=max(x)/xrms;
wi=max(x)/xave;
hist(x); figure(gcf)

clear all
clc


load sygnal
d=0.01
T=1


disp('autokorelacja')
a=xcorr(y)
disp('histogram')
h=hist(y)
disp('wartość średnia')
xave=sum(abs(y))/((T/d)+1)
disp('wartość skuteczna')
xrms=sqrt(sum(y.^2)/((T/d)+1))
disp('wartość międzyszczytowa')
xpp=max(y)-min(y)

2.3

xsrekwa = sum(x.^2)/((T/d)+1); ---------wartość średniokwadratowa
xwspksz = xrms/xave; -----------współczynnik kształtu
xwarszc = max(x); -----------wartość szczytowa największa (warto zrobić jeszcze najmniejszą)
xwspszc = xwarszc/xrms -----------współczynnik szczytu
xwspimp = xwarszc/xave ----------- wpółczynnik impulsowości

---