Forum studentów informatyki UE Wrocław
Forum studentów informatyki UE Wroclaw

Analiza math - [2010L]MatLab

burnold - 2010-02-28, 23:09
Temat postu: [2010L]MatLab
Poszukuję osoby, która zna się na Matlabie na tyle by rozwiązać parę zadań. Poziom trudności pewnie taki jaki będzie pod koniec semestru na laborkach.

Pomoc wynagrodzę walutą studencką w postaci browarków :)

e^x - 2010-03-02, 17:50

Podaj treść zadań to może się coś wymyśli
burnold - 2010-03-05, 00:16

zadanka w załączniku...
e^x - 2010-03-05, 00:45

Zadanie nr 2:

a=[3 2 2 5; 5 -1 2 3; -1 2 3 1; -3 2 -1 0]
b=[-2 3 0 2; 4 0 -1 4; 6 0 3 5]
x=b*inv(a)

Zadanie nr 3:

b=[-2 3 0 2; 4 0 -1 4; 6 0 3 5]
>> inv(b*b')

Zadanie nr 4:

syms x
int((x^4)*(exp(-2*x)),0,2)

Zadanie nr 5:

x=[4.12:-0.09:1.06];
length(x)
suma=(((x(1)+x(35))/2)*35
syms n
ciag_arytmetyczny=4.12-(n-1)*0.09;
limit(ciag_arytmetyczny,n,inf)

(Wersja matlaba 7.7.0, w 7.0 komenda "syms x" nie działa), później może wymyślę resztę

burnold - 2010-03-05, 08:38

super, dzięki szefie!
e^x - 2010-03-05, 13:20

Później dopiszę jeszcze pierwsze zadanie jak znajdę komendy odpowiednie
e^x - 2010-03-05, 22:02

wydaje mi się ,że to może być ten program:

Utwórz nowy Blank M.file :
function[q]=qwerty(arg)
x=arg(1)
y=arg(2)
z=arg(3)
t=arg(4)
q=t^2 - 2*t*y + t*z + 2*x^2 - 2*x*y + 3*x*z - 2*x + 3*y^2 - 2*y + 2*z^2 - 4*z;
zapisz i w matlabie pisz:


fmincon('qwerty',[-inf -inf -inf -inf])
lub
fmaxsearch('qwerty',[-inf -inf -inf inf])

więc nie istnieje minimum
tak mi się wydaje, ale jak będziesz miał prawidłowy kod to chętnie bym go obejrzał

maxdiaz - 2010-03-08, 18:00

eee to ja zarzuce pytaniem z innej beczki a mianowicie
GRANICE FUNKCJI
ma ktoś coś o tym ? jakiekolwiek wytlumaczenie w .pdf albo zeby podał jakies przydklady ktore beda odpowiadały naszemu szanownemu kierownikiowi slużby laboratoryjno analizo matematycznej ?
niby mam jakies tam banały granic ale to chyba go nie zadowoli ... :-)

e^x - 2010-03-08, 21:56

>> help limit
--- help for sym/limit ---

LIMIT Limit of an expression.
LIMIT(F,x,a) takes the limit of the symbolic expression F as x -> a.
LIMIT(F,a) uses findsym(F) as the independent variable.
LIMIT(F) uses a = 0 as the limit point.
LIMIT(F,x,a,'right') or LIMIT(F,x,a,'left') specify the direction
of a one-sided limit.

Examples:
syms x a t h;

limit(sin(x)/x) returns 1
limit((x-2)/(x^2-4),2) returns 1/4
limit((1+2*t/x)^(3*x),x,inf) returns exp(6*t)
limit(1/x,x,0,'right') returns inf //granice jednostronne
limit(1/x,x,0,'left') returns -inf //w punkcie jak się domyślam
limit((sin(x+h)-sin(x))/h,h,0) returns cos(x)
v = [(1 + a/x)^x, exp(-x)];
limit(v,x,inf,'left') returns [exp(a), 0]

Reference page in Help browser
doc sym.limit

Przykłady:

>> syms x //deklaracja po czym liczymy
>> limit((1/x^2),inf) //po x dążącym do +oo

ans =

0

>> limit((1/x^2),-inf) // //po x dążącym do -oo

ans =

0

>> limit((1/x^2),0) //po x dążącym do 0

ans =

Inf

maxdiaz - 2010-03-08, 22:18

ahaaa... INF to jest ustalanie po czym dąży :| no supcio ;] dziekować M ;)
shawnzy - 2010-03-08, 22:25

maxdiaz, inf - oznacza nieskończoność (infinity).
limit((1/x^2),0) - w limicie, drugi argument określa do czego ma dążyć.

e^x - 2010-03-08, 22:46

maxdiaz, Nieprawda że
shawnzy, Prawda ,że

