为什么不心动周期不同电阻和电容的代码运行?

6视图(30天)

显示旧的评论

蒂莫西·珍罗哈斯 2021年11月11日

1
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/1583894-why-doesn-t-the-code-for-cardiac-cycle-with-varying-resistance-and-capacitance-run

评论道: 威廉•罗斯 2021年11月11日

答:接受威廉•罗斯

每当我试图运行代码,出现这个提示:

“函数定义在脚本中必须出现在文件的末尾。

后将所有语句“正弦曲线”函数定义在第一个地方

函数定义”

之前我试图添加函数相等的语句,但我不知道怎么做。

请帮助修复代码。

这是代码:

% %正弦流入

% clc;

%清除所有

R = 1;%血管阻力mmHg-s /毫升为健康的人

C = 1;%电容在ml /毫米汞柱

tc = 0.8;%心动周期的长度(s)

[time_si, pressure_si] =正弦信号(R、C、tc);

% %分析

sys_p = max (pressure_si);%收缩压

dia_p = min (pressure_si);%舒张压

地图= 1/3 * sys_p + 2/3 * dia_p;%平均动脉压

pulse_pressure = sys_p - dia_p;

% %改变电阻

clc

%抗性提高25%

R1 = 1.25;

C = 1;

[time_r1, pressure_r1] =正弦信号(R1, C, tc);

sys_p1 = max (pressure_r1)

dia_p1 = min (pressure_r1);

map1 = 1/3 * sys_p1 + 2/3 * dia_p1;

pulse_pressure_r1 = sys_p1 - dia_p1

%的阻力减少25%

R2 = 0.75;

C = 1;

[time_r2, pressure_r2] =正弦信号(R2、C、tc);

sys_p2 = max (pressure_r2)

dia_p2 = min (pressure_r2)

map2 = 1/3 * sys_p2 + 2/3 * dia_p2;

pulse_pressure_r2 = sys_p2 - dia_p2

lw = 3;

fs = 18;

图

抓住

网格

情节(time_si pressure_si,“g”,“线宽”,lw)

情节(time_r1 pressure_r1,“r”,“线宽”,lw)

情节(time_r2 pressure_r2,“b”,“线宽”,lw)

传奇(“正常R, R增长25%”,“R下降25%”)

集(gca、字形大小,fs)

ylim (150 [0])

包含(“时间[s]”)

ylabel(心室压力(毫米汞柱))

% %改变电容

clc

%电容提高25%

R = 1;

C1 = 1.25;

[time_c1, pressure_c1] =正弦信号(R, C1, tc);

sys_p1 = max (pressure_c1)

dia_p1 = min (pressure_c1);

map1 = 1/3 * sys_p1 + 2/3 * dia_p1;

pulse_pressure_c1 = sys_p1 - dia_p1

%电容减少25%

R = 1;

C2 = 0.75;

[time_c2, pressure_c2] =正弦信号(R, C2, tc);

sys_p2 = max (pressure_c2)

dia_p2 = min (pressure_c2)

map2 = 1/3 * sys_p2 + 2/3 * dia_p2;

pulse_pressure_c2 = sys_p2 - dia_p2

lw = 3;

fs = 18;

图

抓住

网格

情节(time_si pressure_si,“g”,“线宽”,lw)

情节(time_c1 pressure_c1,“r”,“线宽”,lw)

情节(time_c2 pressure_c2,“b”,“线宽”,lw)

传奇(正常的C, C增加了25%,“C下降25%”)

集(gca、字形大小,fs)

ylim (150 [0])

包含(“时间[s]”)

ylabel(心室压力(毫米汞柱))

% %

图

次要情节(121)

抓住

网格

情节(time_si pressure_si,“g”,“线宽”,lw)

情节(time_r1 pressure_r1,“r”,“线宽”,lw)

情节(time_r2 pressure_r2,“b”,“线宽”,lw)

传奇(“正常R, R增长25%”,“R下降25%”)

集(gca、字形大小,fs)

ylim (150 [0])

包含(“时间[s]”)

ylabel(心室压力(毫米汞柱))

标题(电阻对压力的影响)

次要情节(122)

抓住

网格

情节(time_si pressure_si,“g”,“线宽”,lw)

情节(time_c1 pressure_c1,“r”,“线宽”,lw)

情节(time_c2 pressure_c2,“b”,“线宽”,lw)

传奇(正常的C, C增加了25%,“C下降25%”)

集(gca、字形大小,fs)

ylim (150 [0])

包含(“时间[s]”)

ylabel(心室压力(毫米汞柱))

标题(电容对压力的影响)

0评论
显示1年长的评论藏1年长的评论

登录置评。

在回答这个问题。

接受的答案

威廉•罗斯 2021年11月11日

1
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/1583894-why-doesn-t-the-code-for-cardiac-cycle-with-varying-resistance-and-capacitance-run answer_828969

@Timothy Jen罗哈斯 ,

我救了你的程序脚本CVSimMatlabCentral.m。我跑,结果:

> > CVSimMatlabCentral

无法执行任务,因为点索引不支持这种类型的变量。万博1manbetx

错误在正弦信号(第32行)

mstruct。mapparallels = 0;

错误CVSimMatlabCentral(11行)

[time_si, pressure_si] =正弦信号(R、C、tc);

这里的问题是正弦曲线() 是一个内置的Matlab函数。它是一种地图投影。预计输入不同的比你有提供,因此错误。

我怀疑你写自己的sinsoid()例程,但不是在您发布的代码。如果你将正弦信号()你的脚本,最后,它将使用您的程序,而不是内置正弦信号()。

如果你回复这个帖子,请附上您的脚本作为m文件。

2的评论
显示1年长的评论藏1年长的评论

威廉•罗斯 2021年11月11日

CVSimMatlabCentral.m

@Timothy Jen罗哈斯你估计地图,地图= SysP / 3 + 2 * DiaP / 3。

估计会很不准确的,因为各种各样的原因。在你的模拟,估计是非常不准确的。它大大低估了真正的地图。你可以更准确地计算出平均压力为一个节拍时间平均压力的仿真:

mapA =意味着(pressure_si (time_si > = time_si(结束)tc));%准确的地图

试一试。下面的结果(用附加脚本生成),注意到准确的地图和地图之间的赐予。

                                   > > CVSimMatlabCentral
                                  
                                   Normal__________: SysP = 127.5, DiaP = 74.8, = 92.4地图,MAPa = 100.0,页= 52.7
                                  
                                   电阻+ 25%:SysP = 152.4, DiaP = 99.5, = 117.1地图,MAPa = 125.0,页= 52.9
                                  
                                   电阻-25年%:SysP = 102.7, DiaP = 50.4, = 67.8地图,MAPa = 75.0,页= 52.3
                                  
                                   电容+ 25%:SysP = 121.9, DiaP = 79.6, = 93.7地图,MAPa = 100.0,页= 42.4
                                  
                                   电容-25年%:SysP = 136.9, DiaP = 67.2, = 90.4地图,MAPa = 100.0,页= 69.7

登录置评。

类别

数学和优化映射的工具箱地图显示自定义地图

雷达映射的工具箱地图显示自定义地图

找到更多的在自定义地图在帮助中心和文件交换

s manbetx 845

MATLAB

释放

R2021a

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!