matlab实现OFDM信道估计LS算法

function Q6_1_ErrorofChannelEstimation(block) % Q5_1_Spread 用于扩频 setup(block); %endfunction function setup(block) %% Register number of input and output ports block.NumInputPorts = 2; block.NumOutputPorts = 1; %定义两个输入变量，第一个为PN序列周期 %第二个为每个时隙的数据个数 block.NumDialogPrms = 3; %数据载波数 Nd=block.DialogPrm(1).Data; %补零波数 Nz=block.DialogPrm(2).Data; %导频载波数 Np=block.DialogPrm(3).Data; %% Setup functional port properties to dynamically %% inherited. block.SetPreCompInpPortInfoToDynamic; block.SetPreCompOutPortInfoToDynamic; % 输入数据符号 block.InputPort(1).Dimensions = [Nd+Nz,1]; block.InputPort(1).DirectFeedthrough = true; block.InputPort(1).SamplingMode = 'Frame' block.InputPort(2).Dimensions = [Nd+Nz+Np,1]; block.InputPort(2).DirectFeedthrough = true; block.InputPort(2).SamplingMode = 'Frame' block.OutputPort(1).Dimensions = [1,1]; block.OutputPort(1).SamplingMode = 'Frame' %% Set block sample time to inherited block.SampleTimes = [-1 0]; %% Register methods block.RegBlockMethod('Outputs', @Output); block.RegBlockMethod('PostPropagationSetup', @DoPostPropSetup); block.RegBlockMethod('InitializeConditions', @InitConditions); %endfunction function Output(block) He=block.InputPort(1).Data; H=block.InputPort(2).Data; %数据载波数 Nd=block.DialogPrm(1).Data; %补零波数 Nz=block.DialogPrm(2).Data; %导频载波数 Np=block.DialogPrm(3).Data; Nc=Nd+Nz+Np; %提取数据符号（包括零） L=Nc./Np; H1=[]; for k=1:Np index=(k-1).*L; H1=[H1;H(index+1:index+L-1)]; end Error=sum((abs(He)-abs(H1))).^2; block.Dwork(1).Data=block.Dwork(1).Data+Error; block.Dwork(2).Data=block.Dwork(2).Data+Nd+Nz; block.OutputPort(1).Data = block.Dwork(1).Data/block.Dwork(2).Data; %endfunction function DoPostPropSetup(block) %% Setup Dwork %定义变量个数 block.NumDworks = 2; block.Dwork(1).Name = 'x1'; block.Dwork(1).Dimensions = 1; block.Dwork(1).DatatypeID = 0; block.Dwork(1).Complexity = 'Real'; block.Dwork(1).UsedAsDiscState = true; block.Dwork(2).Name = 'x2'; block.Dwork(2).Dimensions = 1; block.Dwork(2).DatatypeID = 0; block.Dwork(2).Complexity = 'Real'; block.Dwork(2).UsedAsDiscState = true; %endfunction function InitConditions(block) %% Initialize Dwork block.Dwork(1).Data = 0; block.Dwork(2).Data = 0;