Упражнение 5. Matlab - Revision history

Anko: /* Eigen values lambda = eig(A) */

2011-04-19T14:01:16Z

Eigen values lambda = eig(A)

2011-04-19T13:23:52Z

Eigen values lambda = eig(A)

2011-04-19T13:21:08Z

Eigen values lambda = eig(A)

2011-04-19T13:05:35Z

Eigen values lambda = eig(A)

2011-04-19T13:04:41Z

Eigen values lambda = eig(A)

2011-04-19T13:04:17Z

Eigen values

2011-04-19T12:58:42Z

Решаване системи линейни уравнения

2011-04-19T12:52:18Z

Решаване системи линейни уравнения

2011-04-19T12:41:52Z

Преобразувания на Лаплас

2011-04-19T12:39:08Z

Преобразувания на Лаплас

← Older revision		Revision as of 14:01, 19 April 2011
Line 204:		Line 204:
	-0.0355		-0.0355

		+	== Корени на алгебрични уравнения fzero() ==

		+	f(x) = 0;
		+
		+	x = fzero(@sin, 3) % 3 e началното приближение
		+	x =
		+	3.1416
		+
		+	x - tg(x) = 0 % da se namerqt pyrvite tri polozhitelni korena
		+
		+	x = -5 : 0.01 : 20; plot(x, (x-tan(x)));
		+	x = -5 : 0.01 : 20; plot(x, (x.*cos(x)-sin(x)));
		+	>> x1 = fzero('x.*cos(x)-sin(x)', [4 5])
		+	x1 =
		+	4.4934
		+
		+	f(x) - kato string
		+
		+	f = inline('x.*cos(x)-sin(x)')
		+	f =
		+	Inline function:
		+	f(x) = x.*cos(x)-sin(x)
		+
		+	== Минимум на една функция ==
		+
		+	fminbnd('x.*cos(x)-sin(x)', -10, 10)
		+	ans =
		+	3.1416
		+
		+	== Числено пресмятане на определени интеграли quad() ==
		+
		+	quad(@function(x), a, b);
		+
		+	quad(@function(x), a, b, eps); % tochnost
		+
		+	quad('1./(x.^3+5)', -0, 2, 1.e-15)

	[[Category:Matlab]]		[[Category:Matlab]]

@@ Line 155: / Line 155: @@
 трептения на системи с няколко степени на свобода
+ A*x = lam*x
-A*x = lam*x
+ (A-lam*I)*x = 0 % homogenna
-(A-lam*I)*x = 0 % homogenna
+ det(A - lam*I) = 0 % togava ima resheniq A[n][m]
+ a[1]lam^n + a[2]lam^(n-1) ... = 0 % edin polinom ot n stepen ima n re[eniq
-det(A - lam*I) = 0 % togava ima resheniq A[n][m]
+ lam[1], lam[2] ... - собствени стойности на матрицата А
 за всяка собствена стойност има собствен вектор от решения

@@ Line 173: / Line 173: @@
 '''Обощена задача на собствените стойности '''
-A*x = lam*B*x % Anxn, Bnxn
+ A*x = lam*B*x % Anxn, Bnxn
-(A-lam*B)x = 0
+ (A-lam*B)x = 0
-det(A-lam*B) = 0
+ det(A-lam*B) = 0
 [[Category:Matlab]]

@@ Line 150: / Line 150: @@
 .666666666666667
    -2.000000000000000
 трептения на системи с няколко степени на свобода
 [[Category:Matlab]]

← Older revision		Revision as of 12:52, 19 April 2011
Line 112:		Line 112:

	== Решаване системи линейни уравнения ==		== Решаване системи линейни уравнения ==
		+
		+	A*x= B
		+	x = A\B % \ може да се прилага при правоъгълни матрици, при повече редове от колони се решава за средно квадратичната стойност
		+
		+	A = [3 7 5; 2 6 9; 4 4 4]
		+	A =
		+
		+	3 7 5
		+	2 6 9
		+	4 4 4
		+	>> b = [2 3 4]'
		+	b =
		+	2
		+	3
		+	4
		+	>> A\b
		+	ans =
		+	1.0500
		+	-0.4500
		+	0.4000
		+
		+	При лошо обоснована матрица може да има грешки
		+	>> A = [1 2 3; 4 5 6; 7 8 9]
		+	A =
		+	1 2 3
		+	4 5 6
		+	7 8 9
		+	>> det(A)
		+	ans =
		+	0 % -> лошо обоснована матрица
		+
		+	>> x = A\b
		+	Warning: Matrix is close to singular or badly scaled.
		+	Results may be inaccurate. RCOND = 1.541976e-18.
		+	x =
		+	-3.333333333333333
		+	5.666666666666667
		+	-2.000000000000000
		+
		+	трептения на системи с няколко степени на свобода

	[[Category:Matlab]]		[[Category:Matlab]]

@@ Line 96: / Line 96: @@
   y = ilaplace(2*p/((p^2+p+1)*((p+1)^2+4)), p, x)
   plot(-13:0.01:3,subs(y,x, -13:0.01:3))
 Проверка
@@ Line 107: / Line 108: @@
   ans =
   sin(2*x)/exp(x)
- [[Image:LPlot.png]]
 [[Category:Matlab]]

@@ Line 65: / Line 65: @@
   fi =
   x*exp(x)
 [[Category:Matlab]]