不少于30行的matlab程序设计,急需一道matlab程序设计题 不少于30行 功能不限 只求注明 包含for和if即可!!!!...

满意答案

hanshugen

2014.02.07

采纳率：44%    等级：12

已帮助：6256人

%CA setup

n=200;

%initialize the arrays

z=zeros(n,n);

cells=z;

sum=z;

%set a few cells to one

cells(100,100)=1;

%index definition for cell update

x=2:1:199;

y=2:1:199;

((x.^2)+(y.^2))<=100^2;

%nearest neighbor sum

sum(x,y)=cells(x,y+1)+cells(x,y-1)+cells(x+1,y)+cells(x-1,y);

% The CA rule

cells=(sum==1);

plot(x,y);

function imm_extend2() %交互式多模型 初步改进

clear all;

close all;

timenum=75;

t=1;

transfer_matric=[0.98 0.02; %马尔科夫转移矩阵 矩阵为j行i列

0.02 0.98];

u_10=0.5; %目标1在模型i在k-1时刻的概率

u_20=0.5; %目标1在模型i在k-1时刻的概率

u=[u_10 u_20];

Q1=100; %目标1的状态噪声协方差阵

Q2=4000; %目标1的状态噪声协方差阵

R=10000; %观测噪声协方差阵

%R2=700;

x00_1=[0 0 0]'; %模型1估计的初始值

x00_2=[0 0 0]'; %模型2估计的初始值

p00_1=[100 0 0;

0 10 0;

0 0 0];

p00_2=[100 0 0;

0 10 0;

0 0 10];

dis=1000; %位移

vel=450; %速度

acc=50; %加速度

xx01=[dis vel 0]'; %状态的初始值

xx02=[dis vel acc]'; %状态的初始值

%xx2=0;

f1=[1 t 0;

0 1 0;

0 0 0]; %状态转移阵

f2=[1 t 0.5*t^2;

0 1 t;

0 0 1];

l1=[0.5*t^2 t 0]'; %状态干扰阵

l2=[0.5*t^2 t 1]';

h=[1 0 0]; %观测转移阵

I=[1 0 0;

0 1 0;

0 0 1];

num=1;

XX1=zeros(3,25);

XX2=XX1;

XX3=XX1;

while num<=timenum

if num<=25

xx1=f1*xx01+l1*sqrt(Q1); %目标的状态值

XX1(:,num)=xx1;

z(num)=h*xx1+sqrt(R)*randn(1); %目标的观测值

else

if 25

xx2=f2*xx02+l2*sqrt(Q2); %目标的状态值

XX2(:,num-25)=xx2;

z(num)=h*xx2+sqrt(R)*randn(1); %目标的观测值

else

if 50

xx3=f1*xx01+l1*sqrt(Q1); %目标的状态值

XX3(:,num-50)=xx3;

z(num)=h*xx3+sqrt(R)*randn(1); %目标的观测值

end;

end;

end;

00分享举报