Export grafů do souboru + grafické formáty
> animate( sin(x*t),x=0..10,t=1..2,frames=50,numpoints=100,color=red);
> animate( [u*t,t,t=1..8*Pi], u=1..4,coords=polar,frames=60,numpoints=100);
> animate(10000/99*sqrt(99)*exp(tl*t)*sin(100*t*sqrt(99)),t=0..0.05,tl=-500..-50,frames=15,numpoints=100,title=`Impulzní charakteristika`,labelfont=[HELVETICA,8],axesfont=[HELVETICA,8]);
> restart;with(plots):
Warning, the name changecoords has been redefined
> Pk := K/(p^2-p*omega[0]/Q+omega[0]^2);
>
Pk_abs:=20*log10(evalc(abs(subs({K=1e6,omega[0]=1000,p=I*omega},Pk)))):
Pkf:=unapply(Pk_abs,omega,Q):
>
bb:=proc(x,y,n)
local i,k;
i:=(log10(y)-log10(x))/(n-1);
> seq(evalhf(10^(log10(x)+(k-1)*i)),k=1..n);
> end:
>
for k in 0.1,0.2,0.5,0.707,1,2,5,10 do
for i from 1 to 50 do
a:=bb(10,1e5,50)[i]:
body[i]:=[a,evalhf(Pkf(a,k))]:
od:
bodyk[k]:=[seq(body[l],l=1..10)]:
od:
>
for k in 0.1,0.2,0.5,0.707,1,2,5,10 do
body[k]:=[seq([bb(100,1e4,100)[i],evalf(Pkf(bb(100,1e4,100)[i],k))],i=1..100)]:od:
> body[];
> semilogplot(body[10],color=red);
>
>
>