QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 25292|回复: 94
打印 上一主题 下一主题

关于长江水质预测与评价模型的

[复制链接]
字体大小: 正常 放大

16

主题

3

听众

55

积分

升级  52.63%

  • TA的每日心情
    开心
    2012-9-7 07:14
  • 签到天数: 10 天

    [LV.3]偶尔看看II

    自我介绍
    我很踏实刚刚搞数学建模,希望大家互相帮助

    新人进步奖

    群组Matlab讨论组

    群组小草的客厅

    群组岩土力学与地下工程

    群组学术交流A

    跳转到指定楼层
    1#
    发表于 2010-8-16 13:23 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    请各位大虾帮我看一下这个程序问题出在哪?还有怎么解决?我运行后总是出现出错的信息。




    Nihe=zeros(6,2);
    Weilai=zeros(6,20);
    for j=1:6;
      %figure(j)
    grid on
    hold on
    subplot(3,2,j)
    %figure(j)
    hold on
    load f:\data_a4.txt
    Data_ok=data_a4;
    Data_ok=Data_ok(:,3:2:13);
    %枯水期
    plot(Data_ok(1:9:90,j),'-*b')
    plot(Data_ok(2:9:90,j),':*b')
    plot(Data_ok(3:9:90,j),'-.*b')
    %legend('枯水期全流域','枯水期干流','枯水期支流',)
    %丰水期plot(Data_ok(3:9:90,j),'-.*b')
    plot(Data_ok(4:9:90,j),'or')
    plot(Data_ok(5:9:90,j),'r')
    plot(Data_ok(6:9:90,j),'-.or')
    %legend('丰水期全流域','丰水期干流','丰水期支流',)
    %水文年
    plot(Data_ok(7:9:90,j),'-xg')
    plot(Data_ok(8:9:90,j),':xg')
    plot(Data_ok(9:9:90,j),'-.xg')
    switch j
        case 1
            xlabel('I类水')
        case 2
            xlabel('II类水')
        case 3
            xlabel('III类水')
        case 4
            xlabel('IV类水')
        case 5
             xlabel('V类水')
        case 6
              xlabel('VI类水')
    end

    if j>=3;
    x=1:10;
    X=[ones(1,length(x))',x'];
    y=Data_ok(7:9:90,j);
    [b,bint,r,rint,stats]=regress(y,X,0.05);
    t=1:20;
    Nihe(j,=b;
    yt=b(1)+b(2).*t';
    WeiLai(j,=yt(1:20);
    plot(t,yt,'LineWidth',2,'Color','k')
    end
    if j==1;
        x=1:10;
        X=[ones(1,length(x))',x'];
        y=Data_ok(7:9:90,j);
        [b,bint,r,rint,stats]=regress(log(y),X,0.05)
        t=1:20;
        Nihe(j,=b;
        yt=b(1)+b(2).*t';
        ytt=exp(yt);
        WeiLai(j,=ytt(1:20);
        plot(t,ytt,'LineWidth',2,'Color','k')
    end
    if j==2;
        x2=5:10;
        X2=[ones(1,length(x2))',x2'];
        y2=Data_ok(43:9:90,j);
        [b2,bint2,r2,rint2,stats2]=regress(y2,X2,0.05)
        t2=1:20;
        Nihe(j,=b;
        yt2=b2(1)+b2(2).*t2';
        WeiLai(j,=yt2(1:20);
        plot(t2,yt2,'LineWidth',2,'Color','k')
    end
    aa=[174 179 183 189 207 234 220.5 256 270 285];
    tt=1:10;
    plot(tt,aa/10,'--*c')
    Xf=[ones(1,length(tt))',tt'];
    [bf,bintf,rf,rintf,statsf]=regress(aa',Xf,0.05)
    tf=1:20;
    ytf=bf(1)+bf(2).*tf';
    plot(tf,ytf/10,'LineWidth',1,'Color','r')
    grid on
    end
    ???     X=[ones(1,length(x))',x'];
                     |
    Error: Missing variable or function

    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

    13

    主题

    2

    听众

    1074

    积分

    升级  7.4%

  • TA的每日心情
    无聊
    2013-12-11 13:50
  • 签到天数: 49 天

    [LV.5]常住居民I

    回复

    使用道具 举报

    linmatsas 实名认证       

    53

    主题

    13

    听众

    3591

    积分

    逍遥游

  • TA的每日心情
    奋斗
    2014-12-2 09:53
  • 签到天数: 54 天

    [LV.5]常住居民I

    自我介绍
    额。。。。世界上最讨厌的事情就是自我介绍。。。

    邮箱绑定达人 新人进步奖 发帖功臣 最具活力勋章

    群组Matlab讨论组

    群组数学建模

    群组小草的客厅

    群组2012数学一考研交流

    群组C 语言讨论组

    回复

    使用道具 举报

    1

    主题

    4

    听众

    36

    积分

    升级  32.63%

    该用户从未签到

    自我介绍
    诚恳
    回复

    使用道具 举报

    6

    主题

    4

    听众

    1859

    积分

    升级  85.9%

    该用户从未签到

    回复

    使用道具 举报

    紫辰 实名认证       

    12

    主题

    16

    听众

    1304

    积分

    升级  30.4%

  • TA的每日心情
    擦汗
    2013-2-5 09:29
  • 签到天数: 35 天

    [LV.5]常住居民I

    自我介绍
    200 字节以内

    不支持自定义 Discuz! 代码

    群组学术交流A

    群组数学建模保研联盟

    群组Matlab讨论组

    群组湖南大学数学建模

    群组学术交流B

    回复

    使用道具 举报

    双调 实名认证       

    4

    主题

    3

    听众

    64

    积分

    升级  62.11%

  • TA的每日心情
    无聊
    2012-4-20 23:11
  • 签到天数: 1 天

    [LV.1]初来乍到

    回复

    使用道具 举报

    追鸿        

    0

    主题

    2

    听众

    47

    积分

    升级  44.21%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    0

    主题

    2

    听众

    30

    积分

    升级  26.32%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    zhoulb        

    0

    主题

    2

    听众

    99

    积分

    升级  98.95%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2025-8-14 06:18 , Processed in 1.566790 second(s), 105 queries .

    回顶部