function [ ] = wentai( para )
%wentai 判断是否算稳
% 画径向中心处的电子密度和电子温度
path0 = pwd;
listing = dir([path0]); %获取路径下的所有文件(夹)
length_listing = length(listing); %获取文件(夹)数目
count_num = 0;
for k = 1:1:length_listing %对所有文件(夹)进行循环
if ~isempty(strfind(listing(k).name,para)) %对文件(夹)名字中含有MHz的进行操作
count_num = count_num +1 ;
para1(count_num) = str2num( listing(k).name(1:end - length(para)));
%获取放电参数值
end
end
para1 = sort(para1);
for k = 1:1:count_num
path1 = [path0,'\',num2str(para1(k)),'.0',para];
str_loadrhoe =strcat(path1,'\','rhoem','.dat');
rhoe0 = load(str_loadrhoe);
rhoe1 = rhoe0(:,1); %改
str_loadte =strcat(path1,'\','Te','.dat'); %
Te0 = load(str_loadte);
Te1 = Te0(:,1); %改
figure(k)
subplot(2,1)
plot(rhoe1);title('径向中心处的N_{e}');
subplot(2,2)
plot(Te1);title('径向中心处的T_{e}')
end
end
判断程序是否算稳的小程序
最新推荐文章于 2023-09-05 00:10:05 发布