pifubhu
2008-09-05, 10:42
我在一个已知传函的二阶曲线上取点，拟合出新的二阶曲线，画出他们的波特图，程序已经做出。现在想知道拟合曲线的效果，就是求求新曲线与原曲线的误差。不知道用什么方法好，帮忙看看，谢啦。
下面是我的程序：
clear all
% 建立一个2阶的系统，有理分式形式
origin_den=[1 1];
origin_num=[1 3 1];
w_1=logspace(-1,1);
[fit_num1,fit_den1]=invfreqs(origin_res1,w_1,1,2);
origin_res1=freqs(origin_num,origin_den,w_1);
origin=tf(origin_num,origin_den);
fit1=tf(fit_num1,fit_den1);
% 用波特图来表示拟合的效果
figure(1)
bode(fit1,'b*',w_1)
hold on
%g=tf(origin_num,origin_den);
%bode(g);grid on
bode(origin,'k--',fit1,'r',{0.01,1000})
title('fitting range 0.1 to 10')
grid on
hold off;
下面是我的程序：
clear all
% 建立一个2阶的系统，有理分式形式
origin_den=[1 1];
origin_num=[1 3 1];
w_1=logspace(-1,1);
[fit_num1,fit_den1]=invfreqs(origin_res1,w_1,1,2);
origin_res1=freqs(origin_num,origin_den,w_1);
origin=tf(origin_num,origin_den);
fit1=tf(fit_num1,fit_den1);
% 用波特图来表示拟合的效果
figure(1)
bode(fit1,'b*',w_1)
hold on
%g=tf(origin_num,origin_den);
%bode(g);grid on
bode(origin,'k--',fit1,'r',{0.01,1000})
title('fitting range 0.1 to 10')
grid on
hold off;