建模浓度使用颂歌

2视图(30天)

显示旧的评论

鲜明的伏 2021年9月26日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/1461089-modeling-concentration-using-ode

回答: Sulaymon Eshkabilov 2021年9月29日

test_data.zip

我想模型染料浓度,但我不能让我卡住了。它让一个错误但不会告诉什么是错。任何一个能帮助谢谢吗?

                          清晰的所有;
                         
                          data =负载(“test_data.dat”);
                         
                             错误使用负载
找不到文件或目录“test_data.dat”。

                          时间=大小(数据,1);
                         
                          u =数据(:,2);
                         
                          tspan =(1:1:时间);
                         
                          y0 =数据(1、3);
                         
                          f = -0.25;
                         
                          g = 0.25;
                         
                          T = 10.0;
                         
                          [t、y] =数值(@ model, tspan, y0, f, g, t, U);
                         
                          情节(t y(:,结束),“* - b”,“线宽”1);
                         
                          函数dydt =模型(t、y, f, g, t, u)
                         
                          如果T - 1 T < =
                         
                          u = 0;
                         
                          其他的
                         
                          t_back = t t;
                         
                          n =地板(t_back);u = (n);
                         
                          结束
                         
                          dydt = f * y + g * u;
                         
                          返回;
                         
                          结束

9日评论
显示8年长的评论隐藏8年长的评论

图像分析 2021年9月26日

选择是什么?它说它不是定义。如果我使用选项而不是选项,我得到更多的错误。

                                clc;%清除命令窗口。
                               
                                关闭所有;%关闭所有数据(除了imtool。)
                               
                                imtool关闭所有;%关闭所有imtool数字图像处理工具箱。
                               
                                清晰;%擦掉所有现有的变量。或者clearvars如果你想要的。
                               
                                工作空间;%确定工作区面板显示。
                               
                                格式短g;
                               
                                格式紧凑的;
                               
                                字形大小= 20;
                               
                                % (.mlx文件)
                               
                                %染料浓度
                               
                                清晰的所有;
                               
                                data =负载(“test_data.dat”);
                               
                                时间=大小(数据,1);%的仿真时间
                               
                                输入=数据(:,2);%最高浓度测量
                               
                                rslt =数据(:,3);%浓度测量底部
                               
                                tspan =(1:1:时间);%时间变量
                               
                                公司=数据(1、3);%初始浓度下
                               
                                %模型参数
                               
                                fparameter = -0.25;
                               
                                gparameter = 0.25;
                               
                                延迟= 10.0;
                               
                                选择= odeset (“RelTol”1 e - 3);
                               
                                [t、y] =数值(@ model, tspan, inc,选项);
                               
                                情节(tspan obsv,“- b”,“线宽”,0.5);持有在;
                               
                                情节(tspan输入,“g”,“线宽”,0.5);
                               
                                情节(t y(:,结束),“m”,“线宽”1);
                               
                                传奇(“rslt”,“输入”);
                               
                                包含(的时间(秒));ylabel (“浓度”);
                               
                                网格;
                               
                                %预定义函数(m文件)
                               
                                函数dydt =模型(t y参数,输入)
                               
                                %模型常微分方程(ODE)模型:
                               
                                % dy / dt = f * y (t) + g * u (t t)
                               
                                dxdt = 0;%返回变量
                               
                                如果t - 1(3) < =参数%延迟初始浓度
                               
                                u_input = 0;
                               
                                其他的
                               
                                t_back = t-parameter (3);
                               
                                n =地板(t_back);u_input =输入(n);
                               
                                结束
                               
                                dydt = f * y (t) + g * u (t t)
                               
                                返回;
                               
                                结束