[Sidney1]

Koledzy i Koleżanki z #matlab trochę może #matematyka i #fizyka

wyrzuca mi w matlabie błąd jak niżej, jak zrobię to dla stałej wartości t wszystko jest ok, a muszę zrobić plot, jakieś pomysły?
Error using ^ (line 51)
Incorrect dimensions for raising a matrix to a power. Check that the matrix is square and the power is
a scalar. To perform elementwise matrix powers, use '.^'.

Error in Untitled (line 6)
u=(1000/kb)*(((t/ta1)^(nb))/(1+(t/ta1)^(nb)))*exp((-t/ta2));

kod:
ta1=0.0000035;
ta2=0.000051;
nb=1.8;
kb=exp((-ta1/ta2)*(-nb*ta2/ta1)^(1/nb));
t =0.0000000001:0.00000001 : 0.0001;
u=(1000/kb)*(((t/ta1)^(nb))/(1+(t/ta1)^(nb)))*exp((-t/ta2));

plot(t,abs(u));

[motoinzyniere]

@Sidney1: To perform elementwise matrix powers, use '.^' ( ͡º ͜ʖ͡º)

[Sidney1]

@motoinzyniere próbowałem wychodzą dziwne wyniki :p

[kolnay1]

(t/ta1)^(nb)

@Sidney1: to nie ma sensu, bo nie jest kwadratowe (jest takich wymiarów jak t)

[Sidney1]

@kolnay1 chyba nie rozumiem

[kolnay1]

@Sidney1: co to znaczy podnosić do potęgi?

[Sidney1]

@kolnay1 czy Tobie w poprzedniej wiadomości chodzi o to, że to macierz? Ja właśnie nie chce, żeby to była macierz tylko plot z jakaś tak gęstościa próbek

[kolnay1]

@Sidney1: ^ to potęgowanie macierzowe, generalnie w matlabie wiele obiektów często okazuje się być macierzami ( ͡° ͜ʖ ͡°)

[Sidney1]

@kolnay1 czyli potegowac przez power() czy jakoś inaczej zapisać t wartość t, bo rozumiem, że traktuje mi to jako wektor?

[kolnay1]

@Sidney1: t po prostu jest wektorem. Podejrzewam, że chcesz potęgować przez .^, a to że wyniki wychodzą dziwne wynika stąd, że coś innego się tu popieprzyło

[Sidney1]

@kolnay1 nie chce, żeby był wektorem. Chce wykreślić funkcję od 0 do jakiś 0.001 t to argument funkcji.

[kolnay1]

@Sidney1: Jak rozumiesz słowo wektor? Czego oczekujesz po takim napisie:

t =0.0000000001:0.00000001 : 0.0001;

?

[Sidney1]

@kolnay1 coś takiego chciałem zrobić tylko dla innej funkcji, zapis chyba taki sam

[kolnay1]

@Sidney1: nie no rozumiem, tylko że w tym wypadku x jest wektorem, tzn macierzą o wymiarach 1 x n. Ale mniejsza o to, bo jak w tym co wrzuciłeś dopiszę znaki mnożenia (które jak podejrzewam wykop zjadł) i zamiast ^ napiszę .^ (czyli potęgowanie wyraz po wyrazie, a nie macierzowe), to wykres wygląda chyba dobrze, tj. reprezentuję to czego spodziewałem się po wzorze.

Spróbuj na czymś prostszym np. na funkcji y=x^2

[Sidney1]

@kolnay1: ok rozumiem różnię, ale wydaje mi się, że z tym wektorem t jest coś nie tak. Według mnie dla 0 wartość powinna być 0, on wykreśla już jakieś spore wartości. W konsoli natomiast nieważne jak gęsto ustawię wektor pokazuje mi tylko 100 kolumn.

[Sidney1]

@kolnay1: a tak w zasadzie to chodzi mi o to, żeby wykreślić coś takiego

[kolnay1]

@Sidney1: hej, trochę o tym myślałem, ale póki co nic sensownego nie wymyśliłem. Dam znać jakby coś się zmieniło.

[Sidney1]

@kolnay1: z tym sobie poradziłem, chodziło to, żeby kropka była przed każdym znakiem ze zmienną (potęga, dzielenie i mnożenie). Ale przy okazji mam następny problem ( ͡° ͜ʖ ͡°) robię z tego FFT i to nawet wychodzi, ale jak chce to fft wyplotować skala nie zależy od częstotliwości jaką ustawie. W sensie pokazuje mi wzmocnienie i przesunięcie, ale nie pokazuje dla jakiej częstotliwości.

[kolnay1]

@Sidney1: no tak, ja wszędzie dałem tę kropkę. To pewnie dlatego ciebie dziwiło to co wychodziło, a mi wychodziło w na oko dobrze ( ͡° ͜ʖ ͡°)