0

この問題は些細なことのように思えますが、解決しようとすると頭を悩ませています。64 QAM コンスタレーションに定数モジュラス手法を使用して分数間隔イコライザーを適用しようとしています。プログラムは QPSK または 4 QAM で動作しますが、64QAM に適用するとエラーがスローされます。

Error using  / 
Matrix dimensions must agree.

Error in Working_FSE_CMA_64QAM (line 68)
sb1=sb/(fh(temp));  % scale the output

通信ツールボックスを持っていないので、前の質問の回答を使用して 64QAM シンボルを生成しました 16 QAM 信号を生成します

誰かがコードを機能させるのを手伝ってくれませんか? ありがとうございました。

% Blind channel estimation/equalization
% adpative CMA method in Fractional space

T=1000;    % total number of data
dB=25;     % SNR in dB value

%%%%%%%%% Simulate the Received noisy Signal  %%%%%%%%%%%
N=5; % smoothing length N+1
Lh=5;  % channel length = Lh+1
Ap=4;  % number of subchannels or receive antennas

h=randn(Ap,Lh+1)+sqrt(-1)*randn(Ap,Lh+1);   % channel (complex)
for i=1:Ap, h(i,:)=h(i,:)/norm(h(i,:));    end        % normalize
s = (randi(8,1,T)*2-5)+j*(randi(8,1,T)*2-5);  %64 QAM
%s=round(rand(1,T))*2-1;  % QPSK or 4 QAM symbol sequence
%s=s+sqrt(-1)*(round(rand(1,T))*2-1);

% generate received noisy signal
x=zeros(Ap,T);    % matrix to store samples from Ap antennas
SNR=zeros(1,Ap);
for i=1:Ap
    x(i,:)=filter(h(i,:),1,s);
    vn=randn(1,T)+sqrt(-1)*randn(1,T);   % AWGN noise (complex)
    vn=vn/norm(vn)*10^(-dB/20)*norm(x(i,:));  % adjust noise power
    SNR(i)=20*log10(norm(x(i,:))/norm(vn));   % Check SNR of the received samples
    x(i,:)=x(i,:)+vn;                        % received signal
end
SNR=SNR    % display and check SNR

%%%%%%%%%%%%% adaptive equalizer estimation via CMA
Lp=T-N;   %% remove several first samples to avoid 0 or negative subscript
X=zeros((N+1)*Ap,Lp);  % sample vectors (each column is a sample vector)
for i=1:Lp
    for j=1:Ap
        X((j-1)*(N+1)+1:j*(N+1),i)=x(j, i+N:-1:i).';
    end
end

e=zeros(1,Lp);  % used to save instant error
f=zeros((N+1)*Ap,1); f(N*Ap/2)=1;    % initial condition
%R2=2;                  % constant modulas of QPSK symbols
R2 = 1.380953; %For 64 QAM http://www.google.com/patents/US7433400
mu=0.001;      % parameter to adjust convergence and steady error
for i=1:Lp
   e(i)=abs(f'*X(:,i))^2-R2;                  % instant error
   f=f-mu*2*e(i)*X(:,i)*X(:,i)'*f;     % update equalizer 
   f(N*Ap/2)=1;
%   i_e=[i/10000 abs(e(i))]             % output information 
end   

%sb=f'*X;   % estimate symbols (perform equalization)
sb = filter(f, 1, X);
% calculate SER
H=zeros((N+1)*Ap,N+Lh+1);  temp=0;
for j=1:Ap
    for i=1:N+1, temp=temp+1; H(temp,i:i+Lh)=h(j,:); end  % channel matrix
end

fh=f'*H;    % composite channel+equalizer response should be delta-like 
temp=find(abs(fh)==max(abs(fh)));   % find the max of the composite response

sb1=sb/(fh(temp));  % scale the output
sb1=sign(real(sb1))+sqrt(-1)*sign(imag(sb1));  % perform symbol detection
start=N+1-temp;  % general expression for the beginning matching point
sb2=sb1(10:length(sb1))-s(start+10:start+length(sb1));  % find error symbols
SER=length(find(sb2~=0))/length(sb2)   % calculate SER

if 1
    subplot(221), 
    plot(s,'o');   % show the pattern of transmitted symbols
    grid,title('Transmitted symbols');  xlabel('Real'),ylabel('Image')
    axis([-2 2 -2 2])

    subplot(222),
    plot(x,'o');  % show the pattern of received samples
    grid, title('Received samples');  xlabel('Real'), ylabel('Image')

    subplot(223),
    plot(sb,'o');   % show the pattern of the equalized symbols
    grid, title('Equalized symbols'), xlabel('Real'), ylabel('Image')

    subplot(224),
    plot(abs(e));   % show the convergence
    grid, title('Convergence'), xlabel('n'), ylabel('Error e(n)')
end
4

1 に答える 1

0

主な問題は、反復アルゴリズムが収束しない (具体的には発散しない) ことです。結果として、値はfcontainsになりますNaN。結果として

temp=find(abs(fh)==max(abs(fh)))

は空のベクトルであり、ライン上の予想されるスカラー サイズと一致しません。

sb1=sb/(fh(temp));  % scale the output

muこの問題を解決するには、使用するとの値がR2、単位分散を持つコンスタレーションに基づいていることに注意してください。このような 64-QAM コンスタレーションは、次の方法で生成できます。

s = ((randi(8,1,T)*2-9)+j*(randi(8,1,T)*2-9))/sqrt(42);  %64 QAM

j=sqrt(-1)また、インデックス変数として使用するために、数行で複素定数を再定義しています。jas index の使用は避けてください。また、sript の先頭ですべての変数をクリアすると ( を使用clear all)、一貫したクリーンな状態で実行を開始できます。

于 2015-11-14T15:49:23.887 に答える