Zaloguj się
Blog na Matlablog
Forum polskich użytkowników
 
UŻYTKOWNICY GRUPY PROFIL Zaloguj się, by sprawdzić wiadomości FAQ
 



Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Badanie sygnału okresowego Idź do strony 1 2  Następny
Forum MATLAB Strona Główna-> MATLAB
Post Wysłany: 16 Marca 2018, Pią 12:55 pm Temat postu: Badanie sygnału okresowego Odpowiedz z cytatem
 
AUTOR:
lucky_
Początki pisania


Dołączył: 16 Mar 2018
Posty: 8


Ogląda profil użytkownika Wyślij prywatną wiadomość
Witam,

dopiero zaczynam przygodę z MatLabem i nie wiem jak poradzić sobie z pewnym problem. Dokładnie chodzi o to, że mam sygnał U1t i muszę do niego wykonać symulację w postaci funkcji czasu i obliczyć wartość średnią, tylko nie wiem za bardzo jak go zapisać w MatLabie. Sygnał jest rozbity na dwa wzory i każdy ma inny zakres i nie wiem jak to zapisać, aby MatLab widział go jako jeden wzór. Sygnał dodany w załączniku.


 

Post Wysłany: 17 Marca 2018, Sob 12:20 am Temat postu: Co zrobić z tą formułą Odpowiedz z cytatem
 
AUTOR:
Jacek Konieczny
Może pisać książki


Dołączył: 31 Mar 2009
Posty: 349
Skąd: Poznań


Ogląda profil użytkownika Wyślij prywatną wiadomość Odwiedź stronę autora
Można to zapisać jako jeden wzór wykorzystując funkcję 'sign' :

2 * U * t * ( T/2 - sign (t) * t ) / T^2

Funkcja `sign` ( funkcja 'signum' w języku matematyki ) :

https://www.mathworks.com/help/matlab/ref/sign.html

Do uśredniania można wykorzystać funkcję 'mean' :

https://www.mathworks.com/help/matlab/ref/mean.html


 

Post Wysłany: 17 Marca 2018, Sob 3:27 pm Temat postu: Re: Co zrobić z tą formułą Odpowiedz z cytatem
 
AUTOR:
lucky_
Początki pisania


Dołączył: 16 Mar 2018
Posty: 8


Ogląda profil użytkownika Wyślij prywatną wiadomość
Jacek Konieczny napisał:
Można to zapisać jako jeden wzór wykorzystując funkcję 'sign' :

2 * U * t * ( T/2 - sign (t) * t ) / T^2

Funkcja `sign` ( funkcja 'signum' w języku matematyki ) :

https://www.mathworks.com/help/matlab/ref/sign.html

Do uśredniania można wykorzystać funkcję 'mean' :

https://www.mathworks.com/help/matlab/ref/mean.html


Czyli wzór będzie wyglądał tak jak napisałeś - 2 * U * t * ( T/2 - sign (t) * t ) / T^2; dla całości?
Potrzebowałbym jeszcze pomocy z zapisaniem tego zakresu dla tego wzoru.


 

Post Wysłany: 17 Marca 2018, Sob 5:56 pm Temat postu: Zakres zmienności dla wzoru Odpowiedz z cytatem
 
AUTOR:
Jacek Konieczny
Może pisać książki


Dołączył: 31 Mar 2009
Posty: 349
Skąd: Poznań


Ogląda profil użytkownika Wyślij prywatną wiadomość Odwiedź stronę autora
Zakres jest zapisany w Twoim skrypcie ( w deklaracji zmiennej `t` ) :

t=(-T/2T/N(5*T/2)-(T/N)); %okres

W oryginalnym zapisie były jakieś emotikony . Moim zdaniem powinien wyglądać raczej tak :

t = [ ( -T/2 ) : dt : ( T/2 ) ]

gdzie "krok" `dt` sam sobie ustalasz ( np. na podstawie N ).


 

Post Wysłany: 4 Kwietnia 2018, Sro 10:40 am Temat postu: Odpowiedz z cytatem
 
AUTOR:
drewbud
Pisze


Dołączył: 16 Paź 2017
Posty: 20


Ogląda profil użytkownika Wyślij prywatną wiadomość Odwiedź stronę autora
czy to nie wpłynie błędnie na zmienną?

https://www.knall.com.pl/c174,Przyciemnienie-i-dekoracja


 

Post Wysłany: 17 Kwietnia 2018, Wto 8:37 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
lucky_
Początki pisania


Dołączył: 16 Mar 2018
Posty: 8


Ogląda profil użytkownika Wyślij prywatną wiadomość
Nie chcę zakładać nowego tematu, więc napiszę tu. Jeśli chodzi o U1 to już zrobiłem i wyszedł wykres jak w załączniku, natomiast teraz mam problem z U2. Nie jestem pewien, ale wykresy powinny wyjść różne dla różnych składowych wyjść, a u mnie wychodzą takie same jak U1. Wklejam dwie składowe, może ktoś wie gdzie tu jest błąd?

