лабораторная работа / лаб.№1
.docx1.задание
>> A=[2 2 6 1; 0 -1 2 1; 4 -3 1 -5; 3 -5 1 1]
A =
2 2 6 1
0 -1 2 1
4 -3 1 -5
3 -5 1 1
>> B=[15;18;37;30]
B =
15
18
37
30
>> x=A\B
x =
-4.1579
-7.9393
7.2834
-4.5061
>> disp(B')
15 18 37 30
>> disp(B*3)
45
54
111
90
>> disp(B+3)
18
21
40
33
>> disp(B.^3)
3375
5832
50653
27000
disp(A')
2 0 4 3
2 -1 -3 -5
6 2 1 1
1 1 -5 1
>> disp(A.^3)
8 8 216 1
0 -1 8 1
64 -27 1 -125
27 -125 1 1
>> disp(inv(A))
0.2105 -0.7368 -0.0526 0.2632
0.1498 -0.4089 -0.0567 -0.0243
0.0324 0.4251 0.0688 -0.1134
0.0850 -0.2591 -0.1943 0.2024
>> disp(expm(A))
66.8127 65.8231 124.0375 -130.3746
12.8811 13.3128 24.4240 -26.0323
5.2554 26.6598 25.6335 -39.2805
47.7319 29.3321 75.4613 -69.4299
>> disp(logm(A))
0.0529 4.5687 1.4325 0.3975
-1.0682 3.1713 0.9196 0.8951
1.4749 -4.2142 0.9782 -1.0059
0.2093 -0.9425 0.7112 1.3070
>> disp(sqrtm(A))
0.8153 2.7520 1.7010 0.6151
-0.6018 2.0362 0.8349 0.7162
1.5424 -2.8966 1.1981 -1.4194
0.5980 -1.4963 0.4514 1.5312
2.задание
>> x=0:0.1:10;
>> y=cos((x.^3)+2)-exp(x);
>> plot(x,y);grid
>> x=0:0.1:1;
>> y=cos((x.^3)+2)-exp(x);
>> [TH,R]=cart2pol(x,y);
>> plot(TH,R),grid
3.задание
>> P1=[3-5];P2=[-6-1-2-8];
>> roots(P1);
>> roots(P2);
>> om0=1e-2;omk=1e2;
>> OM=logspace(-2,2,41); p=i*OM;
>> ch=polyval(P1,p);zn=polyval(P2,p);
>> ACH=abs(ch)/abs(zn);
>> FCH=angle(ch/zn)*180/pi;
>> subplot(2,1,1);
>> loglog(OM,ACH);grid
>> figure,semilogx(OM,FCH);grid