| |
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
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
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 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
|
|