[szczesliwa_patelnia]

#matlab #octave #studbaza #techbaza

Pomoże ktoś z generowaniem funkcji w Matlabie/Octave etc?

[shymon]

@szczesliwa_patelnia: Łapaj i się baw ;)

http://nopaste.info/50d2e90bc1.html

[szczesliwa_patelnia]

@shymon: Dzięki, analizuję powoli :)

[shymon]

@szczesliwa_patelnia: Powiem tak. W ramach ćwiczeń zrób sobie jeszcze sygnał trójkątny (czyli taki jak pilokształny, ale o dwóch zboczach - tu masz tylko rosnące, zrób jeszcze tak, że nie spada od razu do zera, tylko opada od wartości A do 0).

Jak zrobisz to, to dalej jeszcze sygnały bipolarne (czyli nie od 0 do A, lecz od -A do A).

[shymon]

@szczesliwa_patelnia: Nie ;) Bo tam obok jakby "doklejasz" do starego wektora ten stary wektor + ten nowy. Taka rekurencja trochę.

[szczesliwa_patelnia]

@shymon: Tak, doczytałem teraz wytłumacznie.

[szczesliwa_patelnia]

@shymon: Na pewno kod jest OK? Coś prz tym przykładzie ze mną nie współpracuje. Nie chce wyświetlać, próbowałem kilka razy.

http://pbrd.co/HoUkLk

[szczesliwa_patelnia]

@shymon: Wyżej pominąłem N. Po jego dodaniu to nic nie zmienia, program nadal czka na komendy gdy powinien wyświetlać wykres.

[shymon]

@szczesliwa_patelnia: Kurcze, szczerze to nie wiem. W matlabie ładnie chodzi, z Octave nigdy nie korzystałem. Może ten średnik po wektorze t usunąć trzeba?:/ Bo z tego co widzę, to w innych liniach pousuwałeś i przeszło.

[szczesliwa_patelnia]

Spróbuję bez średnika, ale wcześniej nie miał wpływu na takie rzeczy. Wygląda, jakby pętla się nie zamknęła.

[shymon]

@szczesliwa_patelnia: A size(x) zamiast length(x) ?

[szczesliwa_patelnia]

@shymon: Sprawdzę, ale w pętli nic nie zwiększa i a pętla na tym polega. Tu jest błąd na 99%.

[szczesliwa_patelnia]

@shymon: No te zmiany nic nie dają.

[shymon]

@szczesliwa_patelnia: No to niestety, ale nie wiem...

W matlabie działa, więc sama matematyka ok, lecz sposób implementacji pewnie trochę inny będzie. Niestety musisz sam powalczyć :/

[szczesliwa_patelnia]

@shymon: Spróbuję, ale sam działania do końca nie rozumiem + jestem głupi w programowanie. Nie mniej, na logikę - pętla opiera się na zwiększaniu i a te nigdzie później nie występuje, jak to ma działać?

[shymon]

@szczesliwa_patelnia: Ah. To może w ostatniej linii pętli dodaj i=i+1? Bo definiując w matlabie for i=1:4 oznacza to, że pętla wywoła się 4 razy, a 'i' zwiększa się samo z każdym obrotem pętli.

Ta pętla ma na celu po prostu doklejenie kolejnego zbocza do całego wektora x (kolokwialnie w tym przypadku zbocze to x1).

Jeśli wywalisz pętlę i w następnej linii po x1 wstawisz po prostu

x=[x1 x1 x1 x1], to

[szczesliwa_patelnia]

@shymon: Zobaczę za chwilę, tu jest inaczej, krócej opisane:

http://www.teleinformatyk.pl.hostingasp.pl/generacja-sygnałów-okresowych.aspx

[shymon]

@szczesliwa_patelnia: Ano, bo tu jest wykorzystanie gotowej funkcji sawtooth, która mi w ogóle na myśl nie przyszła :)

Zrobić to można na miliard sposobów :P

[szczesliwa_patelnia]

@szczesliwa_patelnia: Ano, bo tu jest wykorzystanie gotowej funkcji sawtooth, która mi w ogóle na myśl nie przyszła :)

@shymon: Ok, w zasadzie skończyłem. Potrzebuję tylko jeszcze zsumować sygnały, podobno jest funkcja subplot.

Liniowy

t=[0:0.00125:1];

y=sawtooth(2pi4*t);

plot(t,y)

Sinus

tp2=0.00125;

f=16;

t2=[0:tp2:1-tp2];

y2=sin(2pif*t2);

plot(t2,y2)

[shymon]

@szczesliwa_patelnia: Co teraz rozumiesz poprzez "zsumować" sygnały?

Wykreślenie ich na jednym wykresie, a zsumowanie sygnałów to inne sprawy.

Sumowanie, to byłoby coś pokroju y3 = y+y2 (pod warunkiem, że mają tyle samo próbek) i wyplotowanie potem (t, y3).

Z kolei subplot tworzy w jednym oknie jakby dwa oddzielne wykresy.

Jeśli chcesz wyplotować je na jednym, to dodaj po prostu komendę "hold on" pomiędzy plotami.

Będzie wyglądało tak:

plot(t,y)

hold on

(blablalba