matlab绘图+计算心得

Matlab 常用于做一些计算实验和绘图，我主要用matlab来绘图件。Matlab 可以用GUI界面来完成绘图，但是既然能用鼠标完成的绘图，那么一定能用代码命令来完成绘图，个人比较喜欢用命令来完成matlab的绘图。

保存图片格式

有很多图片格式可以保存，矢量图片格式(".eps")可以不失真，另外，新版matlab还有“copy figure可用”，但是这项功能是针对小的plot可以，有大数据量的图片或plot好像不能不会成功，画密度很高的等值线时，线密度过高则不会成功

保存命令

imwriteas是保存图片, 比如说jpg,bmp； saveas是保存fig,必须用matlab才能打开

matlab绘制图形图像大小的命令：

figure;
plot(diff, t, \'b\',\'linewidth\',1);
set(gca, \'ydir\', \'reverse\');
set(gca, \'xaxislocation\', \'top\');
set(gcf,\'color\',\'white\');
set(gcf,\'unit\',\'centimeters\',\'position\',[5 3 20 13]);
set(gca,\'Position\',[.105 .04 .87 .82],\'linewidth\',2,\'fontsize\',15);
set(gca,\'linewidth\',1,\'fontsize\',15);

xlim([100 1940]);
xlabel(\'Receiver locations (m)\',\'fontsize\',20);
ylabel(\'Time (s)\',\'fontsize\',20);
shg;

print(\'-depsc2\', \'-r600\', \'shot21\');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

matlab 计算精度显示问题：

如：一个数字32774.48，现在在Workspace里面显示出来都是 3.2774e+04 ，想要的只是 32774.48。

则可用vpa(32774.48,7)，另外在输出文件的时候，可通过format , fprintf来控制输出多少位的，详情查一查help命令

将matlab中数据输出保存为txt或dat格式

第一种方法：save（最简单基本的）
具体的命令是：用save *.txt -ascii x
x为变量
*.txt为文件名,该文件存储于当前工作目录下，再打开就可以打开后,数据有可能是以指数形式保存的.
例子：
a =[17 24 1 8 15;23 5 7 14 16 ;4 6 13 20 22 ;10 12 19 21 3 ;11 18 25 2 9 ]；
save afile.txt -ascii a
afile.txt打开之后，是这样的：
1.7000000e+001 2.4000000e+001 1.0000000e+000 8.0000000e+000 1.5000000e+001
2.3000000e+001 5.0000000e+000 7.0000000e+000 1.4000000e+001 1.6000000e+001
4.0000000e+000 6.0000000e+000 1.3000000e+001 2.0000000e+001 2.2000000e+001
1.0000000e+001 1.2000000e+001 1.9000000e+001 2.1000000e+001 3.0000000e+000
1.1000000e+001 1.8000000e+001 2.5000000e+001 2.0000000e+000 9.0000000e+000
第二种方法：dlmwrite
dlmwrite(\'a.txt\',a,\'precision\',\'%10.0f\')
或者是dlmwrite(\'a.txt\',a,\'delimiter\', \'\t\')
对于只有一行或者一列的数据，很适用，但是多行的，就乱了
网上有很多这一类似的问题，但是都不是很理想
第三种方法：fopen+fprintf
下面主要介绍这种方法，可以解决以上问题：用fprintf命令：以上面的例子为例：
第一种情况：
>> a=[17 24 1 8 15;23 5 7 14 16 ;4 6 13 20 22 ;10 12 19 21 3 ;11 18 25 2 9 ];
>> fid = fopen(\'b.txt\',\'wt\');
fprintf(fid,\'%g\n\',a);       # \n 换行
fclose(fid);
然后用写字板打开b.txt，内容如下：为列向量
17
23
4
10
11
24
5
6
12
18
1
7
13
19
25
8
14
20
21
2
15
16
22
3
9
第二种情况：
对上面的命令做一下改动：# \n 换行改为\t，table键
>> fid = fopen(\'b.txt\',\'w\');
fprintf(fid,\'%g\t\',a);
fclose(fid);
然后用写字板打开b.txt，内容如下：为行向量：

17 23 4 10 11 24 5 6 12 18 1 7 13 19 25 8 14 20 21 2 15 16 22 3 9
第三种情况：
综合上面的两个结果，我们编写以下命令：
<A style="key: ih.baidu.comfid=fopen(\'b.txt\',\'wt\');%写入文件路径
[m,n]=size(a);
for i=1:1:m
    for j=1:1:n
       if j==n
         fprintf(fid,\'%g\n\',a(i,j));
      else
        fprintf(fid,\'%g\t\',a(i,j));
       end
    end
end
fclose(fid);
然后用写字板打开b.txt，内容如下：矩阵
17 24 1 8 15
23 5 7 14 16
4 6 13 20 22
10 12 19 21 3
11 18 25 2 9

说明：以上操作都是在当前的工作目录下完成！下面给出最一般的模型，大家可以试着自己操作,如果需要dat格式，直接把txt换为dat就可以
fid=fopen(\'C:\Documents and Settings\cleantotal.ped\',\'wt\');%写入文件路径
matrix=input_mattrix                        %input_matrix为待输出矩阵
[m,n]=size(matrix);
for i=1:1:m
   for j=1:1:n
      if j==n
        fprintf(fid,\'%g\n\',matrix(i,j));
     else
       fprintf(fid,\'%g\t\',matrix(i,j));
      end
   end
end
fclose(fid);