求助 Matlab函数参数设置问题
在.m文件中如何将F=中的第二个数设为参数，写成F=或F=在VC中都调用失败，不要参数一可以成功，全部代码如下：（因为是毕业设计，恳请大家帮忙，导师要求将这个设为参数在VC中为其赋值，QQ：350892132 谢谢！）
%本程序基于P68(卡尔曼滤波与维纳滤波)
function kalman1(par1,par2)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%_____________________________________________________________
%变量的设置
ST=par1;%存储采样周期
s=50;%控制采样数据点
XG=zeros(2,s);%存储最佳估计值
P=zeros(2,2);%滤波误差方差
Z=zeros(1,s);%存储观测值
Z1=zeros(1,s);%存储观测值
K=zeros(2,s);%存储最优滤波增益
XY=zeros(2,1);%存储最优预测值
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%____________________________________________________________
%常量的初始化
F=;%状态转移矩阵
H=; %观测矩阵
XG(:,1)=';%初始化最佳估计值
P(:,:)=;%初始化滤波误差方差
T=;%
B=';%
%K(:,1)=';%最优滤波增益初始化
E=0;%差值XY(2,2)=XG(:,1);%;%
XY(:,1)=XG(:,1);%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%____________________________________________________________
%产生观测值
t=zeros(50);
l=zeros(50);
Zs=par2;
for k=1:s
l(k)=1+1.5*k; % 用于存储理论值
t(k)=Zs*randn; % 用于存储模拟观测误差
Z(1,k)=1+1.5*k+t(k);%模拟观测值(为了与第５８行对应统一，所以此处设为k+1)
Z1(1,k)=Z(1,k);
end
%Z1(1,20)=-100;
%Z1(1,30)=150;
%Z1(1,40)=-50;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%____________________________________________________________
%产生噪声
Q=zeros(1,s);%存储观测误差方差矩阵
R=zeros(1,s);%存储输入控制方差矩阵
for k=1:s
R(:,k)=Zs*Zs;%观测误差方差矩阵
Q(:,k)=0.1;%输入控制方差矩阵
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%____________________________________________________________
%滤波过程
for k=1:s
Kp=F*K(:,k);%
XY(:,:)=F*XY(:,:)+Kp*E;%最优预测值(注意此处的设计)
E=Z(1,k)-H*XY(:,:);%差值
PP=F*P(:,:)*F'+T*Q(:,k)*T';% 预测误差协方差
K(:,k+1)=PP*H'*inv(H*PP*H'+R(k));%最优滤波增益
P=(eye(2)-K(:,k+1)*H)*PP;%滤波误差方差
XG(:,k+1)=F*XG(:,k)+K(:,k+1)*E;%滤波值
p11(k)=P(1,1);%存储滤波方差的变化（位置）
p22(k)=P(2,2);%存储滤波方差的变化（速度）
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%____________________________________________________________
%显示设置
m=zeros(1,s);%存储横坐标
y=zeros(1,s);%存储观测值
z=zeros(1,s);%存储滤波值
z1=zeros(1,s);
v=zeros(1,s);
for i=1:s
m(i)=i;%产生横坐标
y(1,i)=Z1(1,i);%得到观测值
z(1,i)=XG(1,i+1);%得到滤波值
v(1,i)=XG(2,i+1);
z1(1,i)=0.15;
end
%%%%%%%%%%%______________________________________________________
%%%%%%%%%%%%%信噪比的变化对滤波的影响
Zp=0;%用于存储理论值的平均值
for i=1:s
Zp=l(i)+Zp;
end
Zp=Zp/s;%理论值的平均值
Zx=0;%用于存储信号
for i=1:s
Zx=(l(i)-Zp)*(l(i)-Zp)+Zx;
end
Zx=Zx/s;%信号
xzb=Zx/(Zs*Zs)%信噪比
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%综合显示
subplot(221);
plot(m,y,'k',m,z,'r',m,l,'b');
xlabel('时间(Time)/秒(Sec)'),ylabel('观测值(Measurement)')
subplot(222);
plot(m,z,'r');
xlabel('时间(Time)/秒(Sec)'),ylabel('滤波值(Estimates)')
axis();
subplot(223);
plot(m,v,'k',m,z1,'r');
xlabel('时间(Time)/秒(Sec)'),ylabel('速度(V)/秒(Sec)');
S=zeros(50);
O=zeros(50);
for i=1:s
S(i)=y(i)-l(i);%存储观测误差
O(i)=z(i)-l(i);%存储滤波误差
end
% subplot(234);
% plot(m,S,'k');
% xlabel('时间(Time)/秒(Sec)'),ylabel('测量误差(Measurement Error)');
subplot(224);
plot(m,O,'M');
xlabel('时间(Time)/秒(Sec)'),ylabel('滤波误差(Estimate Eror)');
[ 本帖最后由 feng0830 于 2009-4-22 21:21 编辑 ] pam1 是par1？？？pam2是par2？？要不然你程序怎么可能运行得了。。
回复 沙发 ydlcsu 的帖子
是的，该怎么写？ 那就用par1和par2啊，还整个pam1和pam2出来干嘛？ 不好意思，打错了，是par，在Matlab中参数设置成那样没问题，但是在VC中调用是参数1设成那样无法画图帖子已经改了，麻烦在看下什么问题，.m文件没问题，就是参数一不晓得该怎么写
[ 本帖最后由 ChaChing 于 2009-4-22 22:45 编辑 ]
页:
[1]