请选择 进入手机版 | 继续访问电脑版
设为首页收藏本站

52matlab技术论坛

 找回密码
 立即注册

QQ登录

只需一步,快速开始

扫一扫,访问微社区

搜索
热搜: 活动 交友 discuz
查看: 63|回复: 0

对流-扩散-反应方程的Runge-Kutta-Chebyshev算法

[复制链接]
  • TA的每日心情
    开心
    2015-3-12 11:57
  • 签到天数: 1 天

    [LV.1]初来乍到

    63

    主题

    116

    帖子

    904

    积分

    版主

    Rank: 7Rank: 7Rank: 7

    积分
    904
    发表于 2017-11-24 18:08:42 | 显示全部楼层 |阅读模式
    matlab脚本函数内容如下:

    clear all
    close all
    clc

    global f_cout
    f_cout = 0;

    t_start = cputime;
    %
    % Discretize in space
    %
    % m is the number of segments in x-axis
    m=100;
    xs=[0 1];
    mk=m+1;

    % x is a vector that contains the value of points in x axis
    x=linspace(xs(1),xs(2),mk);
    %

    % Give the initial value
    %        w0 is a vector that contains the initial value of w
    w0=sin(x);

    %        w0(1)=0 and w0(mk)=0 are boundary conditions
    w0(mk)=0;

    f = 'F';
    tspan = linspace(0,0.5,mk);
    tol = 1e-6;
    s = 4;
    [w,t]=rkc1(f,tspan,tol,w0',s,x');
    cputime_t = cputime - t_start;

    disp(['cpu时间:',num2str(cputime_t),'s'])
    disp(['F赋值次数',num2str(f_cout)])


    本帖子中包含更多资源

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

    x
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 立即注册

    本版积分规则

    Archiver|手机版|小黑屋|52matlab技术论坛    

    GMT+8, 2018-2-26 09:55 , Processed in 0.351526 second(s), 28 queries .

    Powered by Discuz! X3.2 Licensed

    © 2001-2013 Comsenz Inc.

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