Kraken波导不变量计算matlab程序

%% ***********Sound Pressure Calculation using KRAKEN*************
%————————————Designed by Okazaki_Jian———————————————
%————————————————07/03/2015——————————————————
%——————————————All Rights Reserved———————————————
%————This programm needs the installation of ‘acoustic toolbox‘————
%% Setting Environment Variables


%% Resetting Workspace
clc;
clear all
close all
lofar_env = ‘KRAKEN_Munk‘;
dos（[‘del ‘ lofar_env ‘.mod‘]）;
dos（[‘del ‘ lofar_env ‘.shd‘]）;
dos（[‘del ‘ lofar_env ‘.prt‘]）;
%% TL Calculation using KRAKEN
rdt = 200; %depth of reciever
    f = 100;
    cd（ ‘D:\deep_sea_waveguide\Main_Programme‘ ）;
    %attach ‘.env‘ to change source frequency
    fir = fopen（ [lofar_env ‘.env‘]  ‘r+‘ ）;
    fseek（fir 15  -1 ）;
    fprintf（ fir ‘%4.6f‘ f ）;
    fclose（‘all‘）;
    runkraken = which（ ‘krakenc.exe‘ ）;
    eval（ [ ‘! “‘ runkraken ‘“ ‘ lofar_env ] ）;
    % pressure calculation using field
    plotssp（lofar_env）;
    sofarmidflp = [lofa