Matlab写一个小程序

%%

clc

clear

m = linspace(1,100,100)

flag = 0

res = 0

for i = 3:length(m)

for j = 2:i-1

if mod(m(i),j) == 0

flag = 1

break

else

flag = 0

end

end

if flag == 0

r = [res,m(i)]

res = r

end

end

r(2:length(r))

1、在matlab *** 作界面，word表格不可能调用，尽量用excel表格。

2、matlab使用excel表格时，可又世脊通过在界面 *** 作窗口的左上方，用创建新表格的方式创建一个空表格，再将excel表格调入即可转化为matlab可使用搜纯渗的数据矩阵。

3、从问题上看，x，y两个变量的维数应该一致，是一裤岁个n*1的n列1行的矩阵。

4、matlab关于函数式的表达，与代数式完全一样：a=((1/y)^(1/2))*x

5、如果x y两个变量的取值是正交方式，那么就要引入循环变量：for i=1:n，j=1:n，上式变为

aij==((1/yij)^(1/2))*xij

6、在matlab窗口下，可以先编程序，也可一条一条命令输入。如果是新手，建议一条一条输入，以便检查每一条的语法是否正确。

刚好这里段悔有一个现成的

%自己编程剔亩昌除数据异常值

clear off

clc

X=[1 3 5 7 4 6 3 8 4 7 12 3 5 7 17 4 6 7 8 5 3]

ls=length(X)

K=3

figure(1)

plot(X)

for i=1:3

A1(i)=abs(X(i)-(X(i+3)+X(i+2)+X(i+1))/3)

Y1=[abs(X(i+3)-X(i+2))abs(X(i+3)-X(i+1))abs(X(i+1)-X(i+2))]

Y(i)=max(Y1')

if(A1(i)/Y(i)>K)

X(i)=(X(i+1)+X(i+2)+X(i+3))/3

else

end

end

for i=ls-2:ls

A2(i)=abs(X(i)-(X(i-1)+X(i-2)+X(i-3))/3)

Z1=[abs(X(i-3)-X(i-2))abs(X(i-3)-X(i-1))abs(X(i-1)-X(i-2))]

Z(i)=max(Z1')

if(A2(i)/Z(i)>K)

X(i)=(X(i-3)+X(i-2)+X(i-1))/3

else

end

end

for i=4:ls-3

A1(i)=abs(X(i)-(X(i-3)+X(i-2)+X(i-1))/3)

A2(i)=abs(X(i)-(X(i+1)+X(i+2)+X(i+3))/3)

Y1=[abs(X(i-3)-X(i-2))abs(X(i-2)-X(i-1))abs(X(i-3)-X(i-2))]

Y(i)=max(Y1')

Z1=[abs(X(i+3)-X(i+2))abs(X(i+3)-X(i+1))abs(X(i+1)-X(i+2))]

Z(i)=max(Z1')

if(A1(i)/Y(i)>K&A2(i)/Z(i)>K)

X(i)=(X(i-3)+X(i-2)+X(i-1)+X(i+1)+X(i+2)+X(i+3))/6

else

end

end

figure(2)

plot(X)

再握耐正来一个简单的

clc

clear all

load X

K=1N=length(X)

x=X

figure(1)

plot(x)

for i=4:N-3

if (x(i)-x(i-1)>K&x(i)-x(i-2)>K&x(i)-x(i-3)>K)

x(i)=(x(i-3)+x(i-2)+x(i-1)+x(i+1)+x(i+2)+x(i+3))/6

end

end

figure(2)

plot(x)

---