如何找到两条曲线的交点(2个交点)- MATLAB答案- MATLAB中心 - 金宝app,下载188bet金宝搏,金宝搏官方网站

明星黾 2021年2月3日

0
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/735247-how-to-find-the-intersection-of-two-curves-2-points-of-intersection#answer_614027

我不能运行你的代码，因为我没有函数。

然而，作为一般规则，找到两条曲线的交点是很简单的，即使计算它们的函数不存在。

用你的向量试试这个:

                        X = linspace(0,10);%的独立变量
                       
                        y = 10 - (x-5) ^2;%因变量#1
                       
                        y2 = 5 * sin(2 *π* x / 5);%因变量#1
                       
                        曲线= y1 - y2;%减去曲线
                       
                        intidx =找到(diff(签署(curvdif)));近似过零指数
                       
                        为k = 1:元素个数(intidx)
                       
                        idxrng =[马克斯(intidx (k) 1): min(元素个数(x) intidx (k) + 1)];%索引范围内插
                       
                        Xv (k) = interp1(curdif (idxrng)， x(idxrng)， 0);%插值估计“x”在交叉点
                       
                        Yv (k) = interp1(x(idxrng)，y1(idxrng)， xv(k));%插值估计“y”在交叉点
                       
                        结束
                       
                        数字
                       
                        图(x, y1, x, y2)
                       
                        持有在
                       
                        图(十五、青年志愿“pg”，“MarkerSize”15)
                       
                        持有从
                       
                        网格
                       
                        传奇(“y_1”，“y_2”，“十字路口”，“位置”，“年代”）

我有许多关于飞行的快乐回忆 C172 年代!玩得开心!

4评论
显示隐藏 3年长的评论

里安农艾略特 2021年2月3日

这就跟你问声好！我不太明白这里发生了什么。我对Matlab比较陌生，不幸的是，由于Covid!在寻找交叉点的过程中到底发生了什么?

我不确定这是否能给我的代码提供更多的上下文，但我已经包含了缺失的函数，以及图的样子:

                             函数[P,rho,T,a] = tms (H, tempDiff, lapseDiff)
                            
                             g = 9.8065;
                            
                             R = 287.05287;
                            
                             γ= 1.4;
                            
                             lNaught = -0.0065 + lapseDiff;
                            
                             taught = 288.15 + tempDiff;
                            
                             tEleven = 216.65 + tempDiff;
                            
                             pNaught = 101325;
                            
                             T = tNaught + (lNaught * H);
                            
                             pEleven = pp0 *(1 + ((lNaught*11000)/ taught))^(-g/(R*lNaught)));
                            
                             pp0 *(1 + ((lNaught*H)/ taught))^(-g/(R*lNaught)));
                            
                             pressureLStrat = pEleven * exp (1) ^ ((g / (R * tEleven)) * (h - 11000));
                            
                             densTrop = pressureTrop / (R (T);
                            
                             densLStrat = pressureLStrat / (R * tEleven);
                            
                             如果(H < 11000)%对流层
                            
                             ρ= densTrop;
                            
                             P = pressureTrop;
                            
                             T = T;
                            
                             elseifH > = 11000%平流层
                            
                             ρ= densLStrat;
                            
                             P = pressureLStrat;
                            
                             T = tEleven;
                            
                             结束
                            
                             一个= sqrt(γ* R * T);
                            
                             结束
                            
                             函数[max_power_hp,max_power_W] = getPAmax (H)
                            
                             hp_sealevel = 180;%从飞行装备的源代码中获得的海平面时的最大能量[hp]
                            
                             [~，rho_0，~，~] = atmosProperties (0,0,0);海平面密度% kgm^-3
                            
                             [~，rho_H，~，~] = atmosProperties (H,0,0);在给定高度的密度% [kgm^-3]
                            
                             prop_eff = 0.85;%螺旋桨效率
                            
                             max_power_hp = prop_eff * hp_seallevel *(1.132*(rho_H/rho_0) - 0.132);
                            
                             max_power_W = max_power_hp * 746;
                            
                             结束

明星黾 2021年2月3日

我的代码需要两个向量，它们是同一个自变量的函数(或者可以计算为)。在那之后，简单地减去他们和使用我的代码找到交叉点。

我修改了我发布的代码，以解决您的问题，并添加或更改了这部分代码:

                             bclr = [c;(大小(c))]。' \ b (:);线性回归系数
                            
                             bvec = [V;(大小(V))]。‘* bclr;线性回归拟合
                            
                             = bvec - P_req(:);%减去曲线
                            
                             curvdif (isnan (curvdif)) = 1 e4;
                            
                             intidx =找到(diff(签署(curvdif)));近似过零指数
                            
                             为k = 1:元素个数(intidx)
                            
                             idxrng =[马克斯(intidx (k) 1): min(元素个数(V), intidx (k) + 1))%索引范围内插
                            
                             xv(k) = interp1(curdif (idxrng)， V(idxrng))， 0，“pchip”，“extrap”）;%插值估计“V”在交叉点
                            
                             yv(k) = interp1(V(idxrng)，P_req(idxrng)， xv(k)); / /计算结果%插值以估计“P_req”在交叉口
                            
                             结束
                            
                             流(1,“\ tIntersections: \ n”）
                            
                             流(1,'\t\tV = %9.4f, P_req = %9.1f\n'[十五(2:结束);青年志愿(2:结束)))
                            
                             图的目标是找到可用功率曲线(平线)和所需功率曲线(曲线)的交点
                            
                             %图例将在稍后添加
                            
                             情节(V, P_req);
                            
                             持有在
                            
                             情节(c, b);
                            
                             图(十五、青年志愿“处方”，“MarkerSize”15)
                            
                             包含(“助教(女士^{1})”）;
                            
                             ylabel (“权力reequired (W)”）
                            
                             传奇(“P_{点播}’，“b”，“十字路口”，“位置”，“西北”）
                            
                             ylim (5 e5 [0])%缩放显示交叉点(删除)

所以完整的代码现在是:

                             确定在给定高度的最大和最小速度的代码
                            
                             %的塞斯纳172便士
                            
                             H = 0;最终，这将成为一个函数的输入
                            
                             %用户只需要调用函数，就可以得到给定高度的最小和最大速度
                            
                             %飞机常数
                            
                             S = 16.16499443;%机翼面积(m ^ 2)
                            
                             g = 9.81;重力加速度% [ms^-2]
                            
                             W = 887.0748299 * g;飞机重量[kgf]
                            
                             C_D0 = 0.026;升力阻力系数为零
                            
                             K = 0.121;%阻力恒定
                            
                             我们希望研究的速度范围[马赫]
                            
                             M = linspace (0, 0.8,100);
                            
                             (ρ,~ ~]= atmosProperties (H, 0, 0);%在指定高度的声速和密度[kgm^-3] [ms^-1]
                            
                             [~, maxPA] = getPAmax (H);%给定高度的最大可用功率[W]
                            
                             V = m *;%速度[^ 1]女士
                            
                             C_L =(2. /ρ)。* (w / S)。* (1. / (v . ^ 2));%给定高度的升力系数
                            
                             C_d = c_d0 + k * c_l ^2;%给定高度阻力系数
                            
                             P_req = 0.5 .* rho .* v ^3 .* S .* C_D;%所需功率曲线
                            
                             b = [maxPA, maxPA];%B和c用于绘制一条线，该线对应于可用的最大功率
                            
                             c = [0272 .235190420871];
                            
                             bclr = [c;(大小(c))]。' \ b (:);线性回归系数
                            
                             bvec = [V;(大小(V))]。‘* bclr;线性回归拟合
                            
                             = bvec - P_req(:);%减去曲线
                            
                             curvdif (isnan (curvdif)) = 1 e4;
                            
                             intidx =找到(diff(签署(curvdif)));近似过零指数
                            
                             为k = 1:元素个数(intidx)
                            
                             idxrng =[马克斯(intidx (k) 1): min(元素个数(V), intidx (k) + 1))%索引范围内插
                            
                             xv(k) = interp1(curdif (idxrng)， V(idxrng)， 0);%插值估计“V”在交叉点
                            
                             yv(k) = interp1(V(idxrng)，P_req(idxrng)， xv(k)); / /计算结果%插值以估计“P_req”在交叉口
                            
                             结束
                            
                             流(1,“\ tIntersections: \ n”）
                            
                             流(1,'\t\tV = %9.4f, P_req = %9.1f\n'[十五(2:结束);青年志愿(2:结束)))
                            
                             图的目标是找到可用功率曲线(平线)和所需功率曲线(曲线)的交点
                            
                             %图例将在稍后添加
                            
                             情节(V, P_req);
                            
                             持有在
                            
                             情节(c, b);
                            
                             图(十五、青年志愿“处方”，“MarkerSize”15)
                            
                             包含(“助教(女士^{1})”）;
                            
                             ylabel (“权力reequired (W)”）
                            
                             传奇(“P_{点播}’，“b”，“十字路口”，“位置”，“西北”）
                            
                             ylim ([0 2 e5)%缩放显示交叉点(删除)

制作有交叉点的地块，并且:

                             十字路口:
                            
                             V = 8.1607, P_req = 114138.0
                            
                             V = 73.3073, P_req = 114138.0

做出适当的改变以得到你想要的结果。

．

明星黾 2021年2月4日

一如既往，我的荣幸!

玩得开心 C172 ！

登录评论。

你现在正关注这个问题

如何求两条曲线的交点(两条交点)

0评论
显示隐藏1以上的评论

接受的答案

4评论
显示隐藏 3年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

你现在正关注这个问题

如何求两条曲线的交点(两条交点)

0评论 显示隐藏1以上的评论

接受的答案

4评论 显示隐藏 3年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

一个错误发生

0评论
显示隐藏1以上的评论

4评论
显示隐藏 3年长的评论