MATLAB程式简化 (用矩证方式)p1=0.8; p2=0.17; p3=0.01r0 = (3+p1-p2+p3)/4s0 = (5+3*p1+p2-p3)/2*h^2s1 = (3*p3-p1+p2-3)/2*h^2输入讯号:y(k+2)=2y(k+1)-y(k)+h^2/2*u(k+1)+h^2/2*u(k)输出讯号:u(k=1)=-r1*u(k)+t0*1-s0*y(k+1)-s1*y(k)h跟k 自己
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/06 05:28:01
MATLAB程式简化 (用矩证方式)p1=0.8; p2=0.17; p3=0.01r0 = (3+p1-p2+p3)/4s0 = (5+3*p1+p2-p3)/2*h^2s1 = (3*p3-p1+p2-3)/2*h^2输入讯号:y(k+2)=2y(k+1)-y(k)+h^2/2*u(k+1)+h^2/2*u(k)输出讯号:u(k=1)=-r1*u(k)+t0*1-s0*y(k+1)-s1*y(k)h跟k 自己
MATLAB程式简化 (用矩证方式)
p1=0.8; p2=0.17; p3=0.01
r0 = (3+p1-p2+p3)/4
s0 = (5+3*p1+p2-p3)/2*h^2
s1 = (3*p3-p1+p2-3)/2*h^2
输入讯号:y(k+2)=2y(k+1)-y(k)+h^2/2*u(k+1)+h^2/2*u(k)
输出讯号:u(k=1)=-r1*u(k)+t0*1-s0*y(k+1)-s1*y(k)
h跟k 自己取值
以下是自己打的程式
clear;clc;
t0=0.1; h=0.5;
y(1)=0; y(2)=0;
u(1)=0; u(2)=0;
for k=1:500;
p1=0.8; p2=0.17; p3=0.01;
r0 = (3+p1-p2+p3)/4;
s0 = (5+3*p1+p2-p3)/2*h^2;
s1 = (3*p3-p1+p2-3)/2*h^2;
y(k+2)=2*y(k+1)-y(k)+h^2/2*u(k+1)+h^2/2*u(k);
u(k+2)=-r0*u(k+1)+t0*1-s0*y(k+2)-s1*y(k+1);
end
plot(y);
听说用矩正来想会更棒 想知道矩证方式该怎麼写?
MATLAB程式简化 (用矩证方式)p1=0.8; p2=0.17; p3=0.01r0 = (3+p1-p2+p3)/4s0 = (5+3*p1+p2-p3)/2*h^2s1 = (3*p3-p1+p2-3)/2*h^2输入讯号:y(k+2)=2y(k+1)-y(k)+h^2/2*u(k+1)+h^2/2*u(k)输出讯号:u(k=1)=-r1*u(k)+t0*1-s0*y(k+1)-s1*y(k)h跟k 自己
下一个计算要用到上一个计算的结果,不能简化