[求助]真诚请教：0.618法的一维搜索

tongjiayan

这是MATLAB中的优化问题，用0.618法做，需要编个程序：
求min y(x)=(x*x-1)(x*x-1) 在区间 [0,2] 上的一个解，要求精度为10^(-1)。
由于在下是第一次遇到这样的问题，所以不怎么懂，还望各位前辈不吝指教！
谢谢！

linda

1. #include
   
2. #include
   
3. /*黄金分割法求最小值的C++程序，部分变量及函数书写并不规范*/
   
4. 
   
5. //δ为题给精度
   
6. int n = (lnδ/ln0.618 + 1) + 1;
   
7. int i;
   
8. float f(float ai, float bi)
   
9. {
   
10. a(i + 1) = ai + 0.618(bi - ai);
    
11. return ai + 1;
    
12. }
    
13. 
    
14. float g(float ai, float bi)
    
15. {
    
16. b(i + 1) = ai + 0.382(bi - ai);
    
17. return b(i + 1);
    
18. }
    
19. 
    
20. float F(float ai, float bi)
    
21. {
    
22. //题给的f(x)函数式;
    
23. return ;
    
24. }
    
25. 
    
26. float A(float ai, float bi)
    
27. {
    
28. int i = 1;
    
29. float result;
    
30. L:do
    
31. {
    
32. a(i + 1) = f(float ai, float bi);
    
33. b(i + 1) = g(float ai, float bi);
    
34. float F1 = F(float ai, float bi);
    
35. float F2 = F(float a(i + 1), float b(i + 1));
    
36. ai = ai, bi = b(i + 1);
    
37. i ++;
    
38. 
    
39. }while(i <= n && F1 >= F2)
    
40. if(i < n)
    
41. {
    
42. B(float ai, float bi);
    
43. }
    
44. else
    
45. result = F2;
    
46. return result;
    
47. 
    
48. }
    
49. 
    
50. float B(float ai, float bi)
    
51. {
    
52. do
    
53. {
    
54. a(i + 1) = f(float ai, float bi);
    
55. b(i + 1) = g(float ai, float bi);
    
56. float F1 = F(float ai, float bi);
    
57. float F2 = F(float a(i + 1), float b(i + 1));
    
58. ai = a(i + 1), bi = bi;
    
59. i ++;
    
60. 
    
61. }while(i <= n && F1 <= F2)
    
62. if(i < n)
    
63. {
    
64. goto L;
    
65. }
    
66. else
    
67. result = F1;
    
68. return result;
    
69. 
    
70. }
    
71. 
    
72. void main()
    
73. {
    
74. int i = 1;
    
75. float A(float ai, float bi);
    
76. cout<<"最小值为："<}

复制代码

这是C语言写的，你改写一下就行了

绿荷雨丝

不用这么麻烦

matlab自带了函数的

用fminbnd就行了

tongjiayan

您好！
感谢您的解答，我会好好研究的。再请教一个问题，只在C++中编写完程序后再导入MATLAB吗？怎么操作比较好？因为我从来没有操作过类似的情况。
非常谢谢！

tongjiayan

您好！
请问您所指的命令该如何使用？这个命令本身就是用0.618法做的吗？还要不要添加别的程序？敢问您有没有现成的例子能让我学习一下的？
非常感谢！

frogfish

写一个fun1.m文件：
function f=fun1(x);
f=(x*x-1)*(x*x-1);

然后才命令窗口输入
[x,y]=fminbnd('fun1',0,2)  就行了

frogfish

二楼的不是这个意思
他的意思是讲C语言的程序翻译成matlab的代码
相当于是自己编程，只是有一个C语言程序作为参照

tongjiayan

您好！
谢谢您的解答，我会试一下的，若再有什么问题再来请教。
另外我想再请教一问：在优化中，是否各种探索法都对应有各自的命令，只要在使用时输入命令和参数就行？
非常感谢！You are my teacher!

tongjiayan

原来如此。看来我太孤陋寡闻了。

tongjiayan

您好！
对不起，又来打扰了，我刚才根据您的指导，把第一行的命令输入，打了回车后系统给了如下的提示：
function f=fun1(x);
??? Strings passed to EVAL cannot contain function declarations.

请教这句话是什么意思？
谢谢！

frogfish

请仔细看6楼，要先写一个.m文件
在命令窗口输入的是[x,y]=fminbnd('fun1',0,2)

frogfish

看你要用什么？用来做什么？

如果你真要用matlab建议先拿本matlab的书先看看，大致了解一下

很简单一般一两天就能初步掌握了

tongjiayan

您好！
谢谢！原来是这个意思，先新建一个M文件，然后再使用。我第一次玩这个东西，多多指教！
虽然您这样批评我感觉很受打击，但还是要谢谢您！

frogfish

没有批评的意思，只是给我的感觉你是刚刚用matlab
很多基本的东西书上都写得很详细了，所以建议你先看看书
有利于你使用matlab

zhj111111

xiexie