当前位置：文档之家› 算法程序

算法程序

clear

clc

N=500;

M=300;

u=0.1;

n=1:(N+1);

a=[-0.98 0.98]

for k=1:2

for j=1:M

vn=randn(1,N);

a1=[1 a(k)];

b1=1;

xn=filter(b1,a1,vn);

d=var(xn);

xn=xn/sqrt(d);

fn(1)=xn(1);

wn(1)=0;

wn(2)=wn(1);

for i=2:length(xn)

fn(i)=xn(i)-wn(i)*xn(i-1);

wn(i+1)=wn(i)+u*xn(i-1)*fn(i);

end;

wm(j,:)=wn;

end;

w(k,:)=wn;

ew(k,:)=mean(wm);

end;

plot(n,w(1,:),n,ew(1,:),n,w(2,:),n,ew(2,:));

clear

clc

N=500;

M=300;

u=0.2;

n=1:N;

a=[0.98 -0.98]

k=1;

for j=1:M

vn=randn(1,N);

a1=[1 a(k)];

b1=1;

xn=filter(b1,a1,vn);

d=var(xn);

xn=xn/sqrt(d);

fn(1)=xn(1);

wn(1)=0;

wn(2)=wn(1);

for i=2:length(xn)

fn(i)=xn(i)-wn(i)*xn(i-1);

wn(i+1)=wn(i)+u*xn(i-1)*fn(i);

f2(i)=fn(i)^2;

end;

fm(j,:)=f2;

end;

fw=mean(fm);

semilogy(n,f2,n,fw);

clear

clc

N=500;

M=300;

u=[0.1 0.05 0.02];

n=1:N;

a=[-0.99 0.99]

k=1;

for l=1:3

for j=1:M

vn=randn(1,N);

a1=[1 a(k)];

b1=1;

xn=filter(b1,a1,vn);

d=var(xn);

xn=xn/sqrt(d);

fn(1)=xn(1);

wn(1)=0;

wn(2)=wn(1);

for i=2:length(xn)

fn(i)=xn(i)-wn(i)*xn(i-1);

wn(i+1)=wn(i)+u(l)*xn(i-1)*fn(i); f2(i)=fn(i)^2;

end;

fm(j,:)=f2;

end;

fw(l,:)=mean(fm);

end;

semilogy(n,fw(1,:),n,fw(2,:),n,fw(3,:));

clear

N = 100;

clear

N = 1000;

t = 0:10/N:10-10/N;

zt = 60*(pdf('Normal', t-8, 0, 0.3) + pdf('Normal', t, 0, 0.8))+rand(1,1000)*3;

figure(1);

x = 0:100/1000:100-100/1000;

plot(x,zt);

[X,Y]=meshgrid(-5:0.1:10,-5:0.1:10);

xy = [X(:) Y(:)];

p=mvnpdf(xy,[],[]); P=reshape(p,size(X));

[X,Y]=meshgrid(-10:0.1:5,-5:0.1:10);

xy = [X(:) Y(:)];

p1=mvnpdf(xy,[],[]); P1=reshape(p1,size(X));

X = 0:160/100:160;

Y = 0:160/150:160;

PS = P/3+P1;

PS = (PS(:,50:150) + randn(151,101)/500)*600;

figure(2);

surf(X,Y,PS);

axis tight

能量检测%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%生成调制信号

fs=100;

%采样频率

fc=40;

%载频

fo=fs/20;

%码率

L=300;

%信号样本

t = (0:L-1)*1/fs;

xn=cos(8*pi*fc*t);

y = xn + randn(size(t));

figure(1)

plot(t,y) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%进行能量检测

NFFT = 2^nextpow2(L);

Y = fft(y,NFFT)/L;

f = fs/2*linspace(0,1,NFFT/2);

figure(2)

plot(f,2*abs(Y(1:NFFT/2)))

程序代码相当全面

clear

N = 100;

clear

N = 1000;

t = 0:10/N:10-10/N;

zt = 60*(pdf('Normal', t-8, 0, 0.3) + pdf('Normal', t, 0, 0.8))+rand(1,1000)*3;

figure(1);

x = 0:100/1000:100-100/1000;

plot(x,zt);

