rundemo.m

来自「MDPSAS工具箱是马里兰大学开发的」· M 代码 · 共 32 行

% Test the catpel class

addpath ../../class
addpath ../../class/double

clear

tm = 20; % Thiele modulus

figure(1)

A = catpel(tm,1);
A = set(A,'solverset',1,'plot');
A = newton(A);

unpack(A)
conc1 = a;

B = catpel(tm+1,1);
B = set(B,'solverset',1,'plot');
B = newton(B);
unpack(B)
conc2 = a;

S = sensitivity(A,'tm');
J = jacobian(A);

dcdtm = scalarfield(R,-J\S);

figure(2)
plot(conc2-conc1,dcdtm)
title('Sensitivity w.r.t. Thiele modulus')