disp('Câu 1');
A=[0 1 0; 0 0 1; -4 5 8]
B=[0; 0; 7]
C=[9.5 0 0]
I=[1 0 0; 0 1 0; 0 0 1]
syms s;
disp('a');
W = C*((s*I-A)^(-1))*B;
disp('W= ');
disp(W);
disp('b');
TBC=mean(W(:));
disp('TBC W= ');
disp(TBC);
disp('c');
D=sum(A(:));
disp('D= ');
disp(D);
disp('d');
minA=min(A(:));
disp('min =');
disp(minA);
%------------------------------
disp('Câu 2');
t = 0:20;
x = 0:2*pi;
w = 0:pi;
f1=8*(1-exp((-1/7)*t));
f2=tan(sqrt(x));
f3=cos(w)+4+j*sin(w);
figure;
subplot(2,2,1);
plot(t, f1, 'b-', 'LineWidth', 2);
grid on;
yticks([0 2 4 6]);
xlabel('x');
ylabel('y');
title('f1');
ylim([0 6]);
xticks([ 0 5 10 ]);
xlim([ 0 10]);
subplot(2,2,2);
plot(x, f2, 'r-', 'LineWidth', 2);
grid on;
yticks ([ -1 0 1 2]);
xlabel('x');
ylabel('y');
title('f2');
ylim([ -1 2]);
xticks([ 0 2 4 6 8]);
xlim([ 0 8 ]);
subplot(3,1,3);
plot(real(f3), imag(f3), 'g-', 'LineWidth', 2);
grid on;
xlabel('x');
ylabel('y');
title('f3');
yticks([ 0 1 2 3 ]);
ylim([0 3]);
xticks([ 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5]);