存储变量后循环(容易)

11视图(30天)

显示旧的评论

保罗•罗杰斯 2021年1月8日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/711948-storing-variables-after-a-loop-easy

评论道: 保罗•罗杰斯2021年1月11日

嗨,我设法使一个循环中运行一个函数。

预期的输出是:

t、y (: 1), y (:, 2)

问题是在迭代结束时我只有最后一个循环的一部分,但是

我将男孩的所有数据,对那些向量,在looop计算。

这是部分代码:

                         为i = 1:(长度(giri) 2)
                        
                         U1 (i) = (D1 *π。*义理(我))/ 60;
                        
                         U2 (i) = (D2 *π。*义理(我))/ 60;
                        
                         B (i) = U1 (i)。/ (2。* wh。* Lc);
                        
                         [t、y] = ode113 (@greitzer (t0(我),特遣部队(i)], [y1 (i), y2 (i)],选项,U1 (i), U2 (i), B (i));%我发现ode113方式更有效率
                        
                         t0 (i + 1) = time_rpm (i + 1);%仿真开始[s]
                        
                         特遣部队(i + 1) = time_rpm (+ 2);%模拟完成的[s]
                        
                         日元(i + 1) = y(结束,1);%初始条件1
                        
                         y2 (i + 1) = y (, 2);%初始条件2
                        
                         情节(t、y (:, 2))
                        
                         网格在
                        
                         网格小
                        
                         包含(“t [s]”)
                        
                         ylabel (“\ψ”)
                        
                         持有在
                        
                         结束

我想叫:

                         t =时间
                        
                         y(: 1) =φ
                        
                         y (:, 2) = psi

8的评论
显示隐藏 7年长的评论

马修•诺伊 2021年1月11日

你好,保罗

我已经开始工作在你的文件从其他帖子,但我有一些错误

我看到你分享一个新版本的代码

