readfile.m

污染物连续衰变浓度变化 基于matlab的编程

function [v,D,num,c0,k,t,xmin,xmax,interval]=Readfile

[filename,pathname] = uigetfile('input.txt','Inputfile');
fid = fopen(strcat(pathname,filename),'rt');
F=load(strcat(pathname,filename));
fclose(fid);
n=size(F,1);
R=F(4);
v=F(1)/R;
D=F(2)/R;
t=F(3);
xmin=F(5);
xmax=F(6);
interval=F(7);
num=F(8);
c0=zeros(num,1);
k=zeros(num,1);
for i=1:num
    c0(i)=F(i+8);
end
for i=1:num
    k(i)=F(num+8+i)/R;
end