Na wstępie chciałem się przywitać :) Wykonuję projekt, w którym jednym z celów jest pomiar kąta przesunięcia fazowego (potrzebuje przetestować kilka metod). Matlabem zajmuję się dopiero od niedawna. Znalazłem książkę z 2000 roku: J. Gajda R. Sroka, Pomiary kąta fazowego. Na końcu książki umieszczone są przykładowe programy napisane w MATLABie, które...
Widzę, że nie bardzo rozumiesz, jak działa Matlab. Matlab to operacje na wektorach. W pętli, z każdą iteracją powinieneś nadawać wartości kolejnym elementom wektora, Ty tymczasem za każdym razem zmieniasz jego wartość. Teraz, jeżeli na jednej osi chcesz n, na drugiej k, a na trzeciej wartości, to potrzeba wykresu trójwymiarowego. Poza tym, z tego, co...
Witam, mam problem w poprawnym zapisaniu transmitacji ciągłej. Przy zapisie dyskretnym nie mam problemu. Moja transmitacja to G(s)=5s^2/(s^3+2s+1) Używam funkcji zpk(Z,P,k), przy czym Z oraz P mają być w postaci wektorów. Wektory staram się wyznaczyć w sposób [syntax=matlab]polynomial = [1 2 1][/syntax] następnie [syntax=matlab]poly_roots = roots(polynomial)...
Potrzebuję pomocy z prostym (prosty to pojęcie względne :P ) zadaniem w matlabie . Jestem zielona jesli chodzi o ten program a potrzebuje tego na studia. Bardzo proszę o pomoc i z góry dziękuję! ;) Wyznacz wszystkie (o ile istnieją) kombinacje liniowe wektorów ai dające wektor b a) a1= (-1,1,5) a2=(-2,1,3) a3=(1,0,2), b=(-3,1,4) b) a1= (2,-1,4,3) a2=(0231)...
Trochę inne są te kody. W C masz funkcję, która da Ci wynik tylko dla jednej liczby, podanej jako argument. Twój kod w Matlabie wykona się dla każdej liczby z wektora x, a ostatecznie arch będzie równy wartości logarytmu dla ostatniej liczby z wektora x.
Zrób wektor z wartościami czasu np. co 10 ms a następnie policz wartość funkcji dla każdej wartości czasu w wektorze.
witam! Muszę wykonać układ sterowania wektorowego w matlabie silnika indukcyjnego. Nie wiem jak wykonać taki układ. Bardzo proszę o porady!
Najwygodniej jest skorzystać z symulacji Simulinkowej, która rozwiąże każdą z podanych przez Ciebie kwestii. Taką symulację mógłbyś wywołać z m-pliku za pomocą komendy sim , w modelu możesz też wyprowadzić przebiegi do pamięci za pomocą bloku 'To Workspace' i później wykonać wykresy z poziomu skryptu matlaba. Jeżeli musisz koniecznie wykonać wszystko...
Witam. Rozwiązałem sobie układ równań nieliniowych metodą Newtona-Raphsona w Matlabie. Oto równanie: A oto rozwiązanie: plik fun.m: [syntax=c] function z=fun(x) z=zeros(size(x)); z(1) = x(1)+2*x(2)-2; z(2) = 2*x(1)^2+4*x(2)^2-4; [/syntax] plik JAC.m: [syntax=c] function j = JAC(x) j=[1 2; 4*x(1) 8*x(2)]; [/syntax] [syntax=c] function[xsol, k] = NewRhp(f,j,x0,eps,maxiter)...
Witam dostałem do analizy zadanie ze schematu. Obwód należało rozwiązać metodą oczkową. Wszystko zastało zaimplementowane bilans prądów się zgadza ale nie zgadza się za cholere bilans oczka 2 i bilans mocy, czy ktoś mógłby coś podpowiedziec? Dane do zadania jakie przyjęto: Ra=5 om Rb=3 om Rc=6 om L1= 0.5 H L2=0.1 H C=0.02 F k=0.8 wpółczynnik sprzezenia...
Witam, mam problem z rysowaniem wykresu FFT w Matlabie, otóż mam funkcje: x=(2-400*t).*cos(2000*pi*t) Dla tej funkcji dobrałem częstotliwość próbkowanie f=16kHz, i przyjąłem N=16 próbek na okres, czyli analizując 2 i 4 okresy kolejno otrzymam 32 i 64 próbki. Kod mojego programu w matlabie jest następujący: >> t=linspace(0,2*0.001,10000); >> A=(2-400*t)...
Witam. Proszę o pomoc w napisaniu a dokładnie zmodyfikowaniu następującego kodu: clear all; close all; fs=8000; fc=100; fd=25; N=8; x=randint(300,1,N); y=dmod(x,fc,fd,fs,'psk',N); SZUM=[-100:10]/10; for numer_SZUM=1:length(SZUM); A=SZUM(numer_SZUM); yszum=awgn(y,A); yszum_zdem=ddemod(yszum,fc,fd,fs,'ps...
Genario , jeżeli chodzi o matlaba, to nie za wiele pomogę, bo nie używam. :( A z tym skalowaniem to chodzi o coś takiego, że nie znasz częstotliwości próbkowania, ale mimo wszystko możesz przedstawić na skali częstotliwości poszczególnych prążków, właśnie w postaci znormalizowanej. mając widmo sygnału (rzeczywistego) o N próbkach składające się z N/2+1...
Witajcie, Od niedawna mam laborki z algorytmów obliczeniowych, a na nich matlab. Mieliśmy za zadanie zrobić skrypt pracujący na wprowadzonych przez użytkownika wektorach. Z początku wszystko fajnie, ale przy takim zapisie: n = input('Wprowadz liczbe elementow wektora: '); Wystarczyło podać liczbę niedodatnią lub niecałkowitą(choć w zasadzie...
Witam Otóż ostatnio zabrałem się za pisanie symulatora n ciał w matlabie. Bazując na wzorze z Rozpisałem sobie wektory kombinacji dla masy i położenia tj. w każdym wierszu posiada zapisane wektory cj (pozycje wszystkich ciał poza tym aktualnie liczonym) dla danego ciała. Wszystko niby ładnie wygląda i dobrze się liczy, ale w przy puszczeniu symulacji,...
Mam za zadanie zrobić filtr IIR dla zadanych wartości (w matlabie). Dane: f=1000 fs= Za mało informacji. Co to jest f=1000? Prawdopodobnie częstotliwość próbkowania, którą przyjęło się oznaczać fs (od sampling frequency). Co to zatem jest fs? Może wektory sygnałów testowych o częstotliwościach 100, 300 i 400, którymi masz sprawdzić filtr? Masz zrobić...
Witam serdecznie. Mam takie pytanie, dotyczące transmitancji operatorowej w Matlabie. Powiedzmy, że mam taką transmitancję: L=1 M=[1 2 3] Gdzie L to rzecz jasna licznik a M to mianownik. No i tak, potrafię to zasymulować w simulinku. Mogę z tego stworzyć np. transfer function powidzmy tak: sys=tf(L,M); No i tu się pojawia mój problem, jak obliczyć z...
Witam MATLAB potężnym i wszechstronnym narządem jest i co do tego nikt nie ma wątpliwości. Matlab udostępnia co najmniej kilka (parfor, batch, job/task, itp.) sposobów na zrównoleglenie wykonywanych zadań. Mam do napisania pewien program w jakiś sposób symulujący przetwarzanie danych w czasie rzeczywistym i mam problem z wyborem odpowiedniego sposobu...
Ponieważ Matlab traktuje Twoje Kp jako wektor. Jeżeli Kp=[1.8,2.5,3.9,4.0] to można się "dostać" do elementu np. tak a=Kp(1) -> a=1.8; W Twoim przypadku indeks tablicy nie będzie wartością całkowita i wywal błąd. To wyjaśnienie. Jak po Kp namalujesz "*" to regulator będzie się liczył :-)
Witam, "" Jak k jest wektorem oraz Tmax to może te wektory wprowadzić do "MATLAB Function" i zastosować porównania w celu przełączenia sygnału. Ja bym tak na początek podszedł do tego problemu. Pozdrawiam,
x(j:k)- elementy wektora wierszowego x o numerach od j do k • A(i,:)- wszystkie elementy w wierszu i macierzy A • A(i,j:l)- wszystkie elementy w wierszu i macierzy A o numerach od j do l • A(i:k,j:l)- wszystkie elementy w kolumnach od j do l wierszy od i do l • A(x,j:l)- wszystkie elementy w kolumnach od j do l w wierszach macierzy...
Pomimo długoletniego stażu w elektronice w tej dziedzinie niestety brak doświadczenia ....Potrzebuję pomocy ,wskazówek przy napisaniu programu przedstawiającego użycie matematyki w robotyce za pomocą MATLABA .....Ma on przedstawić np: _ wyświetlić wynik wyznacznika macierzy wejściowej (3*3 , 4*4) _ dla określonych 4 macierzy (wybór 3x3 lub 4x4) wyświetlić...
Witam, Mam zadanie zaimplementowania w Matlabie filtru Kalmana dla modelu wahadła odwróconego. Dostałem gotowy m-plik z Matlaba, który wyznacza filtr Kalmana dla przykładowego modelu(który został wzięty w komentarze). Zapisałem mój model(wahadło odwrócone) w przestrzeni stanu i wrzuciłem go do tego m-plika i zaczął się problem z którym nie mogę sobie...
Ja tam nie wiem, jak to w Matlabie wygląda, ale w MathCADzie by to jakoś tak wyglądało (zakładamy, że wszystkie indeksy tablic zaczynają się od zera) 0. Próbka dźwięku to wektor (V) 1. wybieramy ilość próbek do przetworzenia przez jedną transformatę (N) 2. Tworzymy macierz M o odpowiedniej wielkości: 2.1 liczba wierszy - dla transformaty liczonej z...
W C, bo później muszę to ruszyć na procesorze. Z wavem w Matlabie nie ma problemu. Tylko, że u mnie szału nie ma z C, a w takim wavie to trzeba by jakoś pomijać nagłówek... Nie chciałem komplikować sprawy, a że ten programik co podałem linka chce same próbki sygnału to robiłem to zapisując je z Matlaba tak (btw. nie ma tu na forum opcji listingu Matlabowego...
Ten wzorek można przerobić. Trzeba zauważyć, że cosh(x) = (exp(x)+exp(-x)) / 2 sinh(x) = (exp(x) - exp(-x)) / 2 i prowadź następujące oznaczenia : k = SQRT(15990142) / 1009 w = 150 / 377 a = 4139 /1009 d = 4139*SQRT(15990142) / (2*7995071) Jak to przekształcisz wszystko otrzymasz (kolega Quartz to zauważy również), że wzorek będzie : f(t) = w*(2 - (d+1)*exp((k-a)*t)...
Z chęcią pomogę, ale na gotowca nie licz. ad 1. po kolei: -wczytać parametry a-f: tu się pobaw -wygenerować wektor x z wartościami w przedziale z pewnym krokiem: x = xmin : krok : xmax; -w pętli for dla każdego x obliczyć wartość y -machnąć jakiś wektor z zerami (o długości wektora x) - narysować y w zależności od x (funkcja plot) - policzyć pierwiastki(funkcja...
Całość udało się zrobić za pomocą Toolboxa "System Identification Toolbox", jedyny haczyk był w dobraniu odstępów między chwilami próbkowania - domyślnie był "1" i program pobierał jedynie 7 próbek z wektora. Projekt z podstaw automatyki skończony, a w/w doktor był wyraźnie z niego zadowolony...stąd 5.5 :)
Witam chciałem zrobic prosty filtr dolnoprzepustowy za pomocą splotu (zgodnie z twierdzeniem ze splot w dziedzinie czasu odpowiada mnożeniu widm w dziedzinie czestotliwosci). Chciałem sprawdzić to w matlabie- odfiltrowac z sumy dwóch sinusów skladową o wyzszej czestotlowsci za pomoca funkcji spotu conv : %Implementacja filtru FIR clear all; close all;...
Zajrzalem tutaj wczoraj, i tak samo cos mi sie wydawalo, ze te rownanai do konca liniowe nie sa. Rownania liniowe to rownania postaci: ax+by=c dx+ey=f Tutaj natomiast mamyiloczyny zmiennych. Proponuje po prostu przeszukanie przestrzeni w poszukiwaniu rozwiazan. Czlon po lewej stronie jest zawsze taki sam, wiec dla kazdego pktu w pewnej czesci przestrzni...
Dziękuję za odpowiedź. Program działa prawie tak jak chciałem jedyne co jeszcze oczekiwałem to wektor który będzie miał 10 wyników lecz to już zrobiłem sam. Zastanawiam się jeszcze jednak co w przypadku gdy ten drugi wektor byłby dłuższy. Pozdrawiam Dodano po 4 Popełniłem głupi błąd przy wypisywaniu zmiennych. W przypadku gdy drugi wektor byłby większy...
Nie wiem czy ci to cos pomoze ale np jak masz dwa wektory tej samej dlugosci i chcesz pomnozyc wszystkie wyrazy miedzy soba (pierwszy razy pierwszy , drugi razy drugi itd) to zamiast for i=1:n C(i)=A(i)*B(i) end uzywasz C=A.*B Tak samo mozna mnozyc odpowiednie kolumny lub wiersze macierzy
Witam mam prośbę. Zrobi ktoś takie zadańka, są krótkie. POZDRAWIAM!! MATLAB 1. Napisz m-plik, wyniku którego powstanie jeden rysunek z dwoma wykresami funkcji y=asin(wt+fi) w przedziale (-pi, +pi) z opisem osi i tekstem wskazującym na wersje danych a)a=5, w=7, fi=45 B) a=7, w=5, fi=-45 2. Wygeneruj wektor liczb x z przedziału (a,b) używając polecenia...
Super dzieki bardzo za pomoc! Jesli ktos moze dorzucic swoje 3 grosze to prosze bardzo.. I ten efekt Gibbsa... Pozdrawiam Dodano po 4 aproksymuje coś wychodzi ale raczej nie to co potrzeba, jestem raczkujący i nawet nie wiem czy to co pisze(polecenia) matlab "rozumie" Ponizej przedstawiam to co wpisuje: Moze napisze po kolei co robie - bedzie łatwiej...
Ja mam pytanie związane z tematem FFT bo nie udało mi się znaleźć nigdzie odpowiedzi na nurtujący mnie problem. W matlabie troszkę pracowałem i akurat z FFT walczyłem. Gdy przerobiłem sygnał podany przez agent_00 swoim algorytmem nie do końca uzyskałem to co kowales . Problem mój polega na tym że w matlabie jako jeden z parametrów polecenia FFT jest...
Witam, chcę odfiltrować z sygnału pewną częstotliwość przy pomocy fft. Wykonuję filtrowanie w matlabie. Problem pojawia się przy sprawdzeniu. Jeśli analizuję cały filtrowany sygnał wszystko jest w porządku ale jeśli poddam analizie fragment sygnału odfiltrowana częstotliwość pojawia się ponownie. Myślałem że jest to kwestia nieprecycyjnego wyznaczenia...
To co udało mi się do tej pory "osiągnąć" zamieszcze w pliki tekstowym tak aby można było wkleić do Matlaba i zobaczyć wynik. Miszmasz o którym pisałem to to iż co prawda wychodzą jakieś piki lecz niestety zupełnie nie przypominają wyższych harmonicznych tzn nie są wielokrotnościami 50Hz... Nie wykorzystuje czytania z pliku excela tylko kopiuje dane...
Zapewne studenci elektrotechniki lub elektroniki mają rachunek wektorowy opanowany więc może pomożecie ;) Sprawa jest taka jak na obrazku. - mamy wektor "W" który rozkładamy na składowe Wx i Wy - obracamy układem współrzędnych o kąt alfa ( UWAGA, wektor zostaje w miejscu ) - teraz ten wektor rozłoży się już na inne składowe Wx' i Wy' bo zrzutowane na...
Witam, na starcie chciałem przeprosić jeśli umieszczam temat w nieodpowiednim dziale, ale nie znalazłem nic co by odpowiadało mojemu problemowi. Mam do wykonania projekt wyznaczenia modelu kinematycznego trójkołowego robota mobilnego z kołami sferycznymi, a następnie zbudowanie do niego układu sterowania, tak aby jechał po zadanej trajektorii. Model...
Witam Tworzysz funkcję o kształcie "przebiegu prostokątnego" :) o wartościach 0 i 1 w odpowiednich przedziałach. Jeżeli znasz postać analityczną funkcji opisującej "twój sinus" to można wyznaczyć to analitycznie, jeżeli nie, to można odczytać z wykresu. Następnie mnożysz szum przez tą funkcję. Zastosuj operator mnożenia "punktowego" i pamiętaj o równej...
Odświeżam temat! Dotarłem do końca swojego projektu ale mam wrażenie, że coś się tutaj jednak nie zgadza. Temat wciąż ten sam, tylko schemat wykorzystałem nieco trudniejszy... Zamieszczam schemat na podstawie którego udowadniałem swoją tezę: ...proszę o bardziej zaawansowanych urzytkowników o prześledzenie mojego listningu programu. % DANE: w=400...
Nie wiem jak w Matlabie, ale w SciLabie jest tak, że .wav wczytany zostaje jako wektor liczb zmiennoprzecinkowych z zakresu od -1 do +1. W związku z powyższym w celu osiągnięcia żądanej ilości poziomów kwantyzacji można zrobić tak: x=wavread('jakis_plik.wav'); x = round(x*2^(ilość_bitow - 1)) / 2^(ilość_bitow - 1); Po...
Witam. Robię projekt w oparciu o Matlaba. potrzebuję zapisać równanie różniczkowe tak, aby można było wykorzystać funkcje ode... do rozwiązania.. sprowadziłam równanie do takiej postaci: xq2=(-WT*D*W.*(diff(xq,2,t)) + WT*U*W.*(diff(xq,t)) + WT*F)\(WT*E*W); jednak dostaję bład w formie ??? Error using ==> sym.diff at 43 Do not recognize argument number...
każdy wielomian w matlabie mozesz przedstawić jako wektor wspolczynników, w Twoim przypadku polynomial = [0.03 0.4 1] a pozniej użyć komendy poly_roots = roots(polynomial) albo odrazu poly_roots = roots([0.03 0.4 1])
Witam, mam problem z wczytaniem pliku audio i zmianie wartości w wektor binarny. Mój kod to funkcja o następującej formie: [syntax=matlab] function [ y, fp ] = audiot( x ) [sound fp]=audioread(x); dod=abs(sound) dzwiek=(uint16(dec2bin(dod))) [/syntax] Problem polega na tym, że wynik takiego zabiegu to macierz o wartości 48, co chyba nie jest prawidłowym...
Podaj przykład takiego przecięcia , bo nie potrafię sobie wyobrazić na szybko jak z tego wyjdzie trapez :?: W Matlabie figury wypełnione rysuje się poleceniem fill x = [0 1 3 4]; y = [0 2 2 0] fill(x,y,'r') co daje wynik (zdjęcie niżej) Aby zakolorować to pole , musisz znać punkty przecięcia i podać wektory czyli obszar zakolorowania.
Witam Mam za zadanie narysować obiekt w 3d majac 3 wektory współrzędnych x,y,z. Np w tej postaci: x=[0,2,2,0,1] y=[0,0,2,2,1] x=[0,0,0,0,3] funkcja powinna rysować w tym przypadku ostrosłup foremny. Próbowałem użyć funkcji plot3() niestety rysuje mi ona wyłącznie linie od pierwszego do ostatniego punktu. Czy da się to jakoś inaczej zrobić w matlabie...
Witaj. Musisz stworzyć dwa wektory x i y. i Potem użyć polecenia plot(x,y,'k'). Matlab automatycznie dopisuje wartości do wektora. k to kolor linii na wykresie
Napięcie wyjściowe filtru to napięcie z dzielnika szeregowego LC, transmitancja wychodzi: \frac{1}{s^2LC+1} - oczywiście w stanie bez obciążenia. Dla wartości na schemacie kod matlaba będzie taki: [syntax=c] C=5e-3; L=200e-6; sys=tf(1, [L*C 0 1]) [/syntax] Następnie symulujesz układ za pomocą lsim, podając sys, sygnał i wektor czasowy (t) jako argumenty.
aha no tak, impedancje opisałem w matlabie, przepraszam, już dodaje: Z1=(1/((j*w*C1)+1/R1)) Z2=(1/(j*w*C2)) Z3=(R2+(1/(j*w*C3))) Z4=(j*w*L1)+ R3 Z5=(j*w*L2) w-omega To jest raczej dobrze?, Zauważ, iż masz tu na lewo od pola edycji Symbole inżynierskie i Twoje w/w zapisy powinny wyglądać tak: Z1 = 1/((j•ω•C1) + 1/R1) Z2 = 1/(j•ω•C2)...
1. Po pierwsze, to gdzie zadajesz "i"? Ja tam widzę jedynie jego użycie. 2. Jeśli "i" jest wektorem to korzystaj z operacji tablicowych, czyli np. ".*" czy ".^" a nie macierzowych ("*", "^") - efektem tego może być punkt zamiast wektora. 3. Działa Ci w ogóle "linspace"? Według dokumentacji ta funkcja ma kolejno parametry: linspace(punkt_startowy, punkt_końcowy,...
Masz rację. Może być taki kompensator. A wartości odbiornika nie mogą być podane w prost? Bo widzę, że obliczasz je później z impedancji i zamiast wartości dostajesz wektor.
Witam! Chciałabym wykorzystać do klasyfikacji obiektów sieć neuronową (nntool w matlabie). Obiektów jest sześć, i są to różne ruchy dłoni. Poszczególny ruch ma postać wektora składającego się z 10 elementów. Takich wektorów dla każdego z ruchów mam 30 (zostało wykonanych 30 pomiarów dla każdego ruchu). Dane zapisane są w notatniku, gdzie w rzędach znajdują...
A nie lepiej jako x podawać wektor z wartościami dla ktorych chcesz obliczyc ową funkcje? Najpierw zgromadzić wartości w wektorze, potem obliczać. Jako wynik dostaniesz równiez wektor, na ktorym latwo dalej operować. Eval to armata, nie wiem czy jej potrzebujesz >> help eval EVAL Execute string with MATLAB expression. EVAL(s), where s...
Oto kod programu napisany przez mojego znajomoego function [mixedVector, fadeOut, fadeIn] = crossfade(leftSig, rightSig, crossLength, nCase) % Przejście jednego sygnału w drugi za pomocą zmiany amplitud. % autor: Bernard Latanowicz latanb(at)wp.pl % nargin % leftSig: lewy wektor % rightSig: prawy wektor % crossLength:...
Witam Programuje troche (symuluje) w matlabie 6.5 silnik indukcyjny klatkowy dwufazowy przy pomocy skrypcików zawartych w plikach „plik.m” oraz „plikfun.m” umieszczonych w folderze „work”. Plik „plikfun.m” zawiera prawe strony rozwiązywanych równań: %plikfun.m %prawe strony równań function yprime=plikfun(t,y);...
"Na piechotę" możesz to zrobić przez wymnożenie widma przygotowaną maską o długości równej rozmiarowi wektora widma. Wartości elementów maski dla częstotliwości które mają zostać wycięte równe są 0, w przeciwnym wypadku 1. Trzeba tylko wiedzieć w jakim miejscu wektora zaczynają się konkretne częstotliwości. Może coś takiego będzie pomocne: Fs=25000;...
Nie rozumiem jaki jest problem z fft, polecenie robi DFT ciągu wartości. Y = fft(X) Zwraca DFT z wektora X policzone FFT. Jeżeli X jest macierzą fft jest wykonywane dla każdej z kolumn. Y = fft(X,n) Zwraca n punktową DFT, jeżeli długość X jest mniejsza niż n, wektor jest dopełniony zerami, jeżeli n jest większe niż długość X, X jest...
Udało mi się stworzyć wykres 3D (atraktor Rosslera) przy użyciu ode45. Zdołałem też zyskać przekroje tej funkcji, ale jedynie dla płaszczyzn, które są prostopadłe czy też równoległe do osi x, y czy z. Mam jednak teraz duży problem z zapisaniem warunku, gdy powierzchnia będzie ułożona w sposób dowolny. Póki co próbowałem z płaszczyzną rysowaną w oparciu...
Witam serdecznie, mam problem z wyszukiwaniem w XML-u. Mam funkcję szukającą w XML-u: Kod: function wynik = szukanieXML(x) wynikSzukania=[{}];%tw... macierzy z strukturą findLabel = x; findCbk=''; xDoc = xmlread('piosenki.xml'); allListitems = xDoc.getElementsByTagName('piosenka'... for k = 0:allListitems.getLength-1...
Witam! Jeśli dobrze rozumiem twoje intencje chcesz z wektora zawierającego zera zrobić mniejszy wektor bez elementów zerowych, czy tak ? x= find() , która wyszukuje wszystkie elementy nie zerowe w wektorze i zwraca indeksy tych elementów. wystarczy wtedy zapisać jedną linijkę kodu y=x(find(x)); Co do drugiego pytania to nie mam jeszcze pomysłu ale zobaczę...
Możesz powiedzieć jakie masz błedy i dokładnie co chcesz wiedzieć na temat tej funkcji. Na tym co przedstawiłeś widzę 1 błąd a mianowicie masz plot w pętli i tam może być problem z wyrysowaniem tego. Przenieś plot za tego end-a powinno działać. Najogólniej opisze Ci jak działa mój kod a później twój Sprawdzamy wielkość wektora wejściowego Iterujemy...
Jest tutaj dysonans nazewnictwa. Widmem sygnału zwykło się nazywać sygnał w dziedzinie częstotliwości, więc określenia "widmo zespolone" czy "widmo amplitudowe" jest czymś dla mnie dziwnym. Cóż jak zrobisz fft z sygnału to otrzymasz w wyniku wektor liczb zespolonych które możesz przedstawić na wykresie, ale to nie będzie widmo. Jeśli z takiego wektora...
Używam zmiennych double - a to chyba maksymalna precyzja jaką można uzyskac w matlabie :( Odnośnie mathcada, czy ktoś może mi napisać w jaki sposób zapisać wielomian jako transmitancję filtru cyfrowego? Bedę jeszcze próbował dzsia nauczyć się mathcada -ale nie wiem jakie będą efekty. Mam wektor, w którym kolejne elementy oznaczją wielkości przy kolenych...
no tak ale jesli dobrze rozumiem to x jest wektorem 2x1 i w "(x(1)-1)^2" bierzesz pierwszy element z wektora odejmujesz 1 i podnosisz do kwadratu, jesli dobrze rozumuje to nie ma tutaj dziedziny funkcji ktora by sie zmieniala
Analiza sygnału: x(t)=4×sin^3 8000πt×cos3000πt Zadania: Dobór parametrów próbkowania (częstotliwości oraz ilości próbek) Przeprowadzić szybką transformatę Fouriera Przekształcenie wzoru sygnału: x(t)=4×sin^3 8000πt×cos3000πt x(t)=4×sin^2 8000πt×(sin8000πt×cos3000π... x(t)=4×sin^2 8000πt×(1/2[sin(11000πt...
Witam. Mam prośbę o pomoc. Mam do rozwiązania takie zadanie: Stosując metodę stycznych oraz metodę siecznych znaleźć dodatni pierwiastek równania: x^3-0.2*x^2-1.2=0 w przedziale (1,2) z dokładnością 0,002. Porównać obie metody, która z nich zapewnia szybsza zbieżność? Ja zastosowałem taki kod: %zdefiniowanie wektora x odpowiadajacego przedzialowi od...
Witajcie, mam problem aby interpolowac funkcje w Matlab. W moim przypadku numerycznym tylko niektore wartosci funkcji odpowiadajace pewnemu wektorowi sa niezerowe, wiec normalne procedury nie polacza mi tych punktow dajac mi obwiednie tylko uwzgledniaja zerowe wartosci rowniez dajac przebieg interpolowany/aproksymowany o zanizonych wartosciach. Pamietam,...
Cześć! Mam do zrobienia projekt filtrów. Filtry już zaprojektowałem i działają, problem tylko mam z wykreśleniem modułów ich funkcji transmitancji w skali logarytmicznej. Możecie powiedzieć jak to trzeba zrobić? Filtry są tak zrobione: % Tworzymy wektor częstotliwości filtru f4=[0 2950 3000 4000 4050 fNq]; ff4=f4/fNq; % Wektor pożądanych amplitud...
Dobrze jest tak powinno być... jeden prążek to znaczy że tylko jest jedna częstotliwość składowa (jeden sinus) i ten prążek jest na wysokości tej częstotliwości... Spróbuj sobie zrobić np: Tc=1; %czas próbki Fs=100; %czestotliwosc próbkowania t=(0:1/Fs:Tc); %wektor czasu x=sin(2*pi*15*t)+sin(2*pi*40... %t wektor czasu,...
Witam. Mam nadzieję, ze nie trafilem źle, jeśli temat jest w złym dziale, to przepraszam za problemy. Usiłuję przeprowadzić analizę częstotliwościową sygnału dźwiękowego w formacie .wav. Na własną rękę, korzystając z własnej wiedzy nie za bardzo udało mi się to zrobić. Zacząłem więc korzystać z poradników i tutoriali. Jednym z nich jest ten tutaj: On...
Nie napisałeś jakie jest próbkowanie. To znaczy zakładam, że jest wielokrotnością okresu szukanej harmonicznej. Myślę, że tu trzeba by zrobić analizę próbki o długości 1s co odpowiada 1Hz. Wtedy w widmie dostaniesz prążki odsunięte od siebie o wielokrotność 1Hz i tylko trzeba znaleźć ten najwyższy. Analizujesz tylko połowę wektora wyjściowego z funkcji...
Witam. Moje pytanie pewnie dla niektórych wyda się banalne. Mam obiekt o transmitancji Go=0.1/s(2s+1)(s+1) który zamknięty jest w pętli sprzężenia zwrotnego z regulatorem P o wzmocnieniu k w taki sposób Moim zadaniem jest narysować linie pierwiastkowe tego układu przy użyciu Matlaba. Liczę transmitancję układu zamkniętego (jeśli zrobiłem to źle proszę...
Mam problem z minimalizacją pewnej funkcji. Na wstępie krótki opis zagadnienia, a potem problem z jakim się borykam. Muszą wyznaczyć 2 parametry pieca - maksymalną moc grzałki (Pmax) oraz iloczyn współczynnika wymiany ciepła i powierzchni wymiany ciepła (k*A). Model matematyczny układu: dQ/dt = Q*we - Qwy + P - kA*(Twy - Tot) W stanie ustalonym po przekształceniach:...
Bez problemu. Tworzysz sobie nowy wektor, który składa się z tych wektorów x1, x2,..., xn. Wygląda to tak: [syntax=matlab] x=[x1 x2 x3 ... xn];[/syntax] Czyli dokładnie tak, jak napisałeś. Możesz też zrobić z tego macierz, jeżeli rozdzielisz te wektory średnikami(wtedy będą one po prostu w poszczególnych wierszach), ale wtedy te wektory muszą być równej...
Witam koledzy! Mam na DSP laborkę z filtrów IIR i w jednym zadaniu mam napisać własną funkcję do filtrowania sygnału wygenerowanym filtrem IIR, coś na wzór filtfilt() tylko pewnie z większymi błędami w obliczeniach. Po wygenerowaniu filtru mam dwa wektory ze współczynnikami [L, M]. Pierwszy wektor L wystarczy pewnie tylko wymnożyć przez sygnał, ale...
Witam.Mam zrobić projekt sterowania wektorowego silnikami indukcyjnymi w programie matlab w simulinku.Siedzę nad tym problemem już dość długo i sporo rzeczy mi nie wychodzi.Jeśli istnieje taka możliwość że ktoś mógłby mi w tym pomóc (oczywiście nie za darmo),to bardzo proszę.Zależy mi na tym aby się spotkać osobiście ponieważ myślę że tak byśmy szybciej...
Witam od pewnego czasu próbuje uzyskać taki o to wykres Tak samo jak na załączonym wyżej obrazku mam widmo drgań w postaci wektorów przy różnych prędkościach obrotowych. Pozdrawiam WK Załącznik usunąłem. - arnoldziq
Możesz stworzyć wektor x z pomocą funkcji random
na drugim, niestety nie mieliśmy za wiele ;/ idee rozumiem, jednak nie wiem jak to przełożyć do matlaba, ma być to coś tego typu: W = zeros(size(Tanalogowy)); for k = 1:length(Tprobek) A = Sygnal(k) * sinc(Tanalogowy - Tprobek(k)); W = W + A; end i W to mój wynik, ale co z Tanalogowy? to ma być wektor czasów, w zakresie czasu próbkowania ale dużo gęstszy?
Czyli jednak error nie jest skalarem. Tak jest za każdym razem, czy tylko czasami? Rozumiem, że całka to powinien być wektor 1xN z którego bierzesz ostatnią wartość. Jeśli error nie jest skalarem tylko czasami, dobrze byłoby ustawić w funkcji fc breakpoint który złapie taki przypadek, tj. dodać na samym końcu if not(isscalar(error))...
Oczywiście przed napisaniem posta wygooglowałem, że najwięcej osób poleca właśnie autocada jednak spotkałem się też z opinią że narysowanie przebiegu sinusoidalnego w autocadzie jest bardzo trudne i lepiej do tego celu użyć innego programu. Na dobrą sprawę potrzebuje tylko jednego rysunku mającego na celu oddać istotę symetrii napięć w układzie trójfazowym...
Okey, właściwie to to nie jest problem, wystarczy w nawiasie napisać oba parametry oddzielając je przecinkiem czyli np [syntax=matlab] funkcja = '(at)(x_2, x_4)(exp(x_4^2)+(x_2^2-5))'; [/syntax] tylko, że teraz mam problem w rysowaniu tej funkcji. Pytanie: czy matlab udostępnia jakąś metodę do rysowania wykresu 3d, kiedy jest posiadany tylko handler...
witam, pewnie to zadanie jest już przedawnione, ale pozwoliłem sobie zapuścić na nie 'scilaba' (=matlab free) kod: //--------------------------- clear; //czyszcze wszystkie zmienne //ilość punktów do przybliżenia generowania modułu transmitancji N=128; //próbkowanie fpr=480000; //zadany wektor współczynników filtra h=[0.0717 0.0451 -0.0956 -0.2705 0.5587...
witam. Zastanawiam się nad programami umieszczonymi w tytule. Chodzi o to który program byłby lepszy i dlaczego do symulacji sterowania wektorowego silnikiem indukcyjnym. Chyba niektórzy znają te programy i mogli by się wypowiedzieć.
Witam, Mam pewien problem, program zatrzymuje się w pewnym momencie. [syntax=MATLAB]x=[1;0.5]; it=1; A=zeros(50,1); while it<50 f=[3*x(1)^2 + 4*x(2)^2 - 3; x(1)^2+x(2)^2 - sqrt(3)/2]; jac=[6*x(1) 8*x(2); 2*x(1) 2*x(2)]; dx=inv(jac)*f; x=x-dx; it=it+1; if max(abs(dx))<0.001 break; end A(it,1)=it; B(it,1)=x; end plot(A,B);[/syntax] Domyślam się...
Witam, Czy jest mi ktoś w stanie pomóc? Mam za zadanie zaprojektować filtr różniczkujący typu III o zadanej przez użytkownika długości. Zastosuj parametryczne okno Kaisera (o różnych wartościach parametru. Wykreśl: odpowiedzi impulsowe, charakterystyki amplitudowe i fazowe filtra. Wszystko, a praktycznie prawie wszystko pięknie działa tylko wychodzi...
Poprawiłem kod według wskazówek, lecz dopiero po zapisie wektorów przedzial i y0 przed „syms y(t)” błąd zniknął. Nie jestem jednak pewny czy wszystko działa poprawnie.
Mam nastepepujacy problem: Mialem za zadanie w Simulinku zaprojektowac system 2-FSK. zadnie to zrobilem. tylko oczywiscie kochany profesorek wymyslil sobie ze mam opracowac interfejs graficzny, ktory umozliwi ogladanie wykresow, wprowadzanie sygnalu nadanego i ogladanie sygnalu odebranego. w GUIDE udalo mi sie zrobic tak: stworzylem 10 Edit Text w kazdym...
Na początek to Ci się pewnie przyda (matlab, wczytanie i downmix pliku): % ----------------------------------------... function Open_Callback(hObject, eventdata, handles) % hObject handle to Open (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and...
"Solving FUN_S1 requires an initial condition vector of length 5" to oznacza, że potrzebujesz do rozwiązania swojego równania FUN_S1 wektora z warunkami początkowymi o długości 5, więc dlaczego zmniejszasz ilość warunków początkowych ? [/code]
Do wektorow to bedzie ciezko bo nie da sie zaznaczyc kolumny w pliku tekstowym. No chyba ze wszystko do excela i wtedy pojedyncze kolumny do wektora. Ale prosciej wszystko wrzucic do macierzy a potem z macierzy wydzielic wektory lub od razu dodac odpowiednie elementy macierzy. Nie bardzo tez rozumiem o co Ci chodzi. Chcesz dodac dwie pierwsze kolumny,...
patrz -> xcorr periodyczne korelacje ma pan taksamo podliczyc wedlug real(ifft(fft(a).*fft(b))) :-) Dodano po 25 przepraszam real(ifft(fft(a).*fft(b))) to jest convolucje korelacja wedlug real(ifft(fft(a).*conj(fft(b)))) Nooo... Ale wyniki będą inne dla takich samych wektorów próbek. Choćby długość wyniku będzie inna. xcorr(a,b) - długość a+b-1 real(ifft(fft(a).*fft(b)))...
Witam Czy chodzi Ci moze o cos takiego?: rand(X,X) np. rand(20,1) da wektor o dlugosci 20. Pozdrawiam
hm, trochę to strzelanie do wróbla z armaty ale może to rozwiąże problem. Spróbuj wygenerować wektor modułu widma i przefiltrować go filtrem dolnoprzepustowym (eksperymentalnie dobierając częstotliwość odcięcia).
Czyli rozumiem że zmieni się admitancja własna węzła V2 oraz coś w wektorze wymuszenia? A reszta będzie taka sama? A co do tego "pi". Chyba już taki matlaba urok że ciężko mu dojść do zera w obliczeniach bo i tak pokazuje to samo. Przykładowo dla bilansu mocy (moc źródeł minus moc odbiorników) wychodzi (-9.094947017729282e-013 -3.637978807091713e-012i)VA...
cos wyszlo ale tez nie tak jak chcialem:(..jak zastosowac okno hamminga do wektora x??.. bo tych szumow widmo mi nieuwydatnia:( clear all clf x=FFT_x0_1(:,2); x=x-mean(x); % to brakowalo Fs=100; T=1/Fs; L=6009; t=(0:L-1)*T; NFFT=2^nextpow2(L); Y=fft(x,NFFT)/L; f = Fs/2*linspace(0,1,NFFT/2); plot(f,2*abs(Y(1:NFFT/2)))
Witam Pierwsza sprawa, to po co Ci ten ramp? Masz policzyc wartosc srednia pewnej funkcji, wiec to co bedziesz mial na osi x nie ma zadnego znaczenia. Uzyj bloczku clock i nazwij go sobie x zamiast t. Czas symulacji ustaw na 12 i bedzie gralo. Druga sprawa to sposob wyliczenia wartosci sredniej: -mozesz to zrobic tak jak napisal moj poprzednik (ale...
Witam Mało co widać że zmienia się kontrast, polecam bardziej przejść z przestrzeni RGB na HSV i na nich operować. Więcej informacji co to HSV na google. A = imread('1.bmp'); n=0; C = rgb2hsv(A); %żeby przyśpieszyć obliczenia C1 = rgb2hsv(A); akcja = ['n=get(suwak,''value''); C(:,:,3) = C1(:,:,3)+n;...
lsim juz odkrylem. ale i to nie wystarczy bo tam z gory zakladam jakis wektor wymuszenia a ja mam dodac regulator wiec wymuszenie sie zmienia. jak narazie wymyslilem ze bede symulowal za pomoca tego lsim i poprostu co kazda symulacje zwiekszalem ilosc probek tzn dodawalem kolejna wartosc wymuszenia itp zawsze jakas metoda. no ale niestety nie wiem jak...