TeXclipで作成した数式の再編集については
を参照してください.
function [ csys ] = obsvcanon(sys) %OBSVCANON Compute observability canonical form [n,d]=tfdata(tf(sys),'v'); W=hankel(fliplr(d(1:numel(d)-1))); G=ss(sys); U0=obsv(G); S=W*U0; csys=ss(S*G.a*inv(S),S*G.b,G.c*inv(S),G.d); end
function [ csys ] = ctrbcanon(sys) %CTRBCANON Compute controllability canonical form [n,d]=tfdata(tf(sys),'v'); W=hankel(fliplr(d(1:numel(d)-1))); G=ss(sys); UC=ctrb(G); S=UC*W; csys=ss(inv(S)*G.a*S,inv(S)*G.b,G.c*S,G.d); end
\ifpdf
% pdf code
\else
% dvi code
\fi


>> addxtick(2/pi)とすると

>> addxtick(pi/2,'pi/2')と入れれば

となります.便利ですよ.
function [] = addxtick( tk, tklbl )
%ADDXTICK --- Add new tick to x-axis
% Ex. addxtick(3.14)
% Ex2. addxtick(3.14, 'Pi')
if nargin==1
set(gca,'XTick',unique(sort([get(gca,'XTick'), tk ] )));
else
temp = cellstr(get(gca,'XTickLabel'));
temp = {temp{:},tklbl};
[temp2,idx] = sort([get(gca,'XTick'),tk]);
[temp2,idx2]=unique(temp2);
idx=idx(idx2);
set(gca,'XTick',temp2);
set(gca,'XTickLabel',temp(idx));
end
end
あと,ほとんど同じですがaddytick.mも置いておきます.function [] = addytick( tk, tklbl )
%ADDYTICK --- Add new tick to y-axis
% Ex. addytick(3.14)
% Ex2. addytick(3.14, 'Pi')
if nargin==1
set(gca,'YTick',unique(sort([get(gca,'YTick'), tk ] )));
else
temp = cellstr(get(gca,'YTickLabel'));
temp = {temp{:},tklbl};
[temp2,idx] = sort([get(gca,'YTick'),tk]);
[temp2,idx2]=unique(temp2);
idx=idx(idx2);
set(gca,'YTick',temp2);
set(gca,'YTickLabel',temp(idx));
end
end