Matlab画图所用到
循环命名label
Ylabel(['y',num2str(10-i)])
读取TXT文件
datafilename = 'winequality-red.txt'; D = load(datafilename);
读取data文件
datafilename = 'housing.data'; D = load(datafilename);
读取xls文件
[num] = xlsread('Folds5x2_pp.xlsx');
D = num;
读取文件
% %------------------数据 6 DA delta Ailerons
% datafilename = 'delta_ailerons.data';
% D = load(datafilename);
% %------------------数据 7 Red Wine quality
% datafilename = 'winequality-red.txt';
% D = load(datafilename);
%
%-----------------数据 8 Combined Cycle Power Plant Data Set
[num] = xlsread('Folds5x2_pp.xlsx');
D = num;
循环读取文件
Num_c = 16; %需要读取的文件个数
%循环读取数据，
for i=1:1:Num_c
n = 8 + (i-1)*4;
fileName = "experiments_for_different_alpha_beta_J_results_c"+n+".txt";
All_data{i} = load(fileName);
end
循环保存TXT文件
%保存文件
%filename=sprintf('syth_h%d_c%d.txt',h,c)
%fid = fopen(filename, 'w+');
fprintf(fid,'%.4f %.4f %.4f \r\n', a1,Q_tr,Q_te) ; %写入TXT文件中
画竖直的线 stem
xx = [21 12 30 6 -10 1]; x_x = -3:2; figure stem(x_x,xx)
设置画图句柄
set (gcf,'Position',[400,240,400,240], 'color','w'); %设置图片属性，位置，背景图颜色，默认为灰色，设置为白色 set(gca,'linewidth',2,'fontsize',30,'fontname','Times'); %设置坐标轴属性，线宽4，刻度字号30，刻度字体罗马，
xlabel('\itx\rm_{1}','FontName','Times New Roman','FontSize',15);
画三维柱状图
datafilename = 'experiments_for_different_alpha_beta_J_results_c70.txt';
full_data = load(datafilename);
J_plot = reshape(full_data(:,3),8,8);
figure
bar3(J_plot,0.25,'detached')
%s设置坐标轴
xlabel('\it \alpha','FontName','Times New Roman');
ylabel('\it \beta','FontName','Times New Roman');
zlabel('\it J','FontName','Times New Roman');
set (gcf, 'color','w')
%set (gcf,'Position',[400,100,299,299], 'color','w')
set(gca,'FontName','Times New Roman');
set(gca,'YTickLabel',{'1';'3';'5';'7';'9';'11';'13';'15'})
set(gca,'XTickLabel',{'begin';'3';'5';'7';'9';'11';'13';'15'})
%title('Width = 0.25')
其中用到的柱状改色，链接参考上一篇博文
http://blog.sina.com.cn/s/blog_739bfa860102uwa5.html
自定义柱状图
将表格插入图片中
f = figure(1);
set(f,'Position',[500 500 300 150]);
dat = {' a', 1, ' units';...
' b', 2, ' units';...
' c', 3, ' units';...
' d', 4, ' units';...
' e', 5, ' units';...
' f', 6, ' units';};
columnname = {'Parameter', 'Value', 'Units'};
columnformat = {'char', 'numeric', 'char'};
t = uitable('Units','normalized','Position',...
[0.05 0.05 0.755 0.87], 'Data', dat,...
'ColumnName', columnname,...
'ColumnFormat', columnformat,...
'RowName',[]);
图片中，局部放大的图片
clear, clc
% 加载图像到内存，如读取matlab自带的图像，并将其显示出来
I = imread('pout.tif');
imshow(I, []);
f = gcf;
% 调用magnify，magnify 接收的参数是 figure
magnify(f); % magnify 函数的处理是一种交互式的
这一部分使用的是自带的画中画
t = 0:0.1:2*pi;
sig = sin(t);
plot(t,sig,'b.-');
t2 = 1.2:0.05:1.99;
sig2 = sin(t2);
ax = axes('Position',[0.5 0.8 0.2 0.075]);
plot(ax,t2,sig2,'b.-');
axis([1.2 1.99 0.9 1.01]);
x = [0.33 0.5]
y = [0.91 0.83]
annotation('textarrow',x,y)
读取文件，保存图片或者TXT文件
%文件的读取
datafilename = 'airfoil_c_vs_avg_rho.txt';
original_data = load(datafilename);
%自动保存图片
figure
set (gcf,'Position',[400,400,400,400], 'color','w');
plot(c_list,Q_list,'k--',c_list,Q_TS_list,'k-','LineWidth',2)
xlabel('\it c','FontName','Times New Roman');
ylabel('\it Q','FontName','Times New Roman');
legend('\it Q_{\rm proposed model}','\it Q_{\rm TS model}');
xlim([c_value_min c_value_max])
%自动保存图片
%
% fig_name=sprintf('experiments_for_different_alpha_beta_J_results_c%d',fig_num);
fig_name = sprintf('airfoil_Qlist_vs_clist')
print(gcf,'-dpng',[fig_name,'.pdf'])
print(gcf,'-dpng',[fig_name,'.fig'])
save('airfoil_c_vs_avg_eps.txt','eps_record','-ascii');
3D图片的设置
figure1 = figure('Color',[1 1 1]);
axes1 = axes('Parent',figure1,'PlotBoxAspectRatio',[1 0.5 0.5],'FontName','Times New Roman');
view(axes1,[22 32]);
%axis([0 1 0 1 0 1]);
axis square;
grid on;
for i=1:1:c_value
x0 = input_granules(i,1);
y0 = input_granules(i,2);
radius = input_granules(i,end);
z0 = output_granules(i,1);
height = output_granules(i,end) *2;
my_plot_cylinder(x0, y0, z0, radius, height);
%my_plot_sphere(x0, y0, z0, radius, height);
end
%同时，将原数据点画到其中
for i=1:1:total_rows
hold on; plot3(training_data(i,1), training_data(i,2), training_data(i,3),'k.','Markersize',6);
end
xlabel('\it x_{\rm 1}','FontName','Times New Roman');
ylabel('\it x_{\rm 2}','FontName','Times New Roman');
zlabel('\it y','FontName','Times New Roman');
%set(gcf,'Position',[400 100 299 299],'color','w');
box on;ax = gca;ax.BoxStyle = 'full';
set(gca)
%set(gca,'xtick',[-0.5 0 0.5 1]);set(gca,'ytick',[-1 -0.5 0 0.5 1.2]);set(gca,'ztick',[-0.2 -0.1 0 0.1 0.2]);
set (gcf,'Position',[400,400,400,400], 'color','w')
句柄设置
set (gcf, 'color','w')
set(gca,'FontName','Times New Roman');
set(gca,'YTickLabel',{'0';'3';'6';'9';'12';'15';'18';'21';'24'})
set(gca,'XTickLabel',{'0';'3';'6';'9';'12';'15';'18';'21';'24'})
set (gcf,'Position',[400,400,400,400], 'color','w');
画三维柱状图
bar3(J_plot,0.25,'detached')
%设置坐标轴
%xlabel('\it \alpha','FontName','Times New Roman');
%ylabel('\it \beta','FontName','Times New Roman');
%zlabel('\it J','FontName','Times New Roman');
set (gcf, 'color','w')
set(gca,'FontName','Times New Roman');
set(gca,'YTickLabel',{'0';'3';'6';'9';'12';'15';'18';'21';'24'})
set(gca,'XTickLabel',{'0';'3';'6';'9';'12';'15';'18';'21';'24'})
xlabel('\it \alpha');ylabel('\it \beta');zlabel('\it J')
%set(gca,'YTickLabel',['0';'12';'24'])
%set(gca,'Xtick',[ ]);set(gca,'ytick',[ ]);set(gca,'ztick',[0 max_J])%图画的比较多的时候，用这个命令取消坐标轴示数
%set (gcf,'Position',[800,700,299,299], 'color','w')
%自动保存图片
fig_name=sprintf('experiments_for_different_alpha_beta_J_results_c%d',fig_num);
print(gcf,'-dpng',[fig_name,'.png'])
画三维柱状图，并保存图片其中，图按照一定顺序命名
%单独画一个图
for i=14%1:1:N
fig_num = 8 + 4*(i-1);
figure
%读取数据
full_data = All_data{i};
max_J = max(full_data(:,3));
max_J = roundn(max_J,-2);%取小数点后两位有效数字的
J_plot = reshape(full_data(:,3),9,9)';%转置后，新的矩阵按照行，从新排列
bar3(J_plot,0.25,'detached')
%设置坐标轴
%xlabel('\it \alpha','FontName','Times New Roman');
%ylabel('\it \beta','FontName','Times New Roman');
%zlabel('\it J','FontName','Times New Roman');
set (gcf, 'color','w')
set(gca,'FontName','Times New Roman');
set(gca,'YTickLabel',{'0';'3';'6';'9';'12';'15';'18';'21';'24'})
set(gca,'XTickLabel',{'0';'3';'6';'9';'12';'15';'18';'21';'24'})
xlabel('\it \alpha');ylabel('\it \beta');zlabel('\it J')
%set(gca,'YTickLabel',['0';'12';'24'])
%set(gca,'Xtick',[ ]);set(gca,'ytick',[ ]);set(gca,'ztick',[0 max_J])%图画的比较多的时候，用这个命令取消坐标轴示数
%set (gcf,'Position',[800,700,299,299], 'color','w')
%自动保存图片
fig_name=sprintf('experiments_for_different_alpha_beta_J_results_c%d',fig_num);
print(gcf,'-dpng',[fig_name,'.png'])
end
MATLAB画图去除白边
addpath('D:\Matlab_my_additional_path\altmany-export_fig-a83b407')%去白边路径函数
openfig('D2_inf_grans_inputspace_beta0.fig')
%x = [-10:0.1:10];
%y = sin(x);
%plot(x,y); % 画图
export_fig( gcf , '-pdf' , '-r256 ', '-painters' , './Test.pdf'); % 测试xpdf模块，能不能正常生成pdf
export_fig( gcf , '-eps' , '-r256' , '-painters' , './Test.pdf'); % 测试ghostscript模块，能不能生成eps
MATLAB去除白边，从源程序
set(gca,'LooseInset',get(gca,'TightInset'))%只要这一个命令就可以了
figure
set (gcf,'Position',[400,280,400,280], 'color','w');
set(gca,'FontName','Times New Roman','linewidth',1.0,'fontsize',14);
x = 0:pi/100:2*pi;
y = sin(x);
plot(x,y);
xlabel('$\beta$','Interpreter','latex','fontsize',18);
ylabel('$Sp$','Interpreter','latex','fontsize',18);
set(gca,'LooseInset',get(gca,'TightInset'))%只要这一个命令就可以了
text(2,10,'$\omega f$','Interpreter','latex')
h=legend('$ex\eta$')
set(h,'Interpreter','latex')
% Gran control
gran_control_data = load('gran_pole_8_12.txt');
N = length(gran_control_data);
gran_rolling_data = [];
for k=1:1:N-window_size
sum_data = sum(gran_control_data(k:k+window_size-1));
gran_rolling_data = [gran_rolling_data;sum_data/window_size];
end
x_var = 1:1:length(gran_rolling_data);
figure
set (gcf,'Position',[400,240,400,240], 'color','w');
set(gca,'FontName','Times New Roman','linewidth',1.0,'fontsize',15);
h1 = plot(x_var,Q_rolling_data);hold on;
h2 = plot(x_var,basic_rolling_data);hold on;
h3 = plot(x_var,random_rolling_data);hold on;
h4 = plot(x_var,gran_rolling_data);hold on;
% h1 = plot(beta_cov_D2(:,1),beta_cov_D2(:,2),'ko-','LineWidth',1.5);hold on;
% h2 = plot(beta_cov_abanlone(:,1),beta_cov_abanlone(:,2),'bv-','LineWidth',1.5);hold on;
% h3 = plot(beta_cov_qsar(:,1),beta_cov_qsar(:,2),'g*-','LineWidth',1.5);hold on;
% h4 = plot(beta_cov_redwine(:,1),beta_cov_redwine(:,2),'mp-','LineWidth',1.5);hold on;
% plot(cov(:,2),sp(:,2),cov_TS0(:,2),sp_TS0(:,2),cov_TS1(:,2),sp_TS1(:,2))
%plot(r2_Q(:,1),r2_Q(:,end-1),'o-',r3_Q(:,1),r3_Q(:,end-1),'+-',r4_Q(:,1),r4_Q(:,end-1),'*-',...
%r5_Q(:,1),r5_Q(:,end-1),'x-',r6_Q(:,1),r6_Q(:,end-1),'.-',r7_Q(:,1),r7_Q(:,end-1),'x-','LineWidth',1.5)
%legend({'Q-table','Basic','Random','G-table'},'NumColumns',2)
%legend('boxoff')
%plot(cov(:,2),sp(:,2),'o-',cov_TS0(:,2),sp_TS0(:,2),'v-',cov_TS1(:,2),sp_TS1(:,2),'*-')
%legend({'synthetic data','abanlone','QSAR','red wine quality'},'NumColumns',1)
xlabel('Trials','FontName','Times New Roman','fontsize',15);
ylabel('Reward','FontName','Times New Roman','fontsize',15);
grid on
%stem(1130,170,'Marker','none');
box1 = [1130 1130 1500 1500];
boxy = [0 200 200 0];
patch(box1,boxy,[0.7 0.7 0.7],'FaceAlpha',0.3);
lgd1 = legend([h1,h2,h3,h4],{'Q-table','Basic','Random','G-table'},'NumColumns',2)
set(lgd1,'FontName','Times New Roman','FontSize',14,'FontWeight','normal','Orientation','horizon')
legend boxoff;
errorbar 和 三维重构图
%% 下面画重构图
%上一步生成的one_one_in_output_granules 对应traning_data 元数据
%这一步，利用one_one_in_output 画出 误差棒状图
figure
% x = one_one_in_output_granules(:,1);
% y=one_one_in_output_granules(:,1);
% variance = one_one_in_output_granules(:,2);
% errorbar(x,y,variance,'.')
plot(check_data(:,end),check_data(:,end),'LineWidth',1.5); %画出target值
hold on
errorbar(check_data(:,end),Y_output(:,1),Y_output(:,2),'.');
%errorbar(one_one_in_output_granules(:,1),one_one_in_output_granules(:,1),one_one_in_output_granules(:,2),'.'); %画出 granule model值
set(gca,'xtick',[0 0.5 1.0]);set(gca,'ytick',[0 0.5 1.0 ]);
set(gca,'FontName','Times New Roman','linewidth',1.0,'fontsize',15);
set (gcf,'Position',[400,400,400,400], 'color','w')
xlabel(' target output','FontName','Times New Roman','FontSize',15);
ylabel(' granule model output','FontName','Times New Roman','FontSize',15);
box off
%画三维重构图
figure
X1 =check_data(:,1);
Y1 = check_data(:,2);
Z1 = check_data(:,3);
Z0 = Y_output(:,1)-Y_output(:,2);%信息粒下面
Z2 = Y_output(:,1)+Y_output(:,2);%对应信息粒上面
[X_grid,Y_grid,Z1_grid]=griddata(X1,Y1,Z1',linspace(min(X1),max(Y1))',linspace(min(Y1),max(Y1)),'v4');
[X_grid,Y_grid,Z0_grid]=griddata(X1,Y1,Z0',linspace(min(X1),max(Y1))',linspace(min(Y1),max(Y1)),'v4');
[X_grid,Y_grid,Z2_grid]=griddata(X1,Y1,Z2',linspace(min(X1),max(Y1))',linspace(min(Y1),max(Y1)),'v4');
%figure,mesh(X_grid,Y_grid,Z_grid)%三维曲面
%colormap hsv
%hold on
mesh(X_grid,Y_grid,Z0_grid,'FaceAlpha',1)%三维曲面
%colormap gray
hold on
mesh(X_grid,Y_grid,Z2_grid,'FaceAlpha',1)%三维曲面
%colormap pink
hold on
plot3(X1,Y1,Z1,'r.','MarkerSize',6)
set(gca,'FontName','Times New Roman','linewidth',1.0,'fontsize',15);
%自定义坐标轴等区域
xlabel('\itx\rm_{1}','FontName','Times New Roman','FontSize',15);
ylabel('\itx\rm_{2}','FontName','Times New Roman','FontSize',15);
zlabel('\it y','FontName','Times New Roman','FontSize',15);
grid on
box off;ax = gca;ax.BoxStyle = 'full';
%set(gca,'xtick',[-0.5 0 0.5 1]);set(gca,'ytick',[-1 -0.5 0 0.5 1.2]);set(gca,'ztick',[-0.2 -0.1 0 0.1 0.2]);
set (gcf,'Position',[400,400,400,400], 'color','w')
% 这一部分画原始数据3D 被切割图
mesh_y = a_list(1:end-1)';%[0;0.33;0.66]
figure
D = original_full_data;
plot3(D(:,1),D(:,2),D(:,3),'.')
%自定义坐标轴等区域
xlabel('$x\rm_{1}$','Interpreter','latex','FontName','Times New Roman','FontSize',18);
ylabel('$x\rm_{2}$','Interpreter','latex','FontName','Times New Roman','FontSize',18);
zlabel('y','FontName','Times New Roman','FontSize',18);
%grid on
%box off;
%grid on
box on
ax = gca;ax.BoxStyle = 'full';
%set(gca,'xtick',[-0.5 0 0.5 1]);set(gca,'ytick',[-1 -0.5 0 0.5 1.2]);set(gca,'ztick',[-0.2 -0.1 0 0.1 0.2]);
%set (gcf,'Position',[400,400,400,400], 'color','w')
set (gcf,'Position',[400,280,400,280], 'color','w');
set(gca,'linewidth',1.0,'fontsize',14,'fontname','Times New Roman');
%设置坐标轴属性，线宽4，刻度字号30，刻度字体罗马，
set(gca,'LooseInset',get(gca,'TightInset'))
for i=1:1:size(mesh_y)
hold on
temp_mesh_y = mesh_y(i);
if temp_mesh_y == 0
continue
else
hold on
[X, Y] = meshgrid(0:0.01:1, 0:0.01:1);
Z = ones(size(X))*temp_mesh_y;
surf(X,Y,Z,'FaceAlpha',0.5,'EdgeColor','none','FaceColor','red') %FaceAlpha定义透明度
end
end
hold off
view([1,1,1])
%------------------------------------------fuzzy
fuzzy_Y = load('fuzzy_standard_testing_Y_output.txt');
fuzzy_y = load('fuzzy_standard_testing_y.txt');
fuzzy_u = load('fuzzy_standard_testing_u.txt');
range_y = max(fuzzy_y)-min(fuzzy_y);
N_data = size(fuzzy_Y,1)
figure
set (gcf,'Position',[400,400,400,400], 'color','w');
for k=1:0.1*N_data
temp_Y = fuzzy_Y(k,:);
temp_y = fuzzy_y(k,:);
temp_u = fuzzy_u(k,:);
Y_min = min(temp_Y);
Y_max = max(temp_Y);
x_coord = [temp_y-y_r,temp_y-y_r,temp_y+y_r,temp_y+y_r,temp_y-y_r];
y_coord = [Y_min, Y_max, Y_max, Y_min, Y_min];
u_coord = temp_u*y_coord;
%画外框
plot(x_coord,y_coord,'k');
%画阴影，填充阴影时，去掉外框,none
patch(x_coord,u_coord,'k','FaceAlpha',0.2,'LineStyle','none');
%画原始点
plot(temp_y,temp_y,'r.','MarkerSize',12)
hold on
end
set(gca,'FontName','Times New Roman','linewidth',1.5,'fontsize',15);
set(gca,'LooseInset',get(gca,'TightInset'))
xlabel('target output','FontName','Times New Roman','fontsize',18);
ylabel('fuzzy granular output','FontName','Times New Roman','fontsize',18);
box off
addpath('C:\Users\Lenovo\Desktop\conditional_clustering\resulting_data_d2')
%parabolic_Y= load('parabolic_standard_testing_Y_output.txt');
%parabolic_y = load('parabolic_standard_testing_y.txt');
parabolic_Y = load('parabolic_standard_testing_Y_output.txt');
parabolic_y = load('parabolic_standard_testing_y.txt');
N_num = 0.2* size(parabolic_Y,1);
N_lines = 5;
Y_up_interval = [];
Y_down_interval= [];
for k=1:N_num
Y_up_interval(k,1) = (parabolic_Y(k,3) - parabolic_Y(k,2))/N_lines;
Y_down_interval(k,1) = (parabolic_Y(k,2) - parabolic_Y(k,1))/N_lines;
end
figure
%set (gcf,'Position',[400,400,400,400], 'color','w');
set (gcf, 'color','w');
Y1 = parabolic_Y(1:N_num,2);
%up part
Y2 = Y1 + Y_up_interval;
Y3 = Y1 + 2*Y_up_interval;
Y4 = Y1 + 3*Y_up_interval;
Y5 = Y1 + 4*Y_up_interval;
Y6 = Y1 + 5*Y_up_interval;
%down part
D2 = Y1 - Y_down_interval;
D3 = Y1 - 2*Y_down_interval;
D4 = Y1 - 3*Y_down_interval;
D5 = Y1 - 4*Y_down_interval;
D6 = Y1 - 5*Y_down_interval;
x = 1:N_num;
x = x';
hold all
plot(x,parabolic_Y(1:N_num,3),'Color',[0.1 0.1 0.1]);
plot(x,Y1,'Color',[0.9 0.9 0.9]);
plot(x,Y2,'Color',[0.9 0.9 0.9]);
plot(x,Y3,'Color',[0.9 0.9 0.9]);
plot(x,Y4,'Color',[0.9 0.9 0.9]);
plot(x,Y5,'Color',[0.9 0.9 0.9]);
plot(x,Y6,'Color',[0.9 0.9 0.9]);
fill([1:N_num,fliplr(1:N_num)],[Y1',fliplr(Y2')],'b','FaceColor',[0.1 0.1 0.1],'EdgeColor','none');
fill([1:N_num,fliplr(1:N_num)],[Y2',fliplr(Y3')],'b','FaceColor',[0.3 0.3 0.3],'EdgeColor','none');
fill([1:N_num,fliplr(1:N_num)],[Y3',fliplr(Y4')],'b','FaceColor',[0.5 0.5 0.5],'EdgeColor','none');
fill([1:N_num,fliplr(1:N_num)],[Y4',fliplr(Y5')],'b','FaceColor',[0.7 0.7 0.7],'EdgeColor','none');
fill([1:N_num,fliplr(1:N_num)],[Y5',fliplr(Y6')],'b','FaceColor',[0.9 0.9 0.9],'EdgeColor','none');
plot(x,D2,'Color',[0.9 0.9 0.9]);
plot(x,D3,'Color',[0.9 0.9 0.9]);
plot(x,D4,'Color',[0.9 0.9 0.9]);
plot(x,D5,'Color',[0.9 0.9 0.9]);
plot(x,D6,'Color',[0.9 0.9 0.9]);
fill([1:N_num,fliplr(1:N_num)],[Y1',fliplr(D2')],'b','FaceColor',[0.1 0.1 0.1],'EdgeColor','none');
fill([1:N_num,fliplr(1:N_num)],[D2',fliplr(D3')],'b','FaceColor',[0.3 0.3 0.3],'EdgeColor','none');
fill([1:N_num,fliplr(1:N_num)],[D3',fliplr(D4')],'b','FaceColor',[0.5 0.5 0.5],'EdgeColor','none');
fill([1:N_num,fliplr(1:N_num)],[D4',fliplr(D5')],'b','FaceColor',[0.7 0.7 0.7],'EdgeColor','none');
fill([1:N_num,fliplr(1:N_num)],[D5',fliplr(D6')],'b','FaceColor',[0.9 0.9 0.9],'EdgeColor','none');
plot(x,parabolic_y(1:N_num),'r-')
set(gca,'FontName','Times New Roman','linewidth',1.5,'fontsize',15);
set(gca,'LooseInset',get(gca,'TightInset'))
xlabel('target output','FontName','Times New Roman','fontsize',18);
ylabel('fuzzy granular output','FontName','Times New Roman','fontsize',18);
box off
hold off
rc_Q_d2 = load('uniform_interval_r_c_cov_sp_Q_D2.txt')
%rc_Q_d2 = load('uniform_fuzzy_r_c_cov_sp_Q_D2.txt')
%rc_Q_d2 = load('uniform_parabolic_r_c_cov_sp_Q_D2.txt')
r1 = rc_Q_d2(1:7,end-1);
r2 = rc_Q_d2(8:14,end-1);
r3 = rc_Q_d2(15:21,end-1);
r4 = rc_Q_d2(22:28,end-1);
r5 = rc_Q_d2(29:35,end-1);
r6 = rc_Q_d2(36:42,end-1);
c_list = 1:7;
c_list = c_list';
figure
set(gcf,'Position',[400,240,400,240], 'color','w');
set(gca,'FontName','Times New Roman','linewidth',1.0,'fontsize',15);
set(gca,'LooseInset',get(gca,'TightInset'))%只要这一个命令就可以了
hold all
plot(c_list,r1,'Color',[0.90 0.90 0.90],"Marker","+","LineWidth",1.5)
plot(c_list,r2,'Color',[0.75 0.75 0.75],"Marker","*","LineWidth",1.5)
plot(c_list,r3,'Color',[0.60 0.60 0.60],"Marker","o","LineWidth",1.5)
plot(c_list,r4,'Color',[0.45 0.45 0.45],"Marker","d","LineWidth",1.5)
plot(c_list,r5,'Color',[0.30 0.30 0.30],"Marker","s","LineWidth",1.5)
plot(c_list,r6,'Color',[0.15 0.15 0.15],"Marker","p","LineWidth",1.5)
xlabel('\rmc','FontName','Times New Roman','fontsize',15);
ylabel('\itQ','FontName','Times New Roman','fontsize',15);
legend({'\rmr=1','\rmr=2','\rmr=3','\rmr=4','\rmr=5','\rmr=6'},'NumColumns',3)
legend boxoff
axis tight
grid off
三维图
A=[2 3 14.2743892 2.1 3.15 14.2642229 2.2 3.3 14.2108503 2.3 3.45 14.2883674 2.4 3.6 13.8931573 2.5 3.75 13.6974583 2.6 3.9 13.4649072 2.7 4.05 13.4445747 2.8 4.2 13.0671555 2.9 4.35 13.4242424 3 4.5 13.142131 3.1 4.65 13.3619744 3.2 4.8 13.1726296 3.3 4.95 13.4750732 3.4 5.1 12.67957 3.5 5.25 12.4635385 3.6 5.4 12.8663735 3.7 5.55 12.8257085 3.8 5.7 12.0454545 3.9 5.85 11.7480939 4 6 12.4622677 4.1 6.15 11.2245356 4.2 6.3 11.3668622 4.3 6.45 12.3059629 4.4 6.6 11.7595308 4.5 6.75 12.2627566 4.6 6.9 11.5371455 4.7 7.05 10.8369502 4.8 7.2 11.4087977 4.9 7.35 10.47478 5 7.5 10.9424243]; x=A(:,1);y=A(:,2);z=A(:,3); scatter(x,y,5,z)%散点图 figure [X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x))',linspace(min(y),max(y)),'v4');%插值 pcolor(X,Y,Z);shading interp%伪彩色图 figure,contourf(X,Y,Z) %等高线图 figure,mesh(X,Y,Z)%三维曲面
x = 0:0.001:4;
y = exp(x) + sin(1./(x-0.5));
% sub region to enlarge
x1 = 0.4:0.001:0.6;
y1 = xp(x1) + sin(1./(x1-0.5));
% main figure
figure
plot(x,y)
grid on
rectangle('Position',[min(x1) 0.5 max(x1)-min(x1) 2.5],'EdgeColor','r')% mark the region
% sub figure
sub = axes('Position',[0.2 0.25 0.3 0.25]) % sub figure location
plot(x1,y1)
xlim([min(x1) max(x1)])
set(sub, 'xtick',[],'ytick',[])
