错误使用内联函数和fzero

1视图(30天)

显示旧的评论

马纳尔M Deshpande 2014年3月20日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/122458-error-in-using-inline-function-and-fzero

评论道: 马纳尔M Deshpande2014年3月23日

答:接受明星黾

你好,有人能帮我吗?我想找出我做程序的错误。唯一的变量X和Cp是一个输入。

                         Cp =双(0 (49792 1));
                        
                         P =大小(x.Time_Series);
                        
                         Q = P (1,1);
                        
                         Cp (1, - 1) = 0.126826003999670;
                        
                         Cp (2, 1) = 0.325943036562508;
                        
                         为k = 3:问;
                        
                         k
                        
                         hx =内联(' 12.175 * (X - 2 * Cp (k - 1, - 1) + Cp (k-2 1)) + 13.525 * ((X - Cp (k - 1, - 1)) * abs (X - Cp (k - 1, - 1))) + X - SINGLE_FINAL_TIME_SERIES_OF_TAPS_NEAR_DOMINANT_OPENING (k, 1)”,“X”,Cp);
                        
                         惠普= fzero (hx 0 Cp SINGLE_FINAL_TIME_SERIES_OF_TAPS_NEAR_DOMINANT_OPENING (k, 1));
                        
                         Cp (k) =惠普(1);
                        
                         结束
                        
                         我我得到这个错误
                        
                         错误:表达式左边的等号不是一个有效的目标为一个
                        
                         任务。

如果在内联函数,我删除Cp后给出X(在最后支架附近),那么它显示错误:太多的输入内联函数。请帮助。谢谢。

接受的答案

明星黾 2014年3月20日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/122458-error-in-using-inline-function-and-fzero answer_129411

它不像 “= =” 在 SINGLE_FINAL_TIME_SERIES_OF_TAPS_NEAR_DOMINANT_OPENING (k, 1) = = 0 。

与使用工作空间变量的匿名函数, 内联函数忽略任何不在他们的参数列表。不知道 SINGLE_FINAL_TIME_SERIES_OF_TAPS_NEAR_DOMINANT_OPENING (k, 1) 或 k 是这样的。

同时, fzero 只接受变量函数。

9日评论
显示8年长的评论隐藏8年长的评论

马纳尔M Deshpande 2014年3月21日

编辑:马纳尔M Deshpande 2014年3月21日

而不是使用内联fzero或匿名函数,我用“解决”所需的解决方案。万博尤文图斯我的代码如下所示。这里的问题是这个项目变得非常缓慢。这是大约10小时解决得到49790的值。SINGLE_FINAL_TIME_SERIES_OF_TAPS_NEAR_DOMINANT_OPENING (k, 1)是49792 * 1的输入数组值。X的值(存储为VD(1))发现每次和这个值从SINGLE_FINAL_TIME_SERIES一起价值……成为下一个的输入方程。有什么方法可以使程序更快呢?我的代码:

                                  抽搐
                                 
                                  Cp (1, - 1) = 0.126931068442331;
                                 
                                  Cp (2, 1) = 0.326161603224261;
                                 
                                  信谊X;
                                 
                                  为k = 3:问;
                                 
                                  k
                                 
                                  VD解决(1)= (12.175 * (X - 2 * Cp (k - 1, - 1) + Cp (k-2 1)) + 13.525 * (X-Cp (k - 1, - 1)) * abs (X-Cp (k - 1, - 1)) + X - SINGLE_FINAL_TIME_SERIES_OF_TAPS_NEAR_DOMINANT_OPENING (k, 1) = = 0, X);
                                 
                                  Cp (k) = VD (1);
                                 
                                  结束
                                 
                                  toc

谢谢。

马纳尔M Deshpande 2014年3月21日

编辑:马纳尔M Deshpande 2014年3月21日

好吧,如果我还不清楚,那么我再解释给你。我第一次手动或使用解决2 *找到第一个2的Cp值。例如,当X是发现在第一步,然后它被用作Cp输入在第二时间步。我使用的值从Single_time ....当我输入。这种方式是重复直到最后。我认为你上面写的是对的。我跑我的程序通过使用解决约10000 Cp值但值是不正确的。以下是前20的价值观Single_time__ 1.88848000000000 1.74232000000000 1.54512000000000 1.46624000000000 1.20408000000000 0.874640000000000 0.570720000000000 1.13216000000000 1.31544000000000 1.43608000000000 1.63328000000000 1.42912000000000 1.22032000000000 1.49640000000000 1.63328000000000 1.71680000000000 1.51728000000000 0.965120000000000 0.468640000000000 0.740080000000000

明星黾 2014年3月22日

编辑:明星黾 2014年3月22日

生产你想要的结果吗?情节看起来很奇怪,但我不知道你在研究什么。(我画它,看看它看起来合理。)

                                  SINGLE_FINAL = [1.88848000000000 1.74232000000000 1.54512000000000 1.46624000000000 1.20408000000000 0.874640000000000 0.570720000000000 1.13216000000000 1.31544000000000 1.43608000000000 1.63328000000000 1.42912000000000 1.22032000000000 1.49640000000000 1.63328000000000 1.71680000000000 1.51728000000000 - 0.965120000000000 0.468640000000000 - 0.740080000000000);
                                 
                                  Cp (1) = 0.126826003999670;
                                 
                                  Cp (2) = 0.325943036562508;
                                 
                                  % hx =内联(12.175 * (X - 2 * Cp (k - 1, - 1) + Cp (k-2 1)) + 13.525 * ((X - Cp (k - 1, - 1)) * abs (X - Cp (k - 1, - 1))) + X - SINGLE_FINAL_TIME_SERIES_OF_TAPS_NEAR_DOMINANT_OPENING (k, 1)’,‘X’, Cp);
                                 
                                  %称“f1”:
                                 
                                  %定义cp = [cp (k - 1);Cp (k-2)]
                                 
                                  % f1 (x, cp,科幻(k))
                                 
                                  f1 = @ (X, cp, SFTSOTNDO) 12.175 . * (X - 2。* cp (1) + (2) cp) + 13.525 . * ((X - cp (1))。* abs (X - cp (1))) + X - SFTSOTNDO;
                                 
                                  为k = 3:长度(SINGLE_FINAL)
                                 
                                  cp = [cp (k - 1);Cp (k-2)];
                                 
                                  科幻小说= SINGLE_FINAL (k);
                                 
                                  X = fzero (@ (X) (f1 (X, cp,科幻小说),1);
                                 
                                  Cp (k) = X;
                                 
                                  结束
                                 
                                  图(1)
                                 
                                  情节(SINGLE_FINAL, Cp)
                                 
                                  网格
                                 
                                  包含(单决赛的)
                                 
                                  ylabel (Cp的)

我一直在你内联函数引用,但说出来。

我也vectorised f1 因为这通常可以避免的问题。

马纳尔M Deshpande 2014年3月23日

尤里卡!它产生的结果非常快。对所有49792个值难以置信的163秒。非常感谢。我想解释一下你这里我试图实现什么尽管它可能不是与你的领域相关。Cp意味着压力系数。Pr.系数是衡量在边界层风洞100年代在500赫兹的频率。这就是为什么50000年价值大约由SFTSOTNDO表示。测量的值是外部Cp和我想获得内部Cp使用上面的方程。没有必要这两个之间的阴谋。 However, the values I have got are relatively higher. They should all be less than 3 or 4. But leaving that, it is major work done for me. Thanks a lot again. I have named this program Star_Strider as a tribute and a sign of respect.