爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6857|回复: 11

[程序设计] matlap程序,出图

[复制链接]

新浪微博达人勋

发表于 2015-4-23 21:18:21 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
各位亲爱的老师,大神们,能不能帮我看一个,我这个程序是出的什么图和怎么用的
clear;
data=xlsread('E:\date\haikou.xlsx');
line=data(:,1);
k=find(line>=254);

for i=1:1:1098
    if i==1098
        n1=k(i,1)+4;
        n2=23670;
    else
        n1=k(i,1)+4;
        n2=k(i+1,1)-1;
        
    end;
    bj{i}=data(n1:n2,:);
end;
for i=3:3:1097
    bj{i}=bj{i+1};
end;
bj{1098}=nan;
% date
load e:\work\station\zdz1.mat
for i=1:1:366*24
    date(i,1)=str2num(zdz1{i,2}(1:8));
end;
date1=reshape(date,24,366);
date2=date1(1,:)';

for i=1096:3:1098
    a1=bj{i};
    a2=bj{i+1};
    a3=bj{i+2};
    m1=find(a1==99999);
    a1(m1)=nan;
    m2=find(a2==99999);
    a2(m2)=nan;
    m3=find(a3==99999);
    a3(m3)=nan;
    f1=find(a1(:,2)==3000);
    f2=find(a2(:,2)==3000);
%     f3=find(a3(:,2)==3000);
    [m1,n1]=size(a1);
    [m2,n2]=size(a2);
    [m3,n3]=size(a3);
    if length(f1)==0
        f1=m1;
    end;
    if f2==0
        f2=m2;
    end;
%     if f3==0
%         f3=m3;
%     end;
    d1=a1([1:f1],:);
    d2=a2([1:f2],:);
%     d3=a3([1:f3],:);
    %     wd
    figure
    plot(d1(:,6),d1(:,2)/10,'-k','LineWidth',2,'MarkerFaceColor','w')
    hold on
    plot(d2(:,6),d2(:,2)/10,'--k','LineWidth',2,'MarkerFaceColor','w')
    hold on
%     plot(d3(:,6),d3(:,2)/10,'-^k','LineWidth',2,'MarkerFaceColor','w')
    set(gca,'xlim',[0 360],'xtick',[0:90:360],'FontSize',15)
    set(gca,'ylim',[300 1100],'ytick',[300:100:1100],'FontSize',15)
    set(gca,'YDir','reverse');
    grid on
%     legend('00','12','00')
    ii=(i-1)/3+1;
    s=strcat(num2str(date2(ii)),'风向')
    title(s)
    saveas(gcf,s,'png')
    close(gcf)
   
%         % ws
%         figure
%         plot(d1(:,7)/10,d1(:,2)/10,'-k','LineWidth',2,'MarkerFaceColor','w')
%         hold on
%         plot(d2(:,7)/10,d2(:,2)/10,'--k','LineWidth',2,'MarkerFaceColor','w')
%         hold on
% %         plot(d3(:,7)/10,d3(:,2)/10,'-^k','LineWidth',2,'MarkerFaceColor','w')
%         set(gca,'xlim',[0 36],'xtick',[0:2:36],'FontSize',15)
%         set(gca,'ylim',[300 1100],'ytick',[300:100:1100],'FontSize',15)
%         set(gca,'YDir','reverse');
%         grid on
% %         legend('00','12','00','Location','NorthEast')
%         ii=(i-1)/3+1;
%         s=strcat(num2str(date2(ii)),'风速')
%         title(s)
%         saveas(gcf,s,'png')
%         close(gcf)
end;

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2015-4-24 15:56:34 | 显示全部楼层
数据都没有你说个毛线……
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2015-4-24 15:57:16 | 显示全部楼层
还有,我没听说过过matlap…你找一个合适的板块问问吧
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-27 09:43:23 | 显示全部楼层
二爷名声在外 发表于 2015-4-24 15:56
数据都没有你说个毛线……

数据是网上下的,说是数据读取程序,能不能我把数据给你,帮我看一下吗??求求你拉,因为一点都看不出来,更不用说拿数据,编程,出图了,求你拉,帮帮忙,老师
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-27 09:44:39 | 显示全部楼层
二爷名声在外 发表于 2015-4-24 15:57
还有,我没听说过过matlap…你找一个合适的板块问问吧

好吧  谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-27 09:47:34 | 显示全部楼层
yx0323 发表于 2015-4-27 09:43
数据是网上下的,说是数据读取程序,能不能我把数据给你,帮我看一下吗??求求你拉,因为一点都看不出来 ...

把数据以附件的形式挂出来,
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2015-4-27 09:53:16 | 显示全部楼层
yx0323 发表于 2015-4-27 09:43
数据是网上下的,说是数据读取程序,能不能我把数据给你,帮我看一下吗??求求你拉,因为一点都看不出来 ...

