2楼: Originally posted by Mr__Right at 2015-08-15 10:41:30
clear all ;
close all;
clc;
x0=40;
r=-3.79263*10^(-5)*x0^4+0.00856*x0^3-0.7236*x0^2+27.35571*x0-387.14149;
lb=x0-r; ub=x0+r;
x=linspace(lb,ub,300);
y1=-sqrt(r^2-x.^2+2*x0*x-x0^2);
y2=sqr ...

clear all ;
close all;
clc;
x0=40;
r=-3.79263*10^(-5)*x0^4+0.00856*x0^3-0.7236*x0^2+27.35571*x0-387.14149;
lb=x0-r; ub=x0+r;
x=linspace(lb,ub,300);
y1=-sqrt(r^2-x.^2+2*x0*x-x0^2);
y2=sqrt(r^2-x.^2+2*x0*x-x0^2);
scrsz = get(0,'ScreenSize');
hb=figure('Position',[150 50 scrsz(3)-200 scrsz(4)-200],'color','white');
h=plot(x,y1,'r',x,y2,'r'); axis([30,80,-12,12]) ;
for x0=40:.5:70, r=-3.79263*10^(-5)*x0^4+0.00856*x0^3-0.7236*x0^2+27.35571*x0-387.14149;
lb=x0-r;
ub=x0+r;
x=linspace(lb,ub,300);
y1=-sqrt(r^2-x.^2+2*x0*x-x0^2);
y2=sqrt(r^2-x.^2+2*x0*x-x0^2);
set(h(1),'xdata',x,'ydata',y1);
set(h(2),'xdata',x,'ydata',y2);     
pause(0.1) ;
end

3楼: Originally posted by anychi at 2015-08-15 14:31:21
谢谢您的回答。我还有一个小问题，就是想在上面的动画中添加一条y=0的水平线，并且一直存在，请问代码应该如何修改一下。...

h=plot(......................);

h=plot(x,y1,'r',x,y2,'r',[0,100],[0,0],'g:');

4楼: Originally posted by Mr__Right at 2015-08-15 16:11:53
把
h=plot(......................);
改成
h=plot(x,y1,'r',x,y2,'r',,,'g:');...