如何改变非标量变量标量变量?gydF4y2Ba

53视图(30天)gydF4y2Ba

光gydF4y2Ba 2023年4月21日gydF4y2Ba

0gydF4y2Ba
链接gydF4y2Ba

这个问题直接联系gydF4y2Ba

//www.tatmou.com/matlabcentral/answers/1951143-how-to-change-non-scalar-variable-to-scalar-variablegydF4y2Ba

评论道:gydF4y2Ba 沃尔特·罗伯森gydF4y2Ba 2023年5月19日在41gydF4y2Ba

PRGA2001阻力系数(2). pdfgydF4y2Ba

美好的一天,我试着在matlab以下步骤解决问题给我们。我找到一首颂歌的积分0的下限和上限是x, x有多个值。年代的解决方案或者(距离),也将多个值。我试图纠正脚本但我得到一个错误,因为x限制不是标量。这是8和im的步骤3最终试图完成步骤8。gydF4y2Ba

这是我迄今为止感谢帮助我收到:gydF4y2Ba

                          % k =空气阻力系数gydF4y2Ba
                         
                          % (theta) = initial angle of elevation of the ball when it was kicked
                         
                          % x和y是水平和垂直空间坐标gydF4y2Ba
                         
                          % t是时间测量从球踢gydF4y2Ba
                         
                          v % | | = mag_v =是速度的大小(速度,随着时间的变化)gydF4y2Ba
                         
                          全球gydF4y2BaggydF4y2Ba
                         
                          m = 0.45gydF4y2Ba%的公斤gydF4y2Ba
                         
                          g = 9.81gydF4y2Ba% m / s ^ 2中重力加速度gydF4y2Ba
                         
                          v0 = 108gydF4y2Ba%在km / hgydF4y2Ba
                         
                          normv =√(v0 * cos(θ)。^ 2)+ (v0 * sin(θ)。^ 2))gydF4y2Ba
                         
                          t = 0: 0.1: 2;gydF4y2Ba
                         
                          θ= 30 *π/ 180;gydF4y2Ba%初始猜测θgydF4y2Ba
                         
                          k = 3gydF4y2Ba%初始猜测N的空气阻力gydF4y2Ba
                         
                          ICs = [0 0 v0 * cos(θ)v0 * sin(θ)]gydF4y2Ba
                         
                          % v = 10gydF4y2Ba
                         
                          % dx / dt = x1gydF4y2Ba
                         
                          % dx1 / dt = x2gydF4y2Ba
                         
                          % (v = (- k) | | (x1)gydF4y2Ba
                         
                          % dy / dt =日元gydF4y2Ba
                         
                          %₁/ dt y2 =gydF4y2Ba
                         
                          % my2 =(毫克)v (- k) | |(日元)gydF4y2Ba
                         
                          % (= @ (t, x) [x (1);x (2);- k * mag_v * x (1)]gydF4y2Ba
                         
                          % my2 = @ (t、y) [y (1);y (2);g - k - m * * mag_v * y (2))gydF4y2Ba
                         
                          % y (3) = dx / dt, y (4) = dy / dtgydF4y2Ba
                         
                          [T、Y] =数值(@ (T, Y)乐趣(T、Y、g、k、m), T, ICs);gydF4y2Ba
                         
                          %的阴谋(T Y (: 1:2))gydF4y2Ba
                         
                          Y (: 1)gydF4y2Ba% X (t)gydF4y2Ba
                         
                          Y (: 2)gydF4y2Ba% Y (t)gydF4y2Ba
                         
                          x = [Y (: 1))”gydF4y2Ba
                         
                          y = [y (:, 2)] 'gydF4y2Ba
                         
                          interp1 (Y (: 1), Y (:, 2), 0.7)gydF4y2Ba
                         
                          f = @ (x) sqrt (1 + (y / y (3) (4))。^ 2)gydF4y2Ba
                         
                          四(f, 0, x)gydF4y2Ba
                         
                          % s =gydF4y2Ba
                         
                          函数gydF4y2Bady =乐趣(t、y、g、k、m)gydF4y2Ba
                         
                          mag_v =√y (3) (4) ^ 2 + y ^ 2);gydF4y2Ba
                         
                          dy = 0 (4,1);gydF4y2Ba
                         
                          dy (1) = y (3);gydF4y2Ba
                         
                          dy (2) = y (4);gydF4y2Ba
                         
                          dy (3) = - k * mag_v * y (3) / m;gydF4y2Ba
                         
                          dy (4) = (- m * g-k * mag_v * y (4)) / m;gydF4y2Ba
                         
                          结束gydF4y2Ba

12个评论gydF4y2Ba
显示11年长的评论gydF4y2Ba隐藏11年长的评论gydF4y2Ba

光gydF4y2Ba 2023年4月24日gydF4y2Ba

编辑:gydF4y2Ba沃尔特·罗伯森gydF4y2Ba 2023年5月19日在38gydF4y2Ba

嗨我有这段代码,之后随着问题表附在第一个评论,我想创建一个while循环的最小二乘法。任何方式可以使用while循环?gydF4y2Ba

                                % compute_drag_scriptgydF4y2Ba
                               
                                D = [0;0.5;1.1;1.7;2.3;2.9;3.5;4.0;4.6;5.2;gydF4y2Ba…gydF4y2Ba
                               
                                5.7;6.3;6.8;7.4;7.9;8.4;9.0;9.5;10.0;10.6;gydF4y2Ba…gydF4y2Ba
                               
                                11.1;11.6;12.1;12.6;13.1;13.6;14.1;14.6;15.1;15.6;gydF4y2Ba…gydF4y2Ba
                               
                                16.1;16.5;17.0;17.5;18.0;18.4;18.9;19.4);gydF4y2Ba
                               
                                V = [108;107;106;105;104;104;103;102;101;100;gydF4y2Ba…gydF4y2Ba
                               
                                99;99;98;97;97;96;95;94;94;93;gydF4y2Ba…gydF4y2Ba
                               
                                92;92;91;91;90;89;89;88;88;87;gydF4y2Ba…gydF4y2Ba
                               
                                87;86;86;85;85;84;83;82);gydF4y2Ba
                               
                                interp_speeds = compute_drag (D, V, 0.0095, 30);gydF4y2Ba
                               
                                函数gydF4y2Bainterp_speeds = compute_drag (gydF4y2Ba…gydF4y2Ba
                               
                                true_distances true_speeds k theta_degrees)gydF4y2Ba
                               
                                参数gydF4y2Ba
                               
                                true_distancesgydF4y2Ba(1):{mustBeNonnegative}gydF4y2Ba
                               
                                true_speedsgydF4y2Ba(1):{mustBeNonnegative}gydF4y2Ba
                               
                                kgydF4y2Ba(1,1){mustBeNonnegative}gydF4y2Ba
                               
                                % drag_coeffgydF4y2Ba
                               
                                theta_degreesgydF4y2Ba(1,1){mustBeInRange (theta_degrees, 0180)}gydF4y2Ba
                               
                                结束gydF4y2Ba
                               
                                % %输入验证gydF4y2Ba
                               
                                如果gydF4y2Ba长度(true_distances) ~ = (true_speeds)gydF4y2Ba
                               
                                错误(gydF4y2Ba的距离和速度向量长度必须相同。gydF4y2Ba);gydF4y2Ba
                               
                                结束gydF4y2Ba
                               
                                如果gydF4y2Ba任何(diff (true_distances) < 0)gydF4y2Ba
                               
                                错误(gydF4y2Ba“距离不应该下去!”gydF4y2Ba);gydF4y2Ba
                               
                                结束gydF4y2Ba
                               
                                % %常数和参数gydF4y2Ba
                               
                                g = 9.81;gydF4y2Ba% m /秒方gydF4y2Ba
                               
                                m = 0.45;gydF4y2Ba%公斤gydF4y2Ba
                               
                                v0_kph = 108;gydF4y2Ba%千米每小时,将转换为米每秒gydF4y2Ba
                               
                                k = 0.00000001;gydF4y2Ba
                               
                                ρ= 1.293gydF4y2Ba%公斤/ (m) ^ 3gydF4y2Ba
                               
                                d = 0.7 /πgydF4y2Ba% mgydF4y2Ba
                               
                                % %转换gydF4y2Ba
                               
                                v0 = v0_kph * 1000/3600;gydF4y2Ba%公里米每秒gydF4y2Ba
                               
                                %我们也可以theta_guess转换为弧度,但我们没有理由gydF4y2Ba
                               
                                %,Matlab以来程度变异的三角函数。gydF4y2Ba
                               
                                % %初始条件gydF4y2Ba
                               
                                x0 = 0;gydF4y2Ba
                               
                                y0 = 0;gydF4y2Ba
                               
                                dxdt0 = v0 * cosd (theta_degrees);gydF4y2Ba
                               
                                dydt0 = v0 *信德(theta_degrees);gydF4y2Ba
                               
                                init_vals = [x0;y0;dxdt0;dydt0];gydF4y2Ba
                               
                                % %步骤2gydF4y2Ba
                               
                                %设置diffeq将数值:gydF4y2Ba
                               
                                %瓦尔(1)x: x”(即瓦尔(3))gydF4y2Ba
                               
                                %瓦尔(2)y: y '(即瓦尔(4))gydF4y2Ba
                               
                                %瓦尔(3)x”: (v - k | | / m) * x 'gydF4y2Ba
                               
                                %瓦尔(4)y ' - g - v (k | | / m) * y 'gydF4y2Ba
                               
                                all_diffeq = @ (t, VALS) (gydF4y2Ba
                               
                                瓦尔斯(3);gydF4y2Ba
                               
                                瓦尔斯(4);gydF4y2Ba
                               
                                瓦尔斯- k *√(3) ^ 2 + VALS(4) ^ 2) *瓦尔(3)/ m;gydF4y2Ba
                               
                                瓦尔斯- g - k *√(3) ^ 2 + VALS (4) ^ 2) * VALS (4) / m);gydF4y2Ba
                               
                                ];gydF4y2Ba
                               
                                %得到一个合适的时间范围,我用最坏的情况:gydF4y2Ba
                               
                                %的速度最慢(例如82公里)球的距离gydF4y2Ba
                               
                                %(例如19.4米)。gydF4y2Ba
                               
                                minspeed = 1000/3600 min (true_speeds) *;gydF4y2Ba
                               
                                max_necessary_time = (true_distances(结束))/ minspeed;gydF4y2Ba
                               
                                tspan = [0, max_necessary_time];gydF4y2Ba
                               
                                [t, vals] =数值(all_diffeq, tspan init_vals);gydF4y2Ba
                               
                                x = vals (: 1);gydF4y2Ba
                               
                                y = vals (:, 2);gydF4y2Ba
                               
                                dxdt = vals (:, 3);gydF4y2Ba
                               
                                dydt = vals (:, 4);gydF4y2Ba
                               
                                % s =√x。^ 2 + y ^ 2);gydF4y2Ba
                               
                                % %步骤3gydF4y2Ba
                               
                                s_segments =√1 + (dydt. / dxdt)。^ 2);gydF4y2Ba
                               
                                s_step = [0; diff (x)];gydF4y2Ba
                               
                                s = cumsum (s_segments。* s_step);gydF4y2Ba
                               
                                v =√dxdt。^ 2 + dydt。^ 2);gydF4y2Ba
                               
                                v_kph = v * 3600/1000;gydF4y2Ba
                               
                                % %步骤4gydF4y2Ba
                               
                                interp_speeds = interp1 (s v_kph true_distances);gydF4y2Ba
                               
                                如果gydF4y2Ba任何(isnan (interp_speeds))gydF4y2Ba
                               
                                警告(gydF4y2Ba“南。你应该增加上限时间。”gydF4y2Ba);gydF4y2Ba
                               
                                结束gydF4y2Ba
                               
                                结束gydF4y2Ba
                               
                                %步骤5gydF4y2Ba
                               
                                而gydF4y2BaV_interpolated > = 0gydF4y2Ba
                               
                                E = cumsum (V_interpolated - V_var) ^ 2gydF4y2Ba
                               
                                结束gydF4y2Ba