Kod:
U2t=((4*U/pi^3)*(sin(w*t)));
        
figure(2);
        
subplot(2,1,1);
        
plot(t,U2t);
        
title('Symulacja dla 1. składowej.');
        
xlabel('czas t [s]');
        
ylabel('amplitua U [V]');
        
grid on;
        
        
y=fft(U2t);
        
m=abs(y);
        
subplot(2,1,2);
        
stem(fk(1:N/2),m(1:N/2),'.');
        
title('Widmo częstotliwościowe');
        
xlabel('Częstotliwość [Hz]');
        
ylabel('amplituda U [V]');
        
grid on;
 
          
U2t=((4*U/pi^3)*(sin(w*t)))+(((4*U/pi^3)*(sin(3*w*t))/(3^3)));
        
figure(3);
        
subplot(2,1,1);
        
plot(t,U2t);
        
title('Symulacja dla 2. składowej.');
        
xlabel('czas t [s]');
        
ylabel('amplitua U [V]');
        
grid on;
        
        
y=fft(U2t);
        
m=abs(y);
        
subplot(2,1,2);
        
stem(fk(1:N/2),m(1:N/2),'.');
        
title('Widmo częstotliwościowe');
        
xlabel('Częstotliwość [Hz]');
        
ylabel('amplituda U [V]');
        
grid on;


 

Post Wysłany: 18 Kwietnia 2018, Sro 6:57 pm Temat postu: Przebieg U2 Odpowiedz z cytatem
 
AUTOR:
Jacek Konieczny
Może pisać książki


Dołączył: 31 Mar 2009
Posty: 349
Skąd: Poznań


Ogląda profil użytkownika Wyślij prywatną wiadomość Odwiedź stronę autora
Proszę zwrócić uwagę na sposób, w jaki sposób zdefiniowany jest przebieg U2 w załączonym skrypcie ; zdefiniowany jest przebieg ( sinusoidalny ! ) o częstotliwości "w" ( podczas gdy powinien zdefiniowany przebieg o składowych nieparzystych ).


 

Post Wysłany: 20 Kwietnia 2018, Pią 4:57 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
lucky_
Początki pisania


Dołączył: 16 Mar 2018
Posty: 8


Ogląda profil użytkownika Wyślij prywatną wiadomość
Nie wiem czy dobrze zrozumiałem, ale chodzi o to "w" w tym wzorze? To jest jako omega.


 

Post Wysłany: 20 Kwietnia 2018, Pią 6:47 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
Jacek Konieczny
Może pisać książki


Dołączył: 31 Mar 2009
Posty: 349
Skąd: Poznań


Ogląda profil użytkownika Wyślij prywatną wiadomość Odwiedź stronę autora
Właśnie ! We wzorze "teoretycznym" były obecne składowe o częstotliwościach : 3 omega oraz 5 omega.


 

Post Wysłany: 20 Kwietnia 2018, Pią 8:49 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
lucky_
Początki pisania


Dołączył: 16 Mar 2018
Posty: 8


Ogląda profil użytkownika Wyślij prywatną wiadomość
O to chodzi?

Kod:
U2t=(4*U*(sin(w*t))/(pi^3))+(4*U*(sin(3*w*t))/(pi^3)*(3^3));
        
figure(2);
        
subplot(2,1,1);
        
plot(t,U2t);
        
title('Symulacja dla 1. składowej.');
        
xlabel('czas t [s]');
        
ylabel('amplitua U [V]');
        
grid on;
        
        
y=fft(U2t);
        
m=abs(y);
        
subplot(2,1,2);
        
stem(fk(1:N/2),m(1:N/2),'.');
        
title('Widmo częstotliwościowe');
        
xlabel('Częstotliwość [Hz]');
        
ylabel('amplituda U [V]');
        
grid on;
 
        
U2t=(4*U*(sin(w*t))/(pi^3))+(4*U*(sin(3*w*t))/(pi^3)*(3^3))+(4*U*(sin(5*w*t))/(pi^3)*(5^3));
        
figure(3);
        
subplot(2,1,1);
        
plot(t,U2t);
        
title('Symulacja dla 2. składowej.');
        
xlabel('czas t [s]');
        
ylabel('amplitua U [V]');
        
grid on;
        
        
y=fft(U2t);
        
m=abs(y);
        
subplot(2,1,2);
        
stem(fk(1:N/2),m(1:N/2),'.');
        
title('Widmo częstotliwościowe');
        
xlabel('Częstotliwość [Hz]');
        
ylabel('amplituda U [V]');
        
grid on;


 

Post Wysłany: 20 Kwietnia 2018, Pią 10:08 pm Temat postu: Ostatnia wersja skryptu. Odpowiedz z cytatem
 
