针对某个高等数学问题,用matlab编程解决,要求有调用函数。

2024-12-02 03:55:25
推荐回答(2个)
回答1:

线性代数问题也可以吧? 求方程组的通解 x1+x2-3x3-x4=1, 3x1-x2-3x3+4x4=4, x1+5x2-9x3-8x4=0 解: A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; b=[1 4 0]’; B=[A b]; n=4; R_A=rank(A) R_B=rank(B) format rat if R_A==R_B&R_A==n X=A\b elseif R_A==R_B&R_A X=A\b C=null(A,’r’) else X=’Equation has no solves’ end 结果为: R_A = 2 R_B = 2 Warning: Rank deficient, rank = 2, tol = 8.8373e-015. X = 0 0 -8/15 3/5 C = 3/2 -3/4 3/2 7/4 1 0 0 1 所以原方程组的通解为 (图)



回答2:

clear
N=200;
w(1)=0;
w=randn(1,N)
x(1)=0;
a=1;
for k=2:N;
x(k)=a*x(k-1)+w(k-1);
end

V=randn(1,N);
q1=std(V);
Rvv=q1.^2;
q2=std(x);
Rxx=q2.^2;
q3=std(w);
Rww=q3.^2;
c=0.2;
Y=c*x+V;

p(1)=0;
s(1)=0;
for t=2:N;
p1(t)=a.^2*p(t-1)+Rww;
b(t)=c*p1(t)/(c.^2*p1(t)+Rvv);
s(t)=a*s(t-1)+b(t)*(Y(t)-a*c*s(t-1));
p(t)=p1(t)-c*b(t)*p1(t);
end

t=1:N;
plot(t,s,'r',t,Y,'g',t,x,'b');
卡尔曼滤波程序
咋样?

什么意思啊?
那我给你编写一个矩阵运算的函数,你自己加一个function名字不久完了啊