% --- Executes on button press in calculate. function calculate_Callback(hObject, eventdata, handles) % hObject handle to calculate (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) set(handles.state,'string','程序正在运行......'); ZXJO=str2num(get(handles.ZXJO,'string')); ZXJT=str2num(get(handles.ZXJT,'string')); DZPBJ=str2num(get(handles.DZPBJ,'string')); JYBJ=str2num(get(handles.JYBJ,'string')); PZJO=str2num(get(handles.PZJO,'string')); PZJT=str2num(get(handles.PZJT,'string')); DCDO=str2num(get(handles.DCDO,'string')); DJD=str2num(get(handles.DJD,'string')); DCDT=str2num(get(handles.DCDT,'string')); SXJDO=str2num(get(handles.SXJDO,'string')); SXJDT=str2num(get(handles.SXJDT,'string')); JPJDO=str2num(get(handles.JPJDO,'string')); JPJDT=str2num(get(handles.JPJDT,'string')); JPCDO=str2num(get(handles.JPCDO,'string')); JPCDT=str2num(get(handles.JPCDT,'string')); CPJDO=str2num(get(handles.CPJDO,'string')); CPJDT=str2num(get(handles.CPJDT,'string')); CPCDO=str2num(get(handles.CPCDO,'string')); CPCDT=str2num(get(handles.CPCDT,'string')); SZPBJ=str2num(get(handles.SZPBJ,'string')); syms t; t=0:0.0001:1.0002; tou=0.0001; N=size(t,2); if get(handles.zhengxian,'value')==1 f1=t-sin(180.*t*pi/180)/pi; f2=1-t-(sin(180.*t*pi/180))/pi; else f1=1-cos((180.*t/2)*pi/180); f2=1-sin((180.*t/2)*pi/180); end q1=-ZXJO.*sin((-SXJDO/2+SXJDO.*t-180+DJD)*pi/180)-DZPBJ.*sin((-180-3*SXJDO/4+3*SXJDO.*t/2-180+DJD)*pi/180)-JYBJ.*sin((180-DJD)*pi/180); q2=ZXJO.*cos((-SXJDO/2+SXJDO.*t-180+DJD)*pi/180)+DZPBJ.*cos((-180-3*SXJDO/4+3*SXJDO.*t/2-180+DJD)*pi/180)-JYBJ.*cos((180-DJD)*pi/180); a=atan(q1./q2)*180/pi; q3=DCDO.*sin((180-a)*pi/180); q4=JYBJ+DCDO.*cos((180-a)*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+DCDO.*cos((180-a-b)*pi/180); theta=90-PZJO+SXJDO/2-SXJDO.*t+b; axes(handles.axes1); ax1=handles.axes1; cla(ax1); polar(theta*pi/180,r,'-b'); hold on; hold on; polar([0:0.01:2*pi],ones(size(0:0.01:2*pi))*400,'b-'); Store_1(:,1)=[r.*cos(theta*pi/180)]'; Store_1(:,2)=[r.*sin(theta*pi/180)]'; save Store_1.xls -ascii Store_1 q1=-ZXJO.*sin((SXJDO/2+CPJDO.*t-180+104.93)*pi/180)-(DZPBJ+CPCDO.*f1).*sin((-180+3*SXJDO/4+3*CPJDO.*t/2-180+104.93)*pi/180)+JYBJ.*sin((-180+104.93)*pi/180); q2=ZXJO.*cos((SXJDO/2+CPJDO.*t-180+104.93)*pi/180)+(DZPBJ+CPCDO.*f1).*cos((-180+3*SXJDO/4+3*CPJDO.*t/2-180+104.93)*pi/180)-JYBJ.*cos((-180+104.93)*pi/180); a=atan(q1./q2)*180/pi; q3=DCDO.*sin((180-a)*pi/180); q4=JYBJ+DCDO.*cos((180-a)*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+DCDO.*cos((180-a-b)*pi/180); theta=90-PZJO-SXJDO/2-CPJDO.*t+b; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store_2(:,1)=[r.*cos(theta*pi/180)]'; Store_2(:,2)=[r.*sin(theta*pi/180)]'; save Store_2.txt -ascii Store_2 n=size(t,2); for i=1:n-1 dr(i)=(r(i+1)-r(i))/tou; end for i=1:n-2 d2r(i)=(dr(i+1)-dr(i))/tou; end m=size(t,2); for j=1:m-1 tr(j)=(theta(j+1)-theta(j))/tou; end for j=1:m-2 t2r(j)=(tr(j+1)-tr(j))/tou; end w1=r(end-2); w3=dr(end-1); w5=d2r(end); y1=theta(end-2); y3=tr(end-1); y5=t2r(end); q1=ZXJT.*sin((-SXJDT/2-JPJDT+JPJDT.*t)*pi/180)+(SZPBJ+JPCDT.*f2).*sin((180-3*SXJDT/4-3*JPJDT/2+3*JPJDT.*t/2)*pi/180); q2=ZXJT.*cos((-SXJDT/2-JPJDT+JPJDT.*t)*pi/180)+(SZPBJ+JPCDT.*f2).*cos((180-3*SXJDT/4-3*JPJDT/2+3*JPJDT.*t/2)*pi/180)-JYBJ; a=atan(q1./q2)*180/pi; q3=DCDO.*sin((DJD+a)*pi/180); q4=JYBJ+DCDO.*cos((DJD+a)*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+DCDO.*cos((DJD+a-b)*pi/180); theta=-90-PZJT+SXJDT/2+JPJDT-JPJDT.*t+b; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(2*N+1:3*N,1)=[r.*cos(theta*pi/180)]'; Store(2*N+1:3*N,2)=[r.*sin(theta*pi/180)]'; n=size(t,2); for i=1:n-1 dr(i)=(r(i+1)-r(i))/tou; end for i=1:n-2 d2r(i)=(dr(i+1)-dr(i))/tou; end m=size(t,2); for j=1:m-1 tr(j)=(theta(j+1)-theta(j))/tou; end for j=1:m-2 t2r(j)=(tr(j+1)-tr(j))/tou; end w2=r(1); w4=dr(1); w6=d2r(1); y2=theta(1); y4=tr(1); y6=t2r(1); syms x C1=[x.^5 x.^4 x.^3 x.^2 x 1]; C2=[5*x.^4 4*x.^3 3*x.^2 2*x 1 0]; C3=[5*4*x.^3 4*3*x.^2 3*2*x 2 0 0]; A=[subs(C1,0);subs(C1,1);subs(C2,0);subs(C2,1);subs(C3,0);subs(C3,1)]; w=[w1 w2 w3 w4 w5 w6]'; e=inv(A)*w; r=e(1).*t.^5+e(2).*t.^4+e(3).*t.^3+e(4).*t.^2+e(5).*t+e(6); e_r=e; y=[y1 y2 y3 y4 y5 y6]'; e=inv(A)*y; theta=e(1).*t.^5+e(2).*t.^4+e(3).*t.^3+e(4).*t.^2+e(5).*t+e(6); e_theta=e; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(3*N+1:4*N,1)=[r.*cos(theta*pi/180)]'; Store(3*N+1:4*N,2)=[r.*sin(theta*pi/180)]'; q1=-ZXJT.*sin((-SXJDT/2+SXJDT.*t-180+DJD)*pi/180)-SZPBJ.*sin((180-3*SXJDT/4+3*SXJDT.*t/2-180+DJD)*pi/180)+JYBJ.*sin((-180+DJD)*pi/180); q2=ZXJT.*cos((-SXJDT/2+SXJDT.*t-180+DJD)*pi/180)+SZPBJ.*cos((180-3*SXJDT/4+3*SXJDT.*t/2-180+DJD)*pi/180)-JYBJ.*cos((-180+DJD)*pi/180); a=atan(q1./q2)*180/pi; q3=DCDO.*sin((180-a)*pi/180); q4=JYBJ+DCDO.*cos((180-a)*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+DCDO.*cos((180-a-b)*pi/180); theta=-90-PZJT+SXJDT/2-SXJDT.*t+b; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(4*N+1:5*N,1)=[r.*cos(theta*pi/180)]'; Store(4*N+1:5*N,2)=[r.*sin(theta*pi/180)]'; q1=-ZXJT.*sin((SXJDT/2+CPJDT.*t-180+DJD)*pi/180)-SZPBJ.*sin((180+3*SXJDT/4+3*CPJDT.*t/2-180+DJD)*pi/180)+JYBJ.*sin((-180+DJD)*pi/180); q2=ZXJT.*cos((SXJDT/2+CPJDT.*t-180+DJD)*pi/180)+SZPBJ.*cos((180+3*SXJDT/4+3*CPJDT.*t/2-180+DJD)*pi/180)-JYBJ.*cos((-180+DJD)*pi/180); a=atan(q1./q2)*180/pi; q3=DCDO.*sin((180-a)*pi/180); q4=JYBJ+DCDO.*cos((180-a)*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+DCDO.*cos((180-a-b)*pi/180); theta=-90-PZJT-SXJDT/2-CPJDT.*t+b; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(5*N+1:6*N,1)=[r.*cos(theta*pi/180)]'; Store(5*N+1:6*N,2)=[r.*sin(theta*pi/180)]'; n=size(t,2); for i=1:n-1 dr(i)=(r(i+1)-r(i))/tou; end for i=1:n-2 d2r(i)=(dr(i+1)-dr(i))/tou; end m=size(t,2); for j=1:m-1 tr(j)=(theta(j+1)-theta(j))/tou; end for j=1:m-2 t2r(j)=(tr(j+1)-tr(j))/tou; end w1=r(end-2); w3=dr(end-1); w5=d2r(end); y1=theta(end-2); y3=tr(end-1); y5=t2r(end); q1=-ZXJO.*sin((-SXJDO/2-JPJDO+JPJDO.*t-180+DJD)*pi/180)-DZPBJ.*sin((-180-3*SXJDO/4-3*JPJDO/2+3*JPJDO.*t/2-180+DJD)*pi/180)+JYBJ.*sin((-180+DJD)*pi/180); q2=ZXJO.*cos((-SXJDO/2-JPJDO+JPJDO.*t-180+DJD)*pi/180)+DZPBJ.*cos((-180-3*SXJDO/4-3*JPJDO/2+3*JPJDO.*t/2-180+DJD)*pi/180)-JYBJ.*cos((-180+DJD)*pi/180); a=atan(q1./q2)*180/pi; q3=DCDO.*sin((180-a)*pi/180); q4=JYBJ+DCDO.*cos((180-a)*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+DCDO.*cos((180-a-b)*pi/180); theta=90-PZJO+SXJDO/2+JPJDO-JPJDO.*t+b-360; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(6*N+1:7*N,1)=[r.*cos(theta*pi/180)]'; Store(6*N+1:7*N,2)=[r.*sin(theta*pi/180)]'; n=size(t,2); for i=1:n-1 dr(i)=(r(i+1)-r(i))/tou; end for i=1:n-2 d2r(i)=(dr(i+1)-dr(i))/tou; end m=size(t,2); for j=1:m-1 tr(j)=(theta(j+1)-theta(j))/tou; end for j=1:m-2 t2r(j)=(tr(j+1)-tr(j))/tou; end w2=r(1); w4=dr(1); w6=d2r(1); y2=theta(1); y4=tr(1); y6=t2r(1); syms x C1=[x.^5 x.^4 x.^3 x.^2 x 1]; C2=[5*x.^4 4*x.^3 3*x.^2 2*x 1 0]; C3=[5*4*x.^3 4*3*x.^2 3*2*x 2 0 0]; A=[subs(C1,0);subs(C1,1);subs(C2,0);subs(C2,1);subs(C3,0);subs(C3,1)]; w=[w1 w2 w3 w4 w5 w6]'; e=inv(A)*w; r=e(1).*t.^5+e(2).*t.^4+e(3).*t.^3+e(4).*t.^2+e(5).*t+e(6); e_r(:,2)=e; y=[y1 y2 y3 y4 y5 y6]'; e=inv(A)*y; theta=e(1).*t.^5+e(2).*t.^4+e(3).*t.^3+e(4).*t.^2+e(5).*t+e(6); e_theta(:,2)=e; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(7*N+1:8*N,1)=[r.*cos(theta*pi/180)]'; Store(7*N+1:8*N,2)=[r.*sin(theta*pi/180)]'; q1=-ZXJO.*sin((-SXJDO/2+SXJDO.*t-180+DJD)*pi/180)-DZPBJ.*sin((-180-3*SXJDO/4+3*SXJDO.*t/2-180+DJD)*pi/180)-JYBJ.*sin((180-DJD)*pi/180); q2=ZXJO.*cos((-SXJDO/2+SXJDO.*t-180+DJD)*pi/180)+DZPBJ.*cos((-180-3*SXJDO/4+3*SXJDO.*t/2-180+DJD)*pi/180)-JYBJ.*cos((180-DJD)*pi/180); a=atan(q1./q2)*180/pi; P=ZXJO.*cos((-SXJDO/2+SXJDO.*t-180+DJD+a)*pi/180)+DZPBJ.*cos((-180-3*SXJDO/4+3*SXJDO.*t/2-180+DJD+a)*pi/180)-JYBJ.*cos((-180+DJD+a)*pi/180); q3=(P-DCDT).*sin((180-a-104.93)*pi/180); q4=JYBJ+(P-DCDT).*cos((180-a-DJD)*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+(P-DCDT).*cos((180-a-DJD-b)*pi/180); theta=90-PZJO+SXJDO/2-SXJDO.*t+b-360; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(8*N+1:9*N,1)=[r.*cos(theta*pi/180)]'; Store(8*N+1:9*N,2)=[r.*sin(theta*pi/180)]'; q1=-ZXJO.*sin((SXJDO/2+CPJDO.*t-180+104.93)*pi/180)-(DZPBJ+CPCDO.*f1).*sin((-180+3*SXJDO/4+3*CPJDO.*t/2-180+104.93)*pi/180)+JYBJ.*sin((-180+104.93)*pi/180); q2=ZXJO.*cos((SXJDO/2+CPJDO.*t-180+104.93)*pi/180)+(DZPBJ+CPCDO.*f1).*cos((-180+3*SXJDO/4+3*CPJDO.*t/2-180+104.93)*pi/180)-JYBJ.*cos((-180+104.93)*pi/180); a=atan(q1./q2)*180/pi; P=ZXJO.*cos((SXJDO/2+CPJDO.*t-180+104.93+a)*pi/180)+(DZPBJ+CPCDO.*f1).*cos((-180+3*SXJDO/4+3*CPJDO.*t/2-180+104.93+a)*pi/180)-JYBJ.*cos((-180+104.93+a)*pi/180); q3=(P-DCDT).*sin((180-DJD-a)*pi/180); q4=JYBJ+(P-DCDT).*cos((180-DJD-a)*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+(P-DCDT).*cos((180-DJD-a-b)*pi/180); theta=90-PZJO-SXJDO/2-CPJDO.*t+b; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(9*N+1:10*N,1)=[r.*cos(theta*pi/180)]'; Store(9*N+1:10*N,2)=[r.*sin(theta*pi/180)]'; n=size(t,2); for i=1:n-1 dr(i)=(r(i+1)-r(i))/tou; end for i=1:n-2 d2r(i)=(dr(i+1)-dr(i))/tou; end m=size(t,2); for j=1:m-1 tr(j)=(theta(j+1)-theta(j))/tou; end for j=1:m-2 t2r(j)=(tr(j+1)-tr(j))/tou; end w1=r(end-2); w3=dr(end-1); w5=d2r(end); y1=theta(end-2); y3=tr(end-1); y5=t2r(end); q1=ZXJT.*sin((-SXJDT/2-JPJDT+JPJDT.*t)*pi/180)+(SZPBJ+JPCDT.*f2).*sin((180-3*SXJDT/4-3*JPJDT/2+3*JPJDT.*t/2)*pi/180); q2=ZXJT.*cos((-SXJDT/2-JPJDT+JPJDT.*t)*pi/180)+(SZPBJ+JPCDT.*f2).*cos((180-3*SXJDT/4-3*JPJDT/2+3*JPJDT.*t/2)*pi/180)-JYBJ; a=atan(q1./q2)*180/pi; P=ZXJT.*cos((-SXJDT/2-JPJDT+JPJDT.*t-a)*pi/180)+(SZPBJ+JPCDT.*f2).*cos((180-3*SXJDT/4-3*JPJDT/2+3*JPJDT.*t/2-a)*pi/180)-JYBJ.*cos(a*pi/180); q3=(P-DCDT).*sin(a*pi/180); q4=JYBJ+(P-DCDT).*cos(a*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+(P-DCDT).*cos((a-b)*pi/180); theta=-90-PZJT+SXJDT/2+JPJDT-JPJDT.*t+b; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(10*N+1:11*N,1)=[r.*cos(theta*pi/180)]'; Store(10*N+1:11*N,2)=[r.*sin(theta*pi/180)]'; n=size(t,2); for i=1:n-1 dr(i)=(r(i+1)-r(i))/tou; end for i=1:n-2 d2r(i)=(dr(i+1)-dr(i))/tou; end m=size(t,2); for j=1:m-1 tr(j)=(theta(j+1)-theta(j))/tou; end for j=1:m-2 t2r(j)=(tr(j+1)-tr(j))/tou; end w2=r(1); w4=dr(1); w6=d2r(1); y2=theta(1); y4=tr(1); y6=t2r(1); syms x C1=[x.^5 x.^4 x.^3 x.^2 x 1]; C2=[5*x.^4 4*x.^3 3*x.^2 2*x 1 0]; C3=[5*4*x.^3 4*3*x.^2 3*2*x 2 0 0]; A=[subs(C1,0);subs(C1,1);subs(C2,0);subs(C2,1);subs(C3,0);subs(C3,1)]; w=[w1 w2 w3 w4 w5 w6]'; e=inv(A)*w; r=e(1).*t.^5+e(2).*t.^4+e(3).*t.^3+e(4).*t.^2+e(5).*t+e(6); e_r(:,3)=e; y=[y1 y2 y3 y4 y5 y6]'; e=inv(A)*y; theta=e(1).*t.^5+e(2).*t.^4+e(3).*t.^3+e(4).*t.^2+e(5).*t+e(6); e_theta(:,3)=e; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(11*N+1:12*N,1)=[r.*cos(theta*pi/180)]'; Store(11*N+1:12*N,2)=[r.*sin(theta*pi/180)]'; q1=-ZXJT.*sin((-SXJDT/2+SXJDT.*t-180+DJD)*pi/180)-SZPBJ.*sin((180-3*SXJDT/4+3*SXJDT.*t/2-180+DJD)*pi/180)+JYBJ.*sin((-180+DJD)*pi/180); q2=ZXJT.*cos((-SXJDT/2+SXJDT.*t-180+DJD)*pi/180)+SZPBJ.*cos((180-3*SXJDT/4+3*SXJDT.*t/2-180+DJD)*pi/180)-JYBJ.*cos((-180+DJD)*pi/180); a=atan(q1./q2)*180/pi; P=ZXJT.*cos((-SXJDT/2+SXJDT.*t-180+DJD+a)*pi/180)+SZPBJ.*cos((180-3*SXJDT/4+3*SXJDT.*t/2-180+DJD+a)*pi/180)-JYBJ.*cos((-180+DJD+a)*pi/180); q3=(P-DCDT).*sin((180-a-DJD)*pi/180); q4=JYBJ+(P-DCDT).*cos((180-a-DJD)*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+(P-DCDT).*cos((180-a-DJD-b)*pi/180); theta=-90-PZJT+SXJDT/2-SXJDT.*t+b; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(12*N+1:13*N,1)=[r.*cos(theta*pi/180)]'; Store(12*N+1:13*N,2)=[r.*sin(theta*pi/180)]'; q1=-ZXJT.*sin((SXJDT/2+CPJDT.*t-180+DJD)*pi/180)-SZPBJ.*sin((180+3*SXJDT/4+3*CPJDT.*t/2-180+DJD)*pi/180)+JYBJ.*sin((-180+DJD)*pi/180); q2=ZXJT*cos((SXJDT/2+CPJDT*t-180+DJD)*pi/180)+SZPBJ*cos((180+3*SXJDT/4+3*CPJDT*t/2-180+DJD)*pi/180)-JYBJ.*cos((-180+DJD)*pi/180); a=atan(q1./q2)*180/pi; P=ZXJT.*cos((SXJDT/2+CPJDT.*t-180+DJD+a)*pi/180)+SZPBJ.*cos((180+3*SXJDT/4+3*CPJDT.*t/2-180+DJD+a)*pi/180)-JYBJ.*cos((-180+DJD+a)*pi/180); q3=(P-DCDT-CPCDT.*f1).*sin((180-a-DJD)*pi/180); q4=JYBJ+(P-DCDT-CPCDT.*f1).*cos((180-a-DJD)*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+(P-DCDT-CPCDT.*f1).*cos((180-a-DJD-b)*pi/180); theta=-90-PZJT-SXJDT/2-CPJDT.*t+b; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(13*N+1:14*N,1)=[r.*cos(theta*pi/180)]'; Store(13*N+1:14*N,2)=[r.*sin(theta*pi/180)]'; n=size(t,2); for i=1:n-1 dr(i)=(r(i+1)-r(i))/tou; end for i=1:n-2 d2r(i)=(dr(i+1)-dr(i))/tou; end m=size(t,2); for j=1:m-1 tr(j)=(theta(j+1)-theta(j))/tou; end for j=1:m-2 t2r(j)=(tr(j+1)-tr(j))/tou; end w1=r(end-2); w3=dr(end-1); w5=d2r(end); y1=theta(end-2); y3=tr(end-1); y5=t2r(end); q1=-ZXJO.*sin((-SXJDO/2-JPJDO+JPJDO.*t-180+DJD)*pi/180)-DZPBJ.*sin((-180-3*SXJDO/4-3*JPJDO/2+3*JPJDO.*t/2-180+DJD)*pi/180)+JYBJ.*sin((-180+DJD)*pi/180); q2=ZXJO.*cos((-SXJDO/2-JPJDO+JPJDO.*t-180+DJD)*pi/180)+DZPBJ.*cos((-180-3*SXJDO/4-3*JPJDO/2+3*JPJDO.*t/2-180+DJD)*pi/180)-JYBJ.*cos((-180+DJD)*pi/180); a=atan(q1./q2)*180/pi; P=ZXJO.*cos((-SXJDO/2-JPJDO+JPJDO.*t-180+DJD+a)*pi/180)+DZPBJ.*cos((-180-3*SXJDO/4-3*JPJDO/2+3*JPJDO.*t/2-180+DJD+a)*pi/180)-JYBJ.*cos((-180+DJD+a)*pi/180); q3=(P-DCDT-JPCDO.*f2).*sin((180-a-DJD)*pi/180); q4=JYBJ+(P-DCDT-JPCDO.*f2).*cos((180-a-DJD)*pi/180); b=atan(q3./q4)*180/pi; r=JYBJ.*cos(b*pi/180)+(P-DCDT-JPCDO.*f2).*cos((180-a-DJD-b)*pi/180); theta=90-PZJO+SXJDO/2+JPJDO-JPJDO.*t+b-360; axes(handles.axes1); polar(theta*pi/180,r,'-b'); Store(14*N+1:15*N,1)=[r.*cos(theta*pi/180)]'; Store(14*N+1:15*N,2)=[r.*sin(theta*pi/180)]'; n=size(t,2); for i=1:n-1 dr(i)=(r(i+1)-r(i))/tou; end for i=1:n-2 d2r(i)=(dr(i+1)-dr(i))/tou; end m=size(t,2); for j=1:m-1 tr(j)=(theta(j+1)-theta(j))/tou; end for j=1:m-2 t2r(j)=(tr(j+1)-tr(j))/tou; end w2=r(1); w4=dr(1); w6=d2r(1); y2=theta(1); y4=tr(1); y6=t2r(1); syms x C1=[x.^5 x.^4 x.^3 x.^2 x 1]; C2=[5*x.^4 4*x.^3 3*x.^2 2*x 1 0]; C3=[5*4*x.^3 4*3*x.^2 3*2*x 2 0 0]; A=[subs(C1,0);subs(C1,1);subs(C2,0);subs(C2,1);subs(C3,0);subs(C3,1)]; w=[w1 w2 w3 w4 w5 w6]'; e=inv(A)*w; r=e(1).*t.^5+e(2).*t.^4+e(3).*t.^3+e(4).*t.^2+e(5).*t+e(6); e_r(:,4)=e; y=[y1 y2 y3 y4 y5 y6]'; e=inv(A)*y; e_theta(:,4)=e; theta=e(1).*t.^5+e(2).*t.^4+e(3).*t.^3+e(4).*t.^2+e(5).*t+e(6); polar(theta*pi/180,r,'-b'); axis([-800 800 -800 800]); hold off; Store(15*N+1:16*N,1)=[r.*cos(theta*pi/180)]'; Store(15*N+1:16*N,2)=[r.*sin(theta*pi/180)]'; save Store.txt -ascii Store set(handles.state,'string','程序运行完成'); set(handles.banjing,'data',e_r'); set(handles.angle,'data',e_theta');
时间: 2025-07-19 22:26:27 浏览: 1
<think>我们正在讨论一个MATLAB GUI回调函数`calculate_Callback`,该函数用于执行复杂的数学计算和极坐标绘图。根据用户提供的错误信息,主要问题是在初始化或数据赋值时出现了“矩形空矩阵的不正确赋值”错误。
错误信息分析:
引用[1]中的错误:
??? Improper assignment with rectangular empty matrix.
Error in ==> chengjichaxun>chengjichaxun_OpeningFcn at 69
chengji(counter-1,:)= [yuwen shuxue yingyu];
这个错误通常发生在尝试将一个矩阵(或向量)赋值给一个数组的某一行(或列),但目标数组的大小与赋值内容不匹配时。具体来说,可能是:
- 目标数组`chengji`在赋值之前没有被正确初始化,或者初始化的维度与赋值内容不一致。
- 赋值时右侧的向量长度与目标数组的列数不一致。
另外,用户还提到了一个文件打开回调函数的例子(引用[2]),但当前问题集中在计算回调函数上。
由于用户希望了解`calculate_Callback`函数,特别是涉及复杂数学计算和极坐标绘图,我们将重点放在如何正确实现这样的回调函数,并避免常见的错误。
步骤:
1. 获取输入数据:从GUI控件(如编辑框)中获取用户输入。
2. 数据验证:确保输入数据是有效的(例如,数字而非字符串,在合理范围内)。
3. 执行计算:进行复杂的数学计算。
4. 绘图:在极坐标下绘图。
5. 错误处理:避免出现类似引用[1]中的错误。
下面是一个示例的`calculate_Callback`函数,用于执行复杂计算(例如,计算一个基于极坐标的数学函数)并绘制极坐标图。
注意:在MATLAB中绘制极坐标图可以使用`polarplot`函数(R2016a及以上版本)或`polar`函数(旧版本)。
示例代码:
```matlab
function calculate_Callback(hObject, eventdata, handles)
% 1. 从编辑框中获取输入参数
% 假设有两个参数:幅度(a)和频率(b)
try
a = str2double(get(handles.edit_amplitude, 'String'));
b = str2double(get(handles.edit_frequency, 'String'));
catch
errordlg('请输入有效的数字!', '输入错误');
return;
end
% 2. 验证输入是否为数字
if isnan(a) || isnan(b)
errordlg('输入必须为数字!', '输入错误');
return;
end
% 3. 定义极坐标的角度变量theta(0到2π)
theta = linspace(0, 2*pi, 1000);
% 4. 计算极径rho(这里是一个示例函数,可以是任意复杂的数学表达式)
% 示例:一个带有幅度和频率的玫瑰曲线
rho = a * sin(b * theta);
% 5. 在极坐标中绘图
% 检查handles中是否有axes对象,如果没有则创建,或者使用已有的axes
if isfield(handles, 'axes_plot')
axes(handles.axes_plot); % 激活该坐标轴
else
handles.axes_plot = axes('Parent', handles.figure1); % 假设figure1是主窗口
% 或者使用gcf,但最好使用handles中的figure句柄
end
% 绘制极坐标图(使用polarplot,适用于新版本)
% 注意:如果使用旧版本MATLAB,则使用polar函数,但polar函数在R2016a之后不推荐使用
% 这里使用polarplot(推荐)
polarplot(handles.axes_plot, theta, rho);
title(handles.axes_plot, '极坐标图');
% 6. 保存数据到handles结构体,以便其他回调使用
handles.theta = theta;
handles.rho = rho;
guidata(hObject, handles); % 更新handles
% 7. 如果没有错误,可以在这里显示成功信息
set(handles.text_status, 'String', '计算完成!');
end
```
避免引用[1]中的错误:
错误`Improper assignment with rectangular empty matrix`通常发生在数组赋值时维度不匹配。在我们的计算回调函数中,如果涉及到数组的构建,需要注意:
- 在初始化数组时,应确保其维度与后续赋值一致。
- 如果使用动态扩展数组(例如,在循环中不断添加数据),应尽量预分配数组以提高效率并避免维度问题。
例如,如果我们需要在回调函数中构建一个数组,可以这样做:
```matlab
% 假设我们要存储多组计算结果,但通常我们不会在回调中这样存储,除非有特殊需求
% 首先检查handles中是否已有该数组
if ~isfield(handles, 'resultsArray')
% 初始化数组
handles.resultsArray = [];
end
% 将新结果添加到数组(假设result是一个行向量)
newResult = [a, b, max(rho)]; % 示例:存储输入参数和rho的最大值
handles.resultsArray = [handles.resultsArray; newResult]; % 垂直连接
% 注意:如果每次回调只产生一组数据(一行),那么这样添加是可行的。但要确保每次添加的行向量长度一致(列数相同)。
```
但是,如果newResult的列数与handles.resultsArray的列数不一致,就会报错。因此,在初始化数组时,最好确定列数,并在每次添加时检查。
另一种方法是使用元胞数组,因为元胞数组可以容纳不同大小的数据。
关于极坐标绘图的注意事项:
- 在旧版MATLAB中,使用`polar`函数,但该函数创建的图形是单独的,不易嵌入到GUI的坐标轴中。新版本中,`polarplot`可以直接在指定坐标轴中绘制。
- 如果使用旧版,可以尝试:`polar(theta, rho);`,但这会创建一个新的图形窗口,所以需要先指定当前坐标轴(axes)并设置其属性为极坐标。
在旧版MATLAB中嵌入极坐标图到GUI坐标轴的方法:
```matlab
% 对于旧版MATLAB(R2014b及之前),可以使用以下技巧:
% 首先,将坐标轴转换为极坐标
axes(handles.axes_plot); % 激活坐标轴
% 然后使用polar,但注意polar会创建新的图形,所以我们可以捕获其输出:
h = polar(theta, rho);
% 将h的父对象设置为handles.axes_plot? 实际上,polar会创建新的坐标轴,所以我们需要将图形对象复制到我们的坐标轴中,这比较复杂。
% 替代方案:使用极坐标到笛卡尔坐标的转换,然后在普通坐标轴中绘制
x = rho .* cos(theta);
y = rho .* sin(theta);
plot(handles.axes_plot, x, y);
axis(handles.axes_plot, 'equal'); % 使坐标轴等比例,以显示圆形
title(handles.axes_plot, '极坐标图(笛卡尔坐标表示)');
```
但建议使用新版本的`polarplot`,因为它更直接。
调试技巧:
- 在可能出现错误的地方(如数组赋值)设置断点。
- 检查赋值语句左右两侧的变量大小(使用`size`函数)。
- 使用`try-catch`捕获错误并显示详细信息。
例如:
```matlab
try
chengji(counter-1,:)= [yuwen shuxue yingyu];
catch ME
disp(ME.message);
disp(['目标数组大小: ', num2str(size(chengji))]);
disp(['赋值数组大小: ', num2str(size([yuwen shuxue yingyu]))]);
end
```
总结:
在`calculate_Callback`中,重点在于安全地获取输入、进行数学计算和绘图。同时,注意数组操作的维度匹配,避免出现引用[1]中的错误。
相关问题:
1. 如何在MATLAB GUI中初始化一个数组,以便在多个回调函数中使用?
2. 在MATLAB GUI中绘制极坐标图有哪些方法?如何兼容不同版本的MATLAB?
3. 如何调试MATLAB GUI回调函数中的数组维度错误?
4. 在回调函数中,如何安全地存储和传递数据,避免数据丢失?</think>### MATLAB GUI `calculate_Callback` 函数详解与调试指南
#### 一、核心功能解析
`calculate_Callback` 是 MATLAB GUI 的核心计算函数,通常关联计算按钮,负责:
1. **数据获取**:从输入控件读取参数
2. **数学计算**:执行复杂数学运算(如极坐标变换)
3. **结果可视化**:在坐标轴绘制极坐标图等
4. **错误处理**:验证输入数据有效性
#### 二、典型错误分析(引用[1]问题)
错误 `Improper assignment with rectangular empty matrix` 通常由以下原因引起:
```matlab
chengji(counter-1,:) = [yuwen shuxue yingyu]; % 错误行
```
- **根本原因**:
- 数组 `chengji` 未预分配或维度不匹配
- 当 `counter-1 ≤ 0` 时索引无效
- 右侧向量长度 ≠ `chengji` 的列数
- **解决方案**:
```matlab
% 1. 预分配数组(确定最大行数和列数)
maxRows = 10; numCols = 3;
chengji = zeros(maxRows, numCols); % 预分配
% 2. 验证索引有效性
if counter > 1 && counter <= maxRows
chengji(counter-1, :) = [yuwen, shuxue, yingyu];
else
errordlg('索引超出范围', '错误');
end
```
#### 三、完整实现框架(含极坐标绘图)
```matlab
function calculate_Callback(hObject, ~, handles)
% 1. 获取输入参数
try
r = str2double(get(handles.edit_radius, 'String'));
theta = deg2rad(str2double(get(handles.edit_angle, 'String')));
catch
errordlg('输入必须为数字', '格式错误');
return;
end
% 2. 参数验证
if any(isnan([r, theta])) || r <= 0
errordlg('半径必须为正数', '输入错误');
return;
end
% 3. 复杂数学计算(极坐标变换)
phi = linspace(0, 2*pi, 500); % 角度范围$[0, 2\pi]$
rho = r * sin(3*phi + theta); % 玫瑰曲线方程
% 4. 极坐标绘图
polarplot(handles.axes_main, phi, rho, 'LineWidth', 2);
title(handles.axes_main, '极坐标函数可视化');
% 5. 保存计算结果
handles.lastCalculation.rho = rho;
handles.lastCalculation.phi = phi;
guidata(hObject, handles); % 关键!保存数据
end
```
#### 四、关键调试技巧
1. **断点设置**:
- 在行号左侧点击设置红点断点
- 使用 `dbstop if error` 自动在错误处暂停
2. **变量检查**:
```matlab
% 在回调函数内添加检查语句
disp(['r = ', num2str(r)]); % 显示半径值
whos chengji % 检查数组维度
```
3. **错误预防**:
- 所有数值输入使用 `str2double()` 而非 `str2num`
- 数组操作前检查维度:`assert(size(A,2)==size(B,2))`
- 使用 `try/catch` 捕获文件操作错误(如引用[2]的文件读取)
4. **极坐标绘图注意事项**:
- 新版本MATLAB使用 `polarplot` 替代旧版 `polar`
- 角度需转换为弧度:`deg2rad()` 或 `*pi/180`
- 设置坐标轴句柄:`polarplot(handles.axes1, ...)`
#### 五、典型问题解决方案
| 问题现象 | 解决方案 |
|--------------------------|------------------------------------------------------------------------|
| 空矩阵赋值错误 | 预分配数组 + 索引范围检查 |
| 极坐标图不显示 | 确认使用 `handles.axes_name` 指定坐标轴 |
| 数据未保存 | 确保调用 `guidata(hObject, handles)` |
| 文件读取失败(引用[2]) | 检查文件路径拼接:`fullfile(filepath, filename)` |
> **最佳实践**:复杂计算分离为独立函数
> ```matlab
> % 独立计算函数
> function [rho, phi] = polar_calc(r, theta, n)
> phi = linspace(0, 2*pi, 500);
> rho = r * sin(n*phi + theta);
> end
>
> % 回调函数中调用
> [rho, phi] = polar_calc(r, theta, 3);
> ```
#### 六、扩展应用:动态参数更新
```matlab
% 添加滑块回调更新计算
function slider_Callback(hObject, ~, handles)
n = get(hObject, 'Value');
set(handles.text_param, 'String', ['n=', num2str(n)]);
calculate_Callback(handles.calc_button, [], handles); % 触发重计算
end
```
阅读全文
相关推荐