你的问题就好像给我一个电饭锅,问我怎么做八宝粥,不告诉我你是用什么米作为原料~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-28 11:06:15 | 显示全部楼层
二爷名声在外 发表于 2015-4-27 09:53
你的问题就好像给我一个电饭锅,问我怎么做八宝粥,不告诉我你是用什么米作为原料~

老师,你好
上次我发的是导师给我用的读取资料的一个参考,我现在是读取不了
我的资料是探空网站上下的下面那种,下面是一天的资料,所以,麻烦老师帮我一下[em:16:]
254      0      1      JAN    2012
      1  99999  57957  25.33N110.30E   166  99999
      2   7000    700    743     23  99999      3
      3          ZGKL                99999     kt
      9  10060    166     84     40     25      2
      4  10000    216     80     37     25      6
      4   9250    855     46     28     50     12
      4   8500   1544     46     33    155      2
      4   7000   3120    -13    -38    280     21
      4   5000   5810    -57   -477    275     35
      4   4000   7520   -173   -553    290     64
      4   3000   9610   -319   -599    285     65
      4   2500  10880   -403   -663    285     71
      8   2500  10880  99999  99999    285     71
      4   2000  12360   -503   -743    285     67
      4   1500  14180   -649  99999    275     64
      4   1000  16560   -777  99999    275     50
      7    743  18231   -834  99999    280     25
      4    700  18580   -761  99999    270     25
      4    500  20550   -673  99999    260     17
      4    300  23680   -595  99999    275     21
      4    200  26280   -501  99999    335     10
      4    100  30800   -489  99999     35     29
    254     12      1      JAN    2012
      1  99999  57957  25.33N110.30E   166  99999
      2    100    700    808     23  99999      3
      3          ZGKL                99999     kt
      9  10050    166    100     40     25      4
      4  10000    204     98     38     25      6
      4   9250    846     52     29     50     13
      4   8500   1537     46     33    115      2
      4   7000   3111    -15    -41    280     21
      4   5000   5800    -33   -463    290     40
      4   4000   7520   -159   -549    290     48
      4   3000   9620   -331   -471    280     60
      6   2710  10322  99999  99999    275     69
      4   2500  10880   -407   -717    275     67
      4   2000  12370   -527   -807    285     62
      4   1500  14170   -659  99999    280     69
      4   1000  16540   -785  99999    275     44
      7    808  17741   -821  99999    260     31
      4    700  18550   -781  99999    265     31
      4    500  20520   -667  99999    260     19
      4    300  23660   -607  99999    265     27
      4    200  26230   -517  99999    335     17
      4    100  30760   -509  99999     65     17
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-28 11:30:01 | 显示全部楼层
denghaijun 发表于 2015-4-27 09:47
把数据以附件的形式挂出来,

老师,你好
上次我发的是导师给我用的读取资料的一个参考,我现在是读取不了
我的资料是探空网站上下的下面那种,下面是一天的资料,所以,麻烦老师帮我一下[em:16:]
254      0      1      JAN    2012
      1  99999  57957  25.33N110.30E   166  99999
      2   7000    700    743     23  99999      3
      3          ZGKL                99999     kt
      9  10060    166     84     40     25      2
      4  10000    216     80     37     25      6
      4   9250    855     46     28     50     12
      4   8500   1544     46     33    155      2
      4   7000   3120    -13    -38    280     21
      4   5000   5810    -57   -477    275     35
      4   4000   7520   -173   -553    290     64
      4   3000   9610   -319   -599    285     65
      4   2500  10880   -403   -663    285     71
      8   2500  10880  99999  99999    285     71
      4   2000  12360   -503   -743    285     67
      4   1500  14180   -649  99999    275     64
      4   1000  16560   -777  99999    275     50
      7    743  18231   -834  99999    280     25
      4    700  18580   -761  99999    270     25
      4    500  20550   -673  99999    260     17
      4    300  23680   -595  99999    275     21
      4    200  26280   -501  99999    335     10
      4    100  30800   -489  99999     35     29
    254     12      1      JAN    2012
      1  99999  57957  25.33N110.30E   166  99999
      2    100    700    808     23  99999      3
      3          ZGKL                99999     kt
      9  10050    166    100     40     25      4
      4  10000    204     98     38     25      6
      4   9250    846     52     29     50     13
      4   8500   1537     46     33    115      2
      4   7000   3111    -15    -41    280     21
      4   5000   5800    -33   -463    290     40
      4   4000   7520   -159   -549    290     48
      4   3000   9620   -331   -471    280     60
      6   2710  10322  99999  99999    275     69
      4   2500  10880   -407   -717    275     67
      4   2000  12370   -527   -807    285     62
      4   1500  14170   -659  99999    280     69
      4   1000  16540   -785  99999    275     44
      7    808  17741   -821  99999    260     31
      4    700  18550   -781  99999    265     31
      4    500  20520   -667  99999    260     19
      4    300  23660   -607  99999    265     27
      4    200  26230   -517  99999    335     17
      4    100  30760   -509  99999     65     17
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2015-4-28 15:17:31 | 显示全部楼层
不全,你还laod了一个文件
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表