严厉的gydF4y2Ba 2023年5月19日在六点四十gydF4y2Ba

你好,gydF4y2Ba

是的,一个最小二乘优化可以实现使用while循环。gydF4y2Ba

在以下几点,我解释的步骤可以利用来实现我们的目标。gydF4y2Ba

设置while循环停止条件:gydF4y2Ba我建议设置while循环停止条件根据误差公差值而在while循环一个零值。下面提供了相同的示例代码:gydF4y2Ba

托尔= 0.00001;gydF4y2Ba%定义一个公差值错误停止优化gydF4y2Ba

而gydF4y2Ba(dE_dtheta > tol | | dE_dk > tol)gydF4y2Ba

%变量实现算法和修正。gydF4y2Ba

结束gydF4y2Ba

实现最小二乘优化算法:gydF4y2Ba附加的PDF文档中指定的算法已经在“最小二乘法”标题下。如果你面对具体问题的实现,请张贴在这里。gydF4y2Ba
纠正变量:gydF4y2Ba修正在初始猜测变量(“k”和“θ”在这种情况下)与每个”,而“循环迭代的优化是一个重要的步骤。您可以使用以下方法执行相同的梯度下降:gydF4y2Ba

α= 0.001;gydF4y2Ba%定义了学习速率或步长梯度下降法。gydF4y2Ba

%在当前迭代计算误差值后,更新如下给出优化变量。gydF4y2Ba

θ=θ-α* dE_dtheta;gydF4y2Ba

k = k -α* dE_dk;gydF4y2Ba

沃尔特·罗伯森gydF4y2Ba 2023年5月19日在41gydF4y2Ba

                                而gydF4y2BaV_interpolated > = 0gydF4y2Ba
                               
                                E = cumsum (V_interpolated - V_var) ^ 2gydF4y2Ba
                               
                                结束gydF4y2Ba

你的while循环gydF4y2Ba 必须gydF4y2Ba 更新至少一个变量被测试的条件,否则你有无限循环。gydF4y2Ba

例外:它是有效循环类似gydF4y2Ba

                                而gydF4y2Ba真正的gydF4y2Ba
                               
                                做gydF4y2Ba一些计算gydF4y2Ba
                               
                                如果gydF4y2Ba条件是gydF4y2Ba见过gydF4y2Ba
                               
                                打破gydF4y2Ba;gydF4y2Ba
                               
                                结束gydF4y2Ba
                               
                                结束gydF4y2Ba