[X,Y]=meshgrid(-5:0.1:10,-5:0.1:10);

xy = [X(:) Y(:)];

p=mvnpdf(xy,[],[]); P=reshape(p,size(X));

[X,Y]=meshgrid(-10:0.1:5,-5:0.1:10);

xy = [X(:) Y(:)];

p1=mvnpdf(xy,[],[]); P1=reshape(p1,size(X));

X = 0:160/100:160;

Y = 0:160/150:160;

PS = P/3+P1;

PS = (PS(:,50:150) + randn(151,101)/500)*600;

figure(2);

surf(X,Y,PS);

axis tight

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%生成调制信号

fs=100;

%采样频率

fc=40;

%载频

fo=fs/20;

%码率

L=300;

%信号样本

t = (0:L-1)*1/fs;

xn=cos(8*pi*fc*t);

y = xn + randn(size(t));

figure(1)

plot(t,y) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%进行能量检测NFFT = 2^nextpow2(L);

Y = fft(y,NFFT)/L;

f = fs/2*linspace(0,1,NFFT/2);

figure(2)

plot(f,2*abs(Y(1:NFFT/2)))

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%生成bpsk调制信号fs=100;

%采样频率

fc=40;

%载频

fo=fs/20;

%码率

a=2*fo;

%循环周期

T=40;

%信号样本

x=rand(40,1);

y=round(x);

a1=y*2-1;

z=ones(200,1);

for i=1:40

for j=1:5

an(5*i-5+j)=z(5*i-5+j)*a1(i);

end

t=1/fs:1/fs:2;

xn=an.*cos(4*pi*fc*t);

figure(1)

plot(t,xn);

