matlab运行答案总是ans=3，请指点原因。在线等

佳琳

请问一下，我的matlab程序运行答案总是ans=3，是怎么回事，请高手指点。
我已经清空过好几次了，都是一个结果。

1. a=linspace(0,0,1000);
   
2. x=zeros(1000,100);
   
3. y=zeros(1000,100);
   
4. w=[];
   
5. m=[];
   
6. k=0;
   
7. i=0;
   
8. j=0;
   
9. t=0;
   
10. d=0;
    
11. s=0;
    
12. g1=0;
    
13. g2=0;
    
14. wlong=0;
    
15. w(1)=0;
    
16. step=0;
    
17. while step<2
    
18. 
    
19. %随机选取顾客
    
20. while t<=600
    
21. kk=round(rand()*999)+1;
    
22. if aa(kk)==0
    
23. aa(kk)=1;
    
24. t=t+1;
    
25. end;
    
26. end;
    
27. k=round(rand()*400);
    
28. t=0;
    
29. while t < =k
    
30. kk = round(rand() * 999) + 1;
    
31. if aa(kk) == 0
    
32. aa(kk) = 2;
    
33. t = t + 1;
    
34. end;
    
35. end;
    
36. for i = 1:1000
    
37. if aa(i) == 0
    
38. aa(i) = 3;
    
39. end;
    
40. end;
    
41. %设置q，表示顾客i对j是否偏好
    
42. for i = 1:1000
    
43. for j = 1:100
    
44. if p(i, j)==0
    
45. q(i, j) = 0;
    
46. else
    
47. q(i,j)=1;
    
48. end;
    
49. end;
    
50. end;
    
51. 
    
52. %开始随即挑选DVD编号
    
53. while s < 100000000
    
54. %给第一组的顾客选DVD，存入X
    
55. for i = 1:1000
    
56. if aa(i) == 1 | aa(i) == 2
    
57. t = 1;
    
58. while t <=3
    
59. kk=round(rand() * 99)+1;
    
60. if x(i, kk) == 0
    
61. x(i, kk) = 1;
    
62. t = t + 1;
    
63. end;
    
64. end;
    
65. end;
    
66. end;
    
67. %给第二组的顾客选DVD，存入y
    
68. for i = 1:1000
    
69. if aa(i) == 3
    
70. t = 1;
    
71. while t <=3
    
72. kk=round(rand() * 99)+1;
    
73. if y(i, kk) == 0
    
74. y(i, kk) = 1;
    
75. t = t + 1;
    
76. end;
    
77. end;
    
78. end;
    
79. end;
    
80. %计算两组顾客分到了想要的DVD的数量
    
81. for r = 1:1000
    
82. for g = 1:100
    
83. if aa(r) == 1
    
84. g1 = g1 + x(r, g) * q(r, g) * y(r, g);
    
85. g2 = g2 + x(r, g) * q(r, g) * y(r, g);
    
86. elseif aa(r) == 2
    
87. g1 = g1 + x(r, g) * q(r, g);
    
88. elseif aa(r) == 3
    
89. g2 = g2 + y(r, g) * q(r, g);
    
90. end;
    
91. end;
    
92. end;
    
93. %看通过率是否达到95%，达到计算此时的满意度和DVD数
    
94. if g1 * (600 + k) / 1600 / (600 + k) + g2 * (1000 - k) / 1600 / (1000 - k) > 0.95
    
95. wlong = wlong + 1;
    
96. for r = 1:1000
    
97. for g = 1:100
    
98. w(wlong) = w(wlong) + x(r, g) * p(r, g);
    
99. m(wlong) = g1 + g2;
    
100. end;
     
101. end;
     
102. end;
     
103. 
     
104. s = s + 1;
     
105. end;
     
106. 
     
107. step=step+1
     
108. end;
     
109. 
     
110. 
     
111. w
     
112. m

复制代码

jiao320123

呵呵 大家好奇嘛 来观看下~~~~  

foreveringxq

a=4；
a也是3吗

厚积薄发

利用clear命令清空下变量空间，然后再试

cp苏薇儿

可能是你以前的变量存储在里面你不知道，可以打开workspace看一下，直接将这个变量删除

machuan1214

利用clear命令清空下变量空间

花道

我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。   

佳琳

回复 厚积薄发 的帖子

初学，可不可以给详细一点的说明呀

   

鱼宝1990

先清空工作区中的变量~~~~~~clear all

study18neu99

变量定义有点乱，没有发挥matlab矩阵运算的优势，建议清空所有变量，并检查循环部分代码