首先,来加速这个过程,我刚刚改变了这两条线,但是您可以很容易地回到自己的调优:

                               % time_rpm = [25:0.01:40];%原
                              
                               time_rpm = [25:1:40];% MN
                              
                               …。
                              
                               sample_frequency = compressor_frequency * 8;%采样频率(赫兹)(至少选择一个频率

然后,我理解,你想连接所有迭代的ode求解器的结果——这就是我(没有什么花哨)和在最后的阴谋策划(寻找行 % MN )

这是整个代码:

                               清晰的
                              
                               clc
                              
                               关闭所有
                              
                               % %气体性质
                              
                               a01 = 340;%声速的天然气(米/秒)
                              
                               p01 = 1 e5;%压缩机入口压力(Pa)
                              
                               T01 = 303.35;%压缩机进气温度(K)
                              
                               cp = 1005;%为恒压比热容,财产的气体(J /(公斤* K))
                              
                               ro1 = 1.15;%气体密度(公斤/米^ 3)
                              
                               Re = 100000;%雷诺数()
                              
                               k = 1.4;%气体常数比率,cp /简历()
                              
                               % %压缩机特性
                              
                               Lc = 0.53;%压缩机和导管的长度(米)
                              
                               Vp = 0.01;%充气卷(m ^ 3)
                              
                               Ac = 0.0026;
                              
                               % kl = 0.0008;%节流,比例节流阀的开口
                              
                               dt₁= 0.074;%叶轮直径在诱导物提示(m)
                              
                               《= 0.032;%叶轮直径套管(m)在中心
                              
                               D1 = 1 /√(2) * sqrt ((dt₁) ^ 2 + (《) ^ 2);%的平均直径(米)
                              
                               Di = 0.02;%的意思是液压叶轮直径(米)
                              
                               Dd = 0.02;%的意思是液压叶轮直径(米)
                              
                               我= 0.001;%压缩机轴的惯性矩(公斤* m ^ 2)
                              
                               alfa1 =π/ 2;在诱导物,%气流角α=π/ 2 rad,假定没有预旋
                              
                               =((π* D1 ^ 2) / 4) * 4;
                              
                               Ai =((π* D1 ^ 2) / 4);%截面面积叶轮(m ^ 2)
                              
                               广告=((π* D1 ^ 2) / 4);%截面面积扩散器(m ^ 2)
                              
                               李= 0.053;%的意思是通道长度(米)
                              
                               ld = 0.053;%的意思是通道长度(米)
                              
                               beta1b = 0.61;%叶轮叶片进口角()
                              
                               Ch = 4 * 0.3164 *(重新)^ -0.25;%摩擦损失系数()
                              
                               kf = (Ch *李)/(2 *迪* ro1 ^ 2 * Ai ^ 2 * (sin (beta1b)) ^ 2) * 4;%摩擦损失叶轮
                              
                               kfd = (Ch * ld) / (2 * Dd * ro1广告^ ^ 2 * 2 * (sin (beta1b)) ^ 2) * 4;%摩擦损失扩散器
                              
                               σ= 0.9;%滑移系数()
                              
                               D2 = 0.128;%叶轮直径提示(m)
                              
                               alfa2b =每股((D1 * tan (beta1b)) /(σ* D2));
                              
                               m = [-0.2:0.001:0.8];%(公斤/ s)质量流量
                              
                               cn = 10;
                              
                               % %的速度
                              
                               giri_medium = 55000;%节流阀孔径
                              
                               delta_giri = 400;% gamma_T_max e gamma_T_min之间的范围
                              
                               giri_max = giri_medium + (delta_giri / 2);%最大光圈
                              
                               giri_min = giri_medium - (delta_giri / 2);%最小孔径
                              
                               A_rpm = (giri_max - giri_min) / 2;%振幅
                              
                               b_rpm = (giri_max + giri_min) / 2;
                              
                               rpm_frequency = 20;%引擎的频率(赫兹)
                              
                               w_rpm = rpm_frequency * 2 *π;%引擎的频率(弧度/秒)
                              
                               % time_rpm = [25:0.01:40];%原
                              
                               time_rpm = [25:1:40];% MN
                              
                               giri = A_rpm * sin (w_rpm * time_rpm) + b_rpm;
                              
                               wh = a01。* (Ac. /副总裁。* Lc)。^ 0.5;% Helmotz的频率
                              
                               % %效率损失
                              
                               deltanbf = 0.02;%回流损失(0.02)环形扩散- 0,03介质体积- 0.05叶片扩压器
                              
                               deltanv = 0.035;%蜗壳损失()
                              
                               deltanc = 0.00;%间隙损失()
                              
                               deltand = 0.00;%扩散损失()
                              
                               德尔塔= deltanbf + deltanv + deltanc + deltand;
                              
                               % %节流阀的参数
                              
                               gamma_T =。;%节流阀孔径
                              
                               Delta_gamma_T = .02点;% gamma_T_max e gamma_T_min之间的范围
                              
                               gamma_T_max = gamma_T + (Delta_gamma_T / 2);%最大光圈
                              
                               gamma_T_min = gamma_T - (Delta_gamma_T / 2);%最小孔径
                              
                               = (gamma_T_max - gamma_T_min) / 2;%振幅
                              
                               b = (gamma_T_max + gamma_T_min) / 2;
                              
                               % %压缩机的频率(频率)
                              
                               compressor_frequency = 10;%引擎的频率(赫兹)
                              
                               w = compressor_frequency * 2 *π;%引擎的频率(弧度/秒)
                              
                               %保存main_parameters.mat
                              
                               % %颂歌的参数
                              
                               t0 (1) = 0;%仿真开始[s]
                              
                               (1)= time_rpm特遣部队(2);%模拟完成的[s]
                              
                               日元(1)= 0;%初始条件1
                              
                               y2 (1) = 0;%初始条件2
                              
                               init = (y1 y2) ';%初始条件向量
                              
                               sample_frequency = compressor_frequency * 8;%采样频率(赫兹)(至少选择一个频率
                              
                               % 10 - 12倍的分析
                              
                               %定义好。只是快速尝试2 - 3倍)
                              
                               maximum_time_step = 1 / sample_frequency;
                              
                               选择= odeset (“MaxStep”,maximum_time_step);%最大时间步大小
                              
                               保存main_parameters.mat
                              
                               t_all = [];% MN
                              
                               y_all = [];% MN
                              
                               为i = 1:(长度(giri) 2)
                              
                               U1 (i) = (D1 *π。*义理(我))/ 60;
                              
                               U2 (i) = (D2 *π。*义理(我))/ 60;
                              
                               B (i) = U1 (i)。/ (2。* wh。* Lc);
                              
                               [t、y] = ode113 (@greitzer (t0(我),特遣部队(i)], [y1 (i), y2 (i)],选项,U1 (i), U2 (i), B (i));%我发现ode113方式更有效率
                              
                               t0 (i + 1) = time_rpm (i + 1);%仿真开始[s]
                              
                               特遣部队(i + 1) = time_rpm (+ 2);%模拟完成的[s]
                              
                               日元(i + 1) = y(结束,1);%初始条件1
                              
                               y2 (i + 1) = y (, 2);%初始条件2
                              
                               %。内格罗蓬特:concat t y为所有迭代
                              
                               t_all = [t_all;t];% MN
                              
                               y_all = [y_all;y);% MN
                              
                               % array_val = 0(长度(giri) 2);
                              
                               %的阴谋(t y (: 1))
                              
                               %网格
                              
                               %网格小
                              
                               %包含(“t [s]”)
                              
                               % ylabel(“\φ”)
                              
                               %抓住
                              
                               图(%)
                              
                               %的阴谋(t, y (:, 2))
                              
                               %网格
                              
                               %网格小
                              
                               %包含(“t [s]”)
                              
                               % ylabel(“\ψ”)
                              
                               %抓住
                              
                               % savefig (10 hz)
                              
                               结束
                              
                               情节(t0, y1)
                              
                               网格在
                              
                               网格小
                              
                               包含(“t [s]”)
                              
                               ylabel (‘\φ)
                              
                               %抓住
                              
                               图()
                              
                               情节(t0, y2)
                              
                               网格在
                              
                               网格小
                              
                               包含(“t [s]”)
                              
                               ylabel (“\ψ”)
                              
                               %抓住
                              
                               图()
                              
                               情节(t_all y_all)
                              
                               网格在
                              
                               网格小
                              
                               包含(“t [s]”)
                              
                               ylabel (“y”)

马修•诺伊 2021年1月11日

确定

很高兴这个小插件解决你的问题

我把相同的代码再下面你可以接受它

为未来好运

                               清晰的
                              
                               clc
                              
                               关闭所有
                              
                               % %气体性质
                              
                               a01 = 340;%声速的天然气(米/秒)
                              
                               p01 = 1 e5;%压缩机入口压力(Pa)
                              
                               T01 = 303.35;%压缩机进气温度(K)
                              
                               cp = 1005;%为恒压比热容,财产的气体(J /(公斤* K))
                              
                               ro1 = 1.15;%气体密度(公斤/米^ 3)
                              
                               Re = 100000;%雷诺数()
                              
                               k = 1.4;%气体常数比率,cp /简历()
                              
                               % %压缩机特性
                              
                               Lc = 0.53;%压缩机和导管的长度(米)
                              
                               Vp = 0.01;%充气卷(m ^ 3)
                              
                               Ac = 0.0026;
                              
                               % kl = 0.0008;%节流,比例节流阀的开口
                              
                               dt₁= 0.074;%叶轮直径在诱导物提示(m)
                              
                               《= 0.032;%叶轮直径套管(m)在中心
                              
                               D1 = 1 /√(2) * sqrt ((dt₁) ^ 2 + (《) ^ 2);%的平均直径(米)
                              
                               Di = 0.02;%的意思是液压叶轮直径(米)
                              
                               Dd = 0.02;%的意思是液压叶轮直径(米)
                              
                               我= 0.001;%压缩机轴的惯性矩(公斤* m ^ 2)
                              
                               alfa1 =π/ 2;在诱导物,%气流角α=π/ 2 rad,假定没有预旋
                              
                               =((π* D1 ^ 2) / 4) * 4;
                              
                               Ai =((π* D1 ^ 2) / 4);%截面面积叶轮(m ^ 2)
                              
                               广告=((π* D1 ^ 2) / 4);%截面面积扩散器(m ^ 2)
                              
                               李= 0.053;%的意思是通道长度(米)
                              
                               ld = 0.053;%的意思是通道长度(米)
                              
                               beta1b = 0.61;%叶轮叶片进口角()
                              
                               Ch = 4 * 0.3164 *(重新)^ -0.25;%摩擦损失系数()
                              
                               kf = (Ch *李)/(2 *迪* ro1 ^ 2 * Ai ^ 2 * (sin (beta1b)) ^ 2) * 4;%摩擦损失叶轮
                              
                               kfd = (Ch * ld) / (2 * Dd * ro1广告^ ^ 2 * 2 * (sin (beta1b)) ^ 2) * 4;%摩擦损失扩散器
                              
                               σ= 0.9;%滑移系数()
                              
                               D2 = 0.128;%叶轮直径提示(m)
                              
                               alfa2b =每股((D1 * tan (beta1b)) /(σ* D2));
                              
                               m = [-0.2:0.001:0.8];%(公斤/ s)质量流量
                              
                               cn = 10;
                              
                               % %的速度
                              
                               giri_medium = 55000;%节流阀孔径
                              
                               delta_giri = 400;% gamma_T_max e gamma_T_min之间的范围
                              
                               giri_max = giri_medium + (delta_giri / 2);%最大光圈
                              
                               giri_min = giri_medium - (delta_giri / 2);%最小孔径
                              
                               A_rpm = (giri_max - giri_min) / 2;%振幅
                              
                               b_rpm = (giri_max + giri_min) / 2;
                              
                               rpm_frequency = 20;%引擎的频率(赫兹)
                              
                               w_rpm = rpm_frequency * 2 *π;%引擎的频率(弧度/秒)
                              
                               % time_rpm = [25:0.01:40];%原
                              
                               time_rpm = [25:1:40];% MN
                              
                               giri = A_rpm * sin (w_rpm * time_rpm) + b_rpm;
                              
                               wh = a01。* (Ac. /副总裁。* Lc)。^ 0.5;% Helmotz的频率
                              
                               % %效率损失
                              
                               deltanbf = 0.02;%回流损失(0.02)环形扩散- 0,03介质体积- 0.05叶片扩压器
                              
                               deltanv = 0.035;%蜗壳损失()
                              
                               deltanc = 0.00;%间隙损失()
                              
                               deltand = 0.00;%扩散损失()
                              
                               德尔塔= deltanbf + deltanv + deltanc + deltand;
                              
                               % %节流阀的参数
                              
                               gamma_T =。;%节流阀孔径
                              
                               Delta_gamma_T = .02点;% gamma_T_max e gamma_T_min之间的范围
                              
                               gamma_T_max = gamma_T + (Delta_gamma_T / 2);%最大光圈
                              
                               gamma_T_min = gamma_T - (Delta_gamma_T / 2);%最小孔径
                              
                               = (gamma_T_max - gamma_T_min) / 2;%振幅
                              
                               b = (gamma_T_max + gamma_T_min) / 2;
                              
                               % %压缩机的频率(频率)
                              
                               compressor_frequency = 10;%引擎的频率(赫兹)
                              
                               w = compressor_frequency * 2 *π;%引擎的频率(弧度/秒)
                              
                               %保存main_parameters.mat
                              
                               % %颂歌的参数
                              
                               t0 (1) = 0;%仿真开始[s]
                              
                               (1)= time_rpm特遣部队(2);%模拟完成的[s]
                              
                               日元(1)= 0;%初始条件1
                              
                               y2 (1) = 0;%初始条件2
                              
                               init = (y1 y2) ';%初始条件向量
                              
                               sample_frequency = compressor_frequency * 8;%采样频率(赫兹)(至少选择一个频率
                              
                               % 10 - 12倍的分析
                              
                               %定义好。只是快速尝试2 - 3倍)
                              
                               maximum_time_step = 1 / sample_frequency;
                              
                               选择= odeset (“MaxStep”,maximum_time_step);%最大时间步大小
                              
                               保存main_parameters.mat
                              
                               t_all = [];% MN
                              
                               y_all = [];% MN
                              
                               为i = 1:(长度(giri) 2)
                              
                               U1 (i) = (D1 *π。*义理(我))/ 60;
                              
                               U2 (i) = (D2 *π。*义理(我))/ 60;
                              
                               B (i) = U1 (i)。/ (2。* wh。* Lc);
                              
                               [t、y] = ode113 (@greitzer (t0(我),特遣部队(i)], [y1 (i), y2 (i)],选项,U1 (i), U2 (i), B (i));%我发现ode113方式更有效率
                              
                               t0 (i + 1) = time_rpm (i + 1);%仿真开始[s]
                              
                               特遣部队(i + 1) = time_rpm (+ 2);%模拟完成的[s]
                              
                               日元(i + 1) = y(结束,1);%初始条件1
                              
                               y2 (i + 1) = y (, 2);%初始条件2
                              
                               %。内格罗蓬特:concat t y为所有迭代
                              
                               t_all = [t_all;t];% MN
                              
                               y_all = [y_all;y);% MN
                              
                               % array_val = 0(长度(giri) 2);
                              
                               %的阴谋(t y (: 1))
                              
                               %网格
                              
                               %网格小
                              
                               %包含(“t [s]”)
                              
                               % ylabel(“\φ”)
                              
                               %抓住
                              
                               图(%)
                              
                               %的阴谋(t, y (:, 2))
                              
                               %网格
                              
                               %网格小
                              
                               %包含(“t [s]”)
                              
                               % ylabel(“\ψ”)
                              
                               %抓住
                              
                               % savefig (10 hz)
                              
                               结束
                              
                               情节(t0, y1)
                              
                               网格在
                              
                               网格小
                              
                               包含(“t [s]”)
                              
                               ylabel (‘\φ)
                              
                               %抓住
                              
                               图()
                              
                               情节(t0, y2)
                              
                               网格在
                              
                               网格小
                              
                               包含(“t [s]”)
                              
                               ylabel (“\ψ”)
                              
                               %抓住
                              
                               图()
                              
                               情节(t_all y_all)
                              
                               网格在
                              
                               网格小
                              
                               包含(“t [s]”)
                              
                               ylabel (“y”)

保罗•罗杰斯 2021年1月11日

noope,贴一个新的答案,而不是第一个留言,或者我不会接受按钮

登录置评。

在回答这个问题。

答案(1)

瑞安Fedeli 2021年1月8日

1
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/711948-storing-variables-after-a-loop-easy answer_593778

这个例子是简化,但为了节省值从每个循环迭代,所有你需要做的就是将其保存到一个数组中。星号的线条是重要的部分:

                             z = linspace(1、50岁,10);%任意的向量
                            
                             b = z;例如函数%任意值
                            
                             array_val = 0(1、长度(z));% preallocate数组循环速度* * *
                            
                             为i = 1:长度(z)
                            
                             a = b(我)* 1.1;%示例函数来创建数据
                            
                             array_val (i) =一个;%保存价值从这个迭代到一个数组中* * *
                            
                             结束
                            
                             disp (array_val)

但如果“a”在这个例子不是一个单一的值,而是一个数组本身,需要将其保存到一个矩阵。在这种情况下,它会看起来像:

                             array_val = 0(长度(z));% 1)preallocate *矩阵* 0,长度(z)的长度(z)的大小。
                            
                             % 2)你需要找出这些维度应该是什么。
                            
                             % 3)这条线是有用的循环运行得更快,但这不是必要的
                            
                             为i = 1:长度(z)
                            
                             a = b (i)。* z;%的例子函数来创建一个行向量的数据(的数据)
                            
                             array_val (:, i) =一个;%保存数据的行向量矩阵的第i个位置
                            
                             结束
                            
                             disp (array_val)

array_val现在是一个矩阵,每一行数据的i迭代循环

我希望这有助于

1评论
显示隐藏 None

保罗•罗杰斯 2021年1月8日

我想跟随你的脑袋,但在你的情况下你也知道会有相同的长度z, infacct我得到一个下标赋值尺寸不匹配。

在我的例子中,我不知道多少分ode发出使用步骤。

登录置评。

在回答这个问题。

类别

万博1manbetx 应用程序通用应用程序

s manbetx 845

MATLAB

释放

R2014b

社区寻宝

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

开始狩猎!

存储变量后循环(容易)

8的评论
显示隐藏 7年长的评论

答案(1)

1评论
显示隐藏 None

另请参阅

类别

标签

s manbetx 845

释放

社区寻宝

存储变量后循环(容易)

8的评论 显示隐藏 7年长的评论

答案(1)

1评论 显示隐藏 None

另请参阅

类别

标签

s manbetx 845

释放

社区寻宝

8的评论
显示隐藏 7年长的评论

1评论
显示隐藏 None