Question

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

Answer 1

%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<num<=50
xx2=f2*xx02+l2*sqrt(Q2); %目标的状态值
XX2(:,num-25)=xx2;
z(num)=h*xx2+sqrt(R)*randn(1); %目标的观测值
else
if 50<num<=timenum
xx3=f1*xx01+l1*sqrt(Q1); %目标的状态值
XX3(:,num-50)=xx3;
z(num)=h*xx3+sqrt(R)*randn(1); %目标的观测值
end;
end;
end;

追问

大哥 没有for啊 而且这是什么功能？？

追答

for可以引用的，这是元胞工程