无法执行赋值，左侧索引与右侧大小不兼容

主程序中涉及到了RSQ代码：

1. % --- Executes on button press in electrodeCalibrate.
   
2. function electrodeCalibrate_Callback(hObject, eventdata, handles)
   
3. % hObject    handle to electrodeCalibrate (see GCBO)
   
4. % eventdata  reserved - to be defined in a future version of MATLAB
   
5. % handles    structure with handles and user data (see GUIDATA)
   
6. global pHmeter slope intercept Rsq efficeny
   
7. 
   
8. [slope, intercept, Rsq, efficeny] = calibration(pHmeter);
   
9. set(handles.slopeWindow,'String',num2str(slope));
   
10. set(handles.interceptWindow,'String',num2str(intercept));
    
11. set(handles.efficencyWindow,'String',num2str(efficeny));
    
12. 
    
13. %%save calibration and to file
    
14. calibdata{1}=num2str(slope);
    
15. calibdata{2}=num2str(intercept);
    
16. calibdata{3}=num2str(efficeny);
    
17. 
    
18. save('electrodeCalibData.mat','calibdata');
    

复制代码

其中RSQ.m的具体代码如下：

1. function [slope, intercept,xModel, yModel,Rsq] = RSQ(x,y)
   
2. 
   
3. if size(x,2)>size(x,1)
   
4. 
   
5.     x=x';
   
6. 
   
7. end
   
8. 
   
9. if size(y,2)>size(y,1)
   
10. 
    
11.     y=y';
    
12. 
    
13. end
    
14. 
    
15. workMatrix=zeros(size(x));
    
16. 
    
17. workMatrix(:,1)=x;
    
18. 
    
19. workMatrix(:,2)=y;
    
20. 
    
21. 
    
22. 
    
23. %remove NaNs
    
24. 
    
25. workMatrix = workMatrix(all(~isnan(workMatrix),2),:);
    
26. 
    
27. %remove infs
    
28. 
    
29. workMatrix = workMatrix(all(~isinf(workMatrix),2),:);
    
30. 
    
31. 
    
32. 
    
33. % calculate slope and intercept
    
34. 
    
35. fitParam=polyfit(workMatrix(:,1),workMatrix(:,2),1);
    
36. 
    
37. %calculate model x and y
    
38. 
    
39. 
    
40. 
    
41. slope=fitParam(1);
    
42. 
    
43. intercept=fitParam(2);
    
44. 
    
45. yModel=workMatrix(:,1)*fitParam(1)+fitParam(2);
    
46. 
    
47. xModel = workMatrix(:,1);
    
48. 
    
49. 
    
50. 
    
51. workMatrix(:,3)=yModel;
    
52. 
    
53. 
    
54. 
    
55. %calculate SSE,SSR,SST
    
56. 
    
57. workMatrix(:,4)=workMatrix(:,2)-workMatrix(:,3); %SSE
    
58. 
    
59. workMatrix(:,5)=workMatrix(:,2)-mean(workMatrix(:,2)); %SSR
    
60. 
    
61. workMatrix(:,6)=workMatrix(:,3)-mean(workMatrix(:,2));%SST
    
62. 
    
63. 
    
64. 
    
65. SSE=sum(workMatrix(:,4).^2);
    
66. 
    
67. SSR=sum(workMatrix(:,5).^2);
    
68. 
    
69. SST=sum(workMatrix(:,6).^2);
    
70. 
    
71. 
    
72. 
    
73. Rsq=1-SSE/SST;
    

复制代码

但是运行后，就卡住了，然后过一段时间以后就开始报错：

1. >> Test
   
2. Unable to perform assignment because the indices on the left side are not compatible with the size of the right side.
   
3. 
   
4. Error in RSQ (line 17)
   
5. workMatrix(:,1)=x;
   
6. 
   
7. Error in granCalc (line 12)
   
8. [slope, intercept,xModel, yModel,rsquare] = RSQ(volumes,ygran);
   
9. 
   
10. Error in Test>StartAutoTitration_Callback (line 937)
    
11. [TA, gamma, rsquare]=granCalc(sampleSize, Cacid, titrationLog(1:titrationStep-1,1)./1000000, titrationLog(1:titrationStep-1,3));
    
12. 
    
13. Error in gui_mainfcn (line 95)
    
14.         feval(varargin{:});
    
15. 
    
16. Error in Test (line 42)
    
17.     gui_mainfcn(gui_State, varargin{:});
    
18. 
    
19. Error in matlab.graphics.internal.figfile.FigFile/read>@(hObject,eventdata)Test('StartAutoTitration_Callback',hObject,eventdata,guidata(hObject))
    
20. 
    
21. Error using Test>emfRead_Callback (line 359)
    
22. Error while evaluating UIControl Callback.

复制代码

请教各位大神，这是什么原因出的错啊？是不是就与RSQ文件中的错误有关啊？应该怎么办啊？