哪位大神有leach-l算法matlab仿真程序

实现这些功能代码可要编不少, 先把以前编过的给你一段试试看

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889
% MSK数字带通调制仿真 n=10000; %总采样点数count=100; N=n/count; %N=100 r=rand(1,N-1); % 产生一个行向量，共9999个 for i=1:N-1, if (r(i)>=0.5) r(i)=1; else r(i)=-1; end; end; %将这9999个数整数化为 1，-1 for i=1:n-count % 取9900个 R(i)=r(((i-1)-mod((i-1),count))/count+1); end; % a(1)=1; for i=2:N if r(i-1)==1 a(i)=-a(i-1); else a(i)=a(i-1); end; end; for i=1:n A(i)=a(((i-1)-mod((i-1),count))/count+1); end; for i=1:N if(mod(i,2)) l(i)=a(i); l(i+1)=a(i); else q(i)=a(i); q(i-1)=a(i); end; end; for i=1:n I(i)=l(((i-1)-mod((i-1),count))/count+1); end; for i=1:n Q(i)=q(((i-1)-mod((i-1),count))/count+1); end; ts=0.0001; Ts=ts*count; fs=1/ts; fc=5/(2*Ts); t0=Ts*N-ts; f=1/(4*Ts); df=0.3; tI=[-Ts:ts:t0-Ts]; tQ=[0:ts:t0]; tQ_R=[0:ts:t0-Ts]; c=2*pi*f; u=I.*cos(c*tI); v=Q.*sin(c*tQ); U=u.*cos(2*pi*fc*tI); V=v.*sin(2*pi*fc*tQ); for i=count+1:n W(i)=U(i)+V(i-count); end; figure subplot(4,1,1); plot(tQ_R,R,'k');axis([-Ts,0.1,-2,2]); title('数字信号'); grid on subplot(4,1,2); plot(tI,A,'k');axis([-Ts,0.1,-2,2]); title('差分编码后的信号'); grid on subplot(4,1,3); plot(tI,I,'k');axis([-Ts,0.1,-2,2]); title('同相调制信号I'); grid on subplot(4,1,4); plot(tQ,Q,'k');axis([-Ts,0.1,-2,2]); title('正交调制信号Q'); grid on figure subplot(3,1,1); plot(tI,U,'k');axis([-Ts,0.1,-2,2]); title('Icos(wct)cos(wt)');grid on subplot(3,1,2); plot(tQ,V,'k');axis([-Ts,0.1,-2,2]); title('Qsin(wct)sin(wt)');grid on subplot(3,1,3); plot(tI,W,'k');axis([-Ts,0.1,-2,2]); title('MSK信号');grid on