exm040401.m

来自「matlab电子书籍」· M 代码 · 共 11 行

format rat					
p1=conv([1,0,2],conv([1,4],[1,1]));		
p2=[1 0 1 1];					
[q,r]=deconv(p1,p2);		
cq='商多项式为  '; cr='余多项式为  ';
disp([cq,poly2str(q,'s')]),disp([cr,poly2str(r,'s')])

qp2=conv(q,p2);				
pp1=qp2+[zeros(1,length(qp2)-length(r)),r];
pp1==p1