noise.m

来自「高频地波雷达完整仿真matlab程序 网上找到的」· M 代码 · 共 67 行

% ***********************************************************************
% *                                                                     *
% *   HFRSIM (RELAESE 2.1   31/07/2006  ).                              *
% *   COPYRIGHT (C)   HIT EIE  2006.                                    *
% *                                                                     *
% ***********************************************************************
% 
% ================================================================
% =   MODIFIED FOR MATLAB Version 7.0.4.365 (R14) Service Pack 2 = 
% =   FOR WINDOWS PERSONAL COMPUTERS                             =
% =   July 2006                                                  =
% =   Mr. Zhou Yang                                              =
% =   Harbin Institute of Technology                             =
% =   mailto:[e-mail]                               =
% ================================================================
function [grFreq,grATMO,grGAL,grMANMADE,grOVERALL]=noise(FMIN,FSTEP,FMAX,season,LAT,LON,BW,TIME)

% Input Variables
% season = l=WINTER, 2=SPRING, 3=SUMMER, 4=AUTUMN, 5=ALL SEASONS
% LAT = LOCATION LATITUDE 
% LON = LOCATION LONGITUDE
% BW = BANDWIDTH
% TIME = SPECIFIC LOCAL MEAN TIME

gr=[];  % 定义一个空矩阵。
sign=floor((FMAX-FMIN)/FSTEP)+1;
flag=0;
progressbar; % 调用自制的进度条

for FREQ=FMIN:FSTEP:FMAX,
    
    flag=flag+1;
%     pause(0.1);
    progressbar(flag/sign); % 更新进度条
    
    % Write the input variables to file(NOISBW.INP).
    fid = fopen('NOISBW.INP','w');
    fprintf(fid,'%f\n%f,%f,''TEST''\n%f,1\n%f\n%f\n1\n0\n',...
        season,LAT,LON,BW,FREQ,TIME);
    fclose(fid);
    
    % Execute NOIS1 command and return result.
    dos NOIS1<NOISBW.INP;
    
    % Read the output data from file(NOISBW.OUT).
    fid = fopen('NOISE.LST', 'r');
    tmp = fscanf(fid, '%g %g %g %g %g %g %g %g %g %g', [10 inf]);   % It has ten rows now.
    tmp = tmp';
    fclose(fid);
    
    % 合并矩阵。
    gr=[gr;tmp];
    
end

% 单独输出矩阵的每一列。
grFreq=gr;
grFreq(:,2)=[];grFreq(:,2)=[];grFreq(:,2)=[];grFreq(:,2)=[];grFreq(:,2)=[];grFreq(:,2)=[];grFreq(:,2)=[];grFreq(:,2)=[];grFreq(:,2)=[];
grATMO=gr;
grATMO(:,1)=[];grATMO(:,2)=[];grATMO(:,2)=[];grATMO(:,2)=[];grATMO(:,2)=[];grATMO(:,2)=[];grATMO(:,2)=[];grATMO(:,2)=[];grATMO(:,2)=[];
grGAL=gr;
grGAL(:,1)=[];grGAL(:,1)=[];grGAL(:,2)=[];grGAL(:,2)=[];grGAL(:,2)=[];grGAL(:,2)=[];grGAL(:,2)=[];grGAL(:,2)=[];grGAL(:,2)=[];
grMANMADE=gr;
grMANMADE(:,1)=[];grMANMADE(:,1)=[];grMANMADE(:,1)=[];grMANMADE(:,2)=[];grMANMADE(:,2)=[];grMANMADE(:,2)=[];grMANMADE(:,2)=[];grMANMADE(:,2)=[];grMANMADE(:,2)=[];
grOVERALL=gr;
grOVERALL(:,1)=[];grOVERALL(:,1)=[];grOVERALL(:,1)=[];grOVERALL(:,1)=[];grOVERALL(:,2)=[];grOVERALL(:,2)=[];grOVERALL(:,2)=[];grOVERALL(:,2)=[];grOVERALL(:,2)=[];