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

[一道高中数学类比推理题急需帮助]f(x)+f(1-x)=1/(2^x-√2)+1/(2^(1-x)-√2)=1/(2^x-√2)+1/[(2/2^)-√2]=1/(2^x-√2)+2^x/(2-2^x*√2)=√2/(2^x*√2-2)+2^x/(2-2^x*√2)=√2/(2^x*√2-2)-2^x/(2^x*√2-2)=（...+阅读

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

%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;

两道MATLAB程序设计题

乘着今天没有事，给你做了做，可以运行：

%%%%第一题%%%%%%%%%%

a(1:1000)=1； %最多可以输入1000个数字

for i=1:1000

disp（'请输入数字后再回车，不要直接回车或输入字符'）；

a(i)=input（'请输入数字，输入0表示结束输入'）；

if a(i)==0

break;

end

if a(1)==0

disp('ERROR！没有输入数字')；

else

disp（'输入的数字为'）；

a(1:i-1)

if i==2

disp（'输入的数字的绝对值为'）

abs(a(1)） %一个数时绝对值

else if i==3

disp（'输入两个数的平均值为'）；

mean(a(1:2)） %两个数时平均值

else if i==4||i==5

disp（'输入的数的和为'）；

sum(a(1:i-1)） %三或四个数时和

else

disp（'输入的数排序为'）；

sort(a(1:i-1)） %五个以上排序

end

%%%%第二题%%%%%%%%%%

x=floor(sqrt(1890)):floor(sqrt(1890+150)）；%假设长寿老人不超过150岁

for i=1:length(x)

if x(i)^2-x(i)>1889&x(i)^2-x(i)disp（'老人出生年份：'）；

x(i)^2-x(i)

break;

end