Kociania - 2010-05-19, 22:13

4) znaleźć rozwiązanie równania sin^2(x)*cos(x) = x*tg(x). Ile jest rozwiązań w przedziale (-2,2)
5) znaleźć takie t, które spełnia równanie: całka oznaczona od 0(dół) do t (góra) e^(-t) = 1
6) całka oznaczona od 0(dół) do t (góra) e^(-t) = 2/3. I podać ten wynik w formie liczby, ponieważ wychodzi tam logarytm.

Czy ktoś mógłby pomóc? Chodzi mi o to, jak zapisać w Matlabie "=" czyli jak zapisać, że coś ma się czemuś równać (tak jak w zadaniu 4)?

e^x - 2010-05-19, 22:29

równość w Matlabie to == ,a nierówność ~=
Kociania - 2010-05-19, 22:35

ale to oddaje tylko wartość 0 lub 1 czyli prawda lub nieprawda, ale nie oblicza wyniku
e^x - 2010-05-19, 22:38

hmm, znajdź komendę rozwiązującą układ równań
killaruna - 2010-05-19, 22:44

e^x napisa³/a:
równość w Matlabie to == ,a nierówność ~=

Jak zastosujesz tu ==, to matlab wykona test logiczny i zwróci ci prawda (jako 1) lub fałsz (jako 0), a chyba nie o to tu chodzi...

Nie wiem, na moje to bym zadeklarowała

>> sym x

a potem

>> y = solve(całe_wyrażenie)

Turinio - 2010-05-19, 22:58

Tak niestety nie działa :/
e^x - 2010-05-19, 23:03

Więc pozostaje przestudiować dwie lektury: google i help
Turinio - 2010-05-19, 23:06

e^x, wielkie dzięki za pomoc, nie wpadliśmy na to wczesniej
e^x - 2010-05-19, 23:08

Może napiszę coś później jak skończę własny skrypt w Matlabie
pit3r - 2010-05-19, 23:10

Turinio, nie obrażaj się na e^x :-P
W niektórych przypadkach to najlepsze wyjście po prostu ;)

Kociania - 2010-05-19, 23:15

Tylko, że my z Turinio, spędziliśmy nad tym sporo czasu :/ i na takie pomysły jak google itd wpadliśmy już na początku.
tomx2 - 2010-05-19, 23:34

x=solve('sin(x)*sin(x)*cos(x)-x*tan(x)=0') ?? nie wiem :shock:
Kociania - 2010-05-19, 23:36

bez tego "=0"
tomx2 - 2010-05-19, 23:37

to w sumie bez znaczenia czy tamto czy x=solve('sin(x)*sin(x)*cos(x)=x*tan(x)') , bo to to samo przeciez, stąd wychodzi że x=0 , ale czy to prawda?
Kociania - 2010-05-19, 23:39

chciałabym wiedzieć -.-
pit3r - 2010-05-19, 23:41

a.zietkowska napisa³/a:
Tylko, że my z Turinio, spędziliśmy nad tym sporo czasu :/ i na takie pomysły jak google itd wpadliśmy już na początku.


no cóż metoda prób i błędów pozostaje w takim razie ;-)

Kociania - 2010-05-19, 23:42

Na jutrzejszym kole to raczej nie będzie czasu an próby,a tym bardziej na błędy :/
tomx2 - 2010-05-19, 23:43

a.zietkowska, chociaż ta metoda wydaje się dobra , bo jak wstawia się inne wartości do tego równania to wylicza całkiem niegłupie rzeczy :idea: jeszcze 2 kolejne zadania i koniec! :-)
tomx2 - 2010-05-20, 00:35

zrobił ktoś 2 ostatnie zadania? wyszlo mi w 5) -0,693 a w 6) -0,5108 .
shawnzy - 2010-05-20, 08:26

tomx2, zad5 masz źle bo ma wyjść brak rozwiązań.
ans =

[ empty sym ]

zad6.
log(3) = 1.0986

denis - 2010-05-20, 08:58

grupa a

>> inv(A)

ans =

-18.0000 14.0000 20.0000 -5.0000
4.0000 -3.0000 -4.0000 1.0000
11.0000 -8.0000 -12.0000 3.0000
4.0000 -3.0000 -5.0000 1.0000

A^2 - (A * A')

ans =

-3 0 3 10
-1 -1 0 -11
0 6 -2 6
8 21 -19 -43

denis - 2010-05-20, 08:58

>> det(A)

ans =

-1

>> det(A)^17

ans =

-1

denis - 2010-05-20, 09:04

>> inv(A)*B

ans =

2.0000 -0.0000
1.0000 1.0000
-1.0000 2.0000
0.0000 5.0000

to jest podpunkt B


Powered by phpBB modified by Przemo © 2011 phpBB Group