axis([0 2 -1.5 1.5]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%进行周期图检测window=kaiser(length(xn),0.5);

%矩形窗

nfft=512;

%512fft

[Pxx,f]=periodogram(xn,window,nfft,fs);

%周期图检测

%%%%%%%%%(生成一段加噪信号）

as=rand(20,1);

b=as*2;

c=floor(b);

d=reshape(c,10,2);

e=bi2de(d,'left-msb');

h=modem.qammod(4);

g=modulate(h,e);

scatterplot(g);

%生成一段随机4QAM信号

N=length(g);

x=real(g);

y=imag(g);

fc=40;

%载波频率，有线电视信号

t0=0.000001;

ts=1e-010;

t1=0:ts:t0;

n1=0.1*randn(1,N);

for i=1:N

s1=x(i)*cos(2*pi*fc*t1)-y(i)*sin(2*pi*fc*t1); end

%生成4QAM调制后信号

%加入噪声nt

t2=0.000001:1e-010:0.000002;

M1=length(t2);

nt=0.1*randn(1,M1);

s2=nt.*cos(2*pi*fc*t2);

%add nt

t=0:ts:0.000002;

s=zeros(1,20001);

for i=1:10001

s(i)=s1(i);

end

for j=1:10000

s(j+10000)=s2(j);

end

%加入间隔无信号时间

%生成调制后信号%%%%%%%%%%%%%%%%%%%%%%（1）取有用信号时间段，计算函数值cl=xcorr(s1);

L=length(cl);

i2=1:2:L;

c2=cl(i2);

%自相关函数

a1=fc;

%循环周期

%构造积分对象pp

p=exp(-2*pi*t1*j);

M=length(p);

pp=c2.*p;

%进行无穷大积分

ts=1e-010;

cxx1=pp*ts;

u1=sum(cxx1(1:10),'double');

%加入统计频率周期a2

N1=length(c2);

N2=length(t);

f1=60e6;

%采样时间

Te=0.03;

i3=1:N1;

a2=i3./(N1*Te);

tj1=exp(a2);

cx1=tj1*u1*1e10;

U1=sum(cx1(1:10));

%%%%%%%%%%%%%%%%%%%%%%（2）取噪声时间段，计算函数值

cl2=xcorr(s2);

L=length(cl2);

i3=1:2:L;

c22=cl2(i3);

p2=exp(-2*pi*t2*j);

M2=length(p);

pp2=c22.*p2;

%进行无穷大积分

ts=1e-010;

cxx2=pp2*ts;

u2=sum(cxx2(1:10),'double');

%加入统计频率周期a2

i3=1:N1;

a2=i3./(N1*Te);

tj2=exp(a2);

cx2=tj2*u2*1e10;

U2=sum(cx2(1:10));

%检测完毕

%%%%%%%%%%%%%%%%%%% 输出检测示意图;

tz=0:ts:0.000002;

sz=zeros(1,20001);

for i=1:10001

sz(i)=U1;

end

for j=1:10000

sz(j+10000)=U2;

end

%%%%%%%%（检测完成，输出图谱）

figure(2)

plot(f,-10*log10(Pxx));

window=kaiser(length(xn),0.5);

nfft=1024; %%%%%%%%%%%%%%%%%%%%%%%%%%%%1024fft

[Pxx,f]=periodogram(xn,window,nfft,fs);

%周期图检测

figure(3)

plot(f,-10*log10(Pxx)); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%对不同类型窗函数进行检测性能分析window1=boxcar(length(xn));

%矩形窗

window2=hamming(length(xn));

%海明窗

window3=blackman(length(xn));

%blackman窗

[Pxx1,f]=periodogram(xn,window1,nfft,fs); [Pxx2,f]=periodogram(xn,window2,nfft,fs); [Pxx3,f]=periodogram(xn,window3,nfft,fs); figure(4)

subplot(221)

plot(f,-10*log10(Pxx));

subplot(222)

plot(f,-10*log10(Pxx1));

subplot(223)

plot(f,-10*log10(Pxx2));

subplot(224)

plot(f,-10*log10(Pxx3));

必修二算法的概念、程序框图(一)

内部资料，请勿外传 1 第三讲算法的概念、程序框图(一) 【考纲要求】： ①了解算法的含义、了解算法的思想． ②理解程序框图的三种基本逻辑结构：顺序、条件分支、循环．一、算法的概念 1．用加减消元法解二元一次方程组2121x y x y ?=-??í?+=?? 的具体步骤是什么？ 2．参照上述思路，一般地，解方程组 1112 22a x b y c a x b y c ì+=??í?+=?? 1221(0)a b a b -≠的基本步骤是什么？ 3．根据上述分析，用加减消元法解二元一次方程组，可以分为五个步骤进行这五个步骤就构成了解二元一次方程组的一个“算法”.我们再根据这一算法编制计算机程序，就可以让计算机来解二元一次方程组.那么解二元一次方程组的算法包括哪些内容？ 4．一般地，算法是由按照一定规则解决某一类问题的基本步骤组成的，你认为这些步骤的个数是有限的还是无限的？每个步骤是否有明确的计算任务？ 5．有人对哥德巴赫猜想“任何一个大于4的偶数都能写成两个质数之和”，设计了如下操作步骤：第一步，检验6=3+3，第二步，检验8=3+5，第三步，检验10=5+5， …… 利用计算机无穷地检验下去！请问：这是一个算法吗？ 6．根据上述分析，归纳出算法的概念：在数学中，按照一定规则解决某一类问题的明确和有限的步骤称为算法. 二、算法的步骤设计不同类型的问题有不同内容的算法，我们以判断一个整数是否为质数为例，一起来探讨算法的步骤设计. 1．如果让计算机判断7是否为质数，如何设计算法步骤？ 2．如果让计算机判断35是否为质数，如何设计算法步骤？ 3．整数89是否为质数？如果让计算机判断89是否为质数，按照上述算法需要设计多少个步骤？ 4．用2～88逐一去除89求余数，需要87个步骤，这些步骤基本是重复操作，我们可以按下面的思路改进这个算法，减少算法的步骤. （1）用i 表示2～88中的任意一个整数，并从2开始取数；（2）用i 除89，得到余数r. 若r=0，则89不是质数；若r≠0，将i 用i+1替代，再执行同样的操作；（3）这个操作一直进行到i 取88为止. 你能按照这个思路，设计一个“判断89是否为质数”的算法步骤吗？ 5．一般地，判断一个大于2的整数是否为质数的算法步骤如何设计？ ① ② ① ②

《算法的含义与流程图》测试1

《算法的含义与流程图》测试1 1.下面的结论正确的是（） A．一个程序的算法步骤是可逆的B、一个算法能够无止境地运算下去的 C、完成一件情况的算法有且只有一种 D、设运算法要本着简单方便的原则 2、早上从起床到出门需要洗脸刷牙(5 min)、刷水壶(2 min)、烧水(8 min)、泡面 (3 min)、吃饭(10 min)、听广播(8 min)几个步骤、从下列选项中选最好的一种算法 ( ) A、S1 洗脸刷牙、S2刷水壶、S3 烧水、S4 泡面、S5 吃饭、S6 听广播 B、S1刷水壶、S2烧水同时洗脸刷牙、S3泡面、S4吃饭、S5 听广播 C、S1刷水壶、S2烧水同时洗脸刷牙、S3泡面、S4吃饭同时听广播 D、S1吃饭同时听广播、S2泡面、S3烧水同时洗脸刷牙、S4刷水壶 3、闻名数学家华罗庚“烧水泡茶的两个算法、算法一：算法二：这两个算法的区不在哪里？哪个算法更高效？什么缘故？ 4、写出求 1+2+3+4+5+6……+100 的一个算法。可运用公式 1+2+3+……+ n= 2)1 ( n n 直截了当运算、 5、已知一个学生的语文成绩为89，数学成绩为96，外语成绩为99。求他的总分和平均成绩的一个算法为：第二步①；第三步② 6、“鸡兔同笼“是我国隋朝时期的数学著作《孙子算经》中的一个有味而具有深远阻碍的题目： “今有雉兔同笼，上有三十五头，下有九十四足，咨询雉兔各几何。用方程组的思想不难解决这一咨询题，请你设计一个这类咨询题的通用算法。 7、已知直角坐标系的两点A(－1，0)，B(3，2)，写出直线AB的方程的一个算法。 8.写出交换两个大小相同的杯子中的液体(A 水、 B 酒) 的两个算法。

2012年高三数学一轮复习资料第十四章算法初步第1讲算法的概念与程序框图

- 1 - 第1讲算法的概念与程序框图 ★知识梳理★ 1.算法：可以用计算机来解决的某一类问题的程序或步骤． 2.算法中的程序和步骤必须是明确和有效的，而且能够在有限步之内完成. 3.算法具有概括性（能解决一类问题），确切性（每一步操作的内容和顺序必须是明确的），有穷性（必须在有限步内结束并返回一个结果），不唯一性（一个问题可以有多个算法，算法有优劣之分），普遍性（很多具体的问题，都可以设计合理的算法去解决）. 4.程序框图又称流程图，是一种用规定的图形，指向线及文字说明来准确地、直观地表示算法的图形； 5.算法的基本逻辑结构（顺序结构、条件结构和循环结构） ①顺序结构表示语句和语句之间，框与框之间是按顺序进行的； ②条件结构是需要先根据条件作出判断，再决定执行哪一种操作的结构； ③循环结构是需要反复执行某一处理步骤的结构，分为当型（WHILE 型）和直到型（UNTIL 型），当型（WHILE 型）循环是指在每次执行循环体前对控制循环条件进行判断，当条件满足时执行循环体，不满足时停止，直到型（UNTIL 型）循环是先执行一次循环体，然后对控制循环条件进行判断，当条件不满足时执行循环体，满足则停止. ★重难点突破★ 1.重点：理解程序框图的三种基本逻辑结构，掌握三种逻辑结构在程序框图中的体现和特点． 2.难点：绘制简单实际问题的流程图，正确理解各种算法语句的实际意义． 3.重难点：设计算法时要综合考虑问题中可能涉及的各种情况：必须能解决一类问题，并且能重复使用；算法过程要一步一步执行，每一步执行的操作，必须确切，不能含糊不清，而且在有限步后得出结果．条件结构主要用在一些需要依据条件进行判断的算法中，如分段函数的求值、参数的讨论等．循环结构主要用在一些有规律的重复计算的算法中，如累加求和、累乘求积等． ★热点考点题型探析★ 考点一算法与程序框图题型1 对算法阅读能力的考查【例1】一个算法如下：第一步：计算2 44ac b m a -= ；第二步：若0>a ，输出最小值m ；第三步：若0