sin_maker_ok.m

mimo2x2天线选择系统的全系统matlab程序,先前的是dsp程序.

clear,clf,close all
% 
fs = 48000;
f1 = 5000;
f2 = 11000;
N = 8*1024;
t = (0:(N/2-1))/fs;
x1 = [sin(2*pi*f1*t), sin(2*pi*f1*t)];
%x2 = [sin(2*pi*f1*t), sin(2*pi*f2*t)];

% f=linspace(-1/2,1/2,length(x1));
% X1=fftshift(fft(x1));
% subplot(211),plot(t,x1), grid
% 
% f=linspace(-1/2,1/2,length(x2));
% X2=fftshift(fft(x2));
% subplot(222),plot(f,abs(X2)), grid

Xout = [x1; x1];

Xout = Xout / (max(max(abs(Xout))));
Xout = Xout * (2^15 - 1);
% xout1=Xout(1,:);
% f=linspace(-1/2,1/2,length(xout1));
% Xout1=fftshift(fft(xout1));
% subplot(212),plot(t,xout1), grid


%Xout =0* randn(2,N);
% figure(2)
% subplot(221)
% plot(Xout(1,:))
% title('Audio Out, left')
% subplot(222)
% plot(Xout(2,:))
% title('Audio Out, right')

[Xin,no_samples_delivered,buffer_size]=DSK6713_rxtx(N,2,2,Xout,fs);

figure(2)
xin1=Xin(1,:);
xin2=Xin(2,:);
% 
% f=linspace(-1/2,1/2,length(xin1));
% Xin1=fftshift(fft(xin1));
subplot(211),plot(xin1), grid
% 
% f=linspace(-1/2,1/2,length(xin2));
% Xin2=fftshift(fft(xin2));
% subplot(212),plot(f,abs(Xin2)), grid

% % plot(Xin')
% figure(1)
% subplot(223)
% plot(Xin(1,:))
% title('Audio In, left')
% figure(1)
% subplot(224)
% specgram(Xin(2,:),256,fs)
% title('Audio In, right')
title('y(n)')
xlabel('v=fo/fs (fs=96000)')
ylabel('power')