24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1347  |  回复: 1

在捂汗

铁虫 (小有名气)

[求助] matlab中怎么给四条曲线作出图例啊?求大神指导 已有1人参与

clear all; clc
disp('插值求取平滑曲线后的数值' )
t = [5 10 15 20 30 40 60 80 120] ;
A = [29.84 40.01 46.71 53.87 63.32 67.01 67.33 67.54 67.87];

knots = 3; K = 3 ;          % 三次B样条
sp = spap2(knots,K,t,A);
pp = fnder(sp);         % 计算B样条函数的导函数
dAdt = fnval(pp,t);    % 计算t处的导函数值
rAm = dAdt;

% 绘制浓度拟合曲线
ti = linspace(t(1),t(end),200);
Ai = fnval(sp,ti);
plot(t,A,'ro',ti,Ai,'b:')
hold on
t = [5 10 15 20 30 40 60 80 120] ;
B = [36.32 48.21 55.81 62.21 69.81 72.39 72.72 72.91 73.92];
knots = 3; K = 3 ;          % 三次B样条
sp = spap2(knots,K,t,B);
pp = fnder(sp);         % 计算B样条函数的导函数
dBdt = fnval(pp,t);    % 计算t处的导函数值
rBm = dBdt;
ti = linspace(t(1),t(end),200);
Bi = fnval(sp,ti);
plot(t,B,'rx',ti,Bi,'b:')
hold on
t = [5 10 15 20 30 40 60 80 120] ;
C = [41.91 55.62 64.13 70.46 77.64 79.80 80.04 80.57 83.92];

knots = 3; K = 3 ;          % 三次B样条
sp = spap2(knots,K,t,C);
pp = fnder(sp);         % 计算B样条函数的导函数
dCdt = fnval(pp,t);    % 计算t处的导函数值
rCm = dCdt;

% 绘制浓度拟合曲线
ti = linspace(t(1),t(end),200);
Ci = fnval(sp,ti);
plot(t,C,'rs',ti,Ci,'b:')
hold on
t = [5 10 15 20 30 40 60 80 120] ;
D = [42.16 62.13 72.91 79.87 84.21 85.61 85.84 86.21 86.22];

knots = 3; K = 3 ;          % 三次B样条
sp = spap2(knots,K,t,D);
pp = fnder(sp);         % 计算B样条函数的导函数
dDdt = fnval(pp,t);    % 计算t处的导函数值
rDm = dDdt;

% 绘制浓度拟合曲线
ti = linspace(t(1),t(end),200);
Di = fnval(sp,ti);
plot(t,D,'rv',ti,Di,'b:')
Xlabel('t')
ylabel('X_A')
legend(o,'100%',x,'90%',s,'80%',v,'70%')
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)

Matlab专家

感谢参与,应助指数 +1
本帖仅楼主可见
2楼2015-03-09 20:48:58
已阅   申请程序强帖   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 在捂汗 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 7/350 2026-05-17 19:42 by Equinoxhua
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-17 19:39 by Equinoxhua
[考博] 找博士生导师 +4 小代想上岸 2026-05-15 5/250 2026-05-17 17:44 by 脱颖而出
[考博] 2026博士还有哪些学校有名额 +6 小王求读研 2026-05-15 7/350 2026-05-17 16:54 by 知音湖畔
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +13 dd921ww 2026-05-12 16/800 2026-05-17 12:38 by 于轩
[考博] 26/27博士推荐 +3 1木头人13949 2026-05-13 3/150 2026-05-17 09:41 by YuY66
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:06 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:31 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 8/400 2026-05-17 07:26 by 11n4dfd8yn
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:16 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 7/350 2026-05-17 06:55 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:45 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-16 17:45 by j6b2pdz07o
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[文学芳草园] 风把牡丹吹跑了 +5 myrtle 2026-05-12 9/450 2026-05-15 15:27 by myrtle
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
请填处理意见