AUTOR:
Jacek Konieczny
Może pisać książki


Dołączył: 31 Mar 2009
Posty: 349
Skąd: Poznań


Ogląda profil użytkownika Wyślij prywatną wiadomość Odwiedź stronę autora
Właśnie o to chodziło. W poprzedniej wersji skryptu była generowana tylko jedna składowa ( o częstotliwości "w" ).


 

Post Wysłany: 21 Kwietnia 2018, Sob 9:00 am Temat postu: Odpowiedz z cytatem
 
AUTOR:
lucky_
Początki pisania


Dołączył: 16 Mar 2018
Posty: 8


Ogląda profil użytkownika Wyślij prywatną wiadomość
Dodałem jeszcze 3 składową o wartości 7 omega i wykresy prezentują się tak: http://wstaw.org/m/2018/04/21/screenshot.179.jpg


 

Post Wysłany: 21 Kwietnia 2018, Sob 10:32 am Temat postu: WyĹźsze składowe Odpowiedz z cytatem
 
AUTOR:
Jacek Konieczny
Może pisać książki


Dołączył: 31 Mar 2009
Posty: 349
Skąd: Poznań


Ogląda profil użytkownika Wyślij prywatną wiadomość Odwiedź stronę autora
Trochę się to nie zgadza dla 2 i 3 skladowej ; widmo pokazuje te same składowe. Może dałoby się to sprawdzić w ten sposób, że wygenerowany przebeg ( np. trójton ) zapisuje się do pliku WAV ( przy pomocy funkcji "wavwrite" ), a następnie sprawdza się widmo tego wygenerowanego i zapisanego przebiegu przy pomocy jakiegoś programu "zewnętrznego" ( np. przy pomocy darmowego edytora dźwięków o nazwie `Audacity` ).


 

Post Wysłany: 21 Kwietnia 2018, Sob 12:23 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
lucky_
Początki pisania


Dołączył: 16 Mar 2018
Posty: 8


Ogląda profil użytkownika Wyślij prywatną wiadomość
Sprawdziłem jeszcze raz kod i miałem błąd, po poprawie wygląda to tak: http://wstaw.org/m/2018/04/21/screenshot.180.jpg


 

Post Wysłany: 21 Kwietnia 2018, Sob 1:55 pm Temat postu: Proporcke amplitud składowych Odpowiedz z cytatem
 
AUTOR:
Jacek Konieczny
Może pisać książki


Dołączył: 31 Mar 2009
Posty: 349
Skąd: Poznań


Ogląda profil użytkownika Wyślij prywatną wiadomość Odwiedź stronę autora
Ciekawe, że im wyższa składowa, tym wyższa amplituda tej składowej. Natomiast we wzorze "definicyjnym" ( we wzorze "teoretycznym" ) amplituda kolejnej, wyższej składowej maleje jak mnożnik częstotliwości podniesiony do potęgi trzeciej.


 

Forum MATLAB Strona Główna-> MATLAB
Wyświetl posty z ostatnich:   
Idź do strony 1 2  Następny

Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Wszystkie czasy w strefie CET (Europa)

Skocz do:  

Statystyki forum:



Od dnia 08.06.2006 forum odwiedzano 45741380
Najwięcej użytkowników 266 było obecnych 19 Lutego 2015, Czw 7:03 pm

Aktualnie online:




Najnowsze posty na forum:
Telewizja online za darmo  (24 Maj 2018, Czw 11:00 pm)
Ja, Godard 2017 Obejrzyj Cały Film Online Lektor PL  (23 Maj 2018, Sro 7:44 pm)
Wykres i widmo zerowej składowej  (23 Maj 2018, Sro 3:09 pm)
(s01e10) Nawiedzona panna młoda 1 Sezon 10 Odcinek Online  (22 Maj 2018, Wto 10:17 pm)
Nawiedzona panna młoda s01e09 transmisja na Ĺźywo  (22 Maj 2018, Wto 10:17 pm)
Nawiedzona panna młoda s01e09 gdzie obejrzeć online?  (22 Maj 2018, Wto 10:16 pm)
(s01e09) Nawiedzona panna młoda 1 Sezon 9 Odcinek Online  (22 Maj 2018, Wto 10:16 pm)
Nawiedzona panna młoda s01e08 gdzie obejrzeć online?  (22 Maj 2018, Wto 10:16 pm)
(s01e08) Nawiedzona panna młoda 1 Sezon 8 Odcinek Online  (22 Maj 2018, Wto 10:15 pm)
(s01e07) Nawiedzona panna młoda 1 Sezon 7 Odcinek Online  (22 Maj 2018, Wto 10:15 pm)
Twoje prawa:
Nie możesz pisać nowych tematów
Nie możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach
Nie możesz załączać plików na tym forum
Nie możesz ściągać plików na tym forum