运用以下代码实现输入和输出功能：function [T,D] = quadraticForm(A) % 输入对称矩阵A，输出矩阵T和元素值为0、1、-1的对角矩阵D，实现T'AT=D为对角阵 n = size(A,1); T = eye(n); D = zeros(n); for k = 1:n-1 % 判断是否需要进行变换 if A(k+1,k) ~= 0 % 计算变换矩阵 [c,s] = givens(A(k,k),A(k+1,k)); G = [c s; -s c]; % 更新A和T A([k k+1],k:n) = G'*A([k k+1],k:n); T(:,[k k+1]) = T(:,[k k+1])*G; end end for i = 1:n D(i,i) = sign(A(i,i)); if D(i,i) == 0 D(i,i) = 1; end end T = T'; D = D*diag(abs(diag(A))); end function [c,s] = givens(a,b) % 计算Givens变换矩阵的元素c和s if b == 0 c = 1; s = 0; elseif abs(b) > abs(a) t = -a/b; s = 1/sqrt(1+t^2); c = s*t; else t = -b/a; c = 1/sqrt(1+t^2); s = c*t; end end

你可以将这个代码保存到一个.m文件中，例如quadraticForm.m。然后在MATLAB命令窗口中输入以下命令：

A = [1 2 3; 2 5 6; 3 6 9];
[T,D] = quadraticForm(A);
disp("T = ");
disp(T);
disp("D = ");
disp(D);

这个例子中，我们输入了一个对称矩阵A，然后调用quadraticForm函数计算出T和D，并使用disp函数将它们输出到命令窗口。你可以修改A的值和输出方式来适应你的需求。