[已答复] 【绝对原创】Matlab【角脉冲响应处理程序及结果曲线分享】

[复制链接]
shatanshzao 发表于 2010-2-28 16:45:35
不知道大家有没有学汽车动力学的,现把汽车瞬态响应之一的角脉冲响应处理程序和大家分享,大家一起交流,共同提高,附件还有运行后的结果-横摆加速度频率响应特性曲线。程序绝对原创。

脉冲程序和结果.rar

8.52 KB, 下载次数: 12923

10 条回复


cherishkf 发表于 2010-2-28 17:42:51
帮顶一下!!!

shatanshzao 发表于 2010-2-28 23:06:35

为了不让大家花钱,特把程序复制出来

clc
clear
%% 下面的一段是加载方向盘转角和横摆角速度数据,txt数据是两列:第一列是时间,第二列是响应值。
%% 在运行程序之前需要把txt文本数据前面的文字去掉,如果不去掉需加上A=A.data命令,A是某个参数
swa_p=importdata('swa_pulse.txt');  %% 加载方向盘转角数据
t=swa_p(:,1);  %% 加载时间,方向盘转角第一列是时间
swa=swa_p(:,2)*pi/180;
r=importdata('r_pulse.txt');
r=r(:,2)*pi/180;  %% 加载横摆角速度数据,取第二列

A1=0;A2=0;A3=0;A4=0;  %% 中间变量赋初值
dt=t(10)-t(9);  %% 时间步长
f=0:0.01:1.8; %% 频率取到1.8Hz
for kk=1:length(f)
    for k=1:length(t)
        A1=A1+r(k)*cos(2*pi*f(kk)*t(k))*dt;   %% 计算a1,A1就是公式中的a1
        A2=A2+r(k)*sin(2*pi*f(kk)*t(k))*dt;   %% 计算a2,A2就是公式中的a2
        A3=A3+swa(k)*cos(2*pi*f(kk)*t(k))*dt;  %% 计算a3,A3就是公式中的a3
        A4=A4+swa(k)*sin(2*pi*f(kk)*t(k))*dt;  %% 计算a4,A4就是公式中的a4
    end
    B1=(A1*A3+A2*A4)/(A3^2+A4^2);  %% 计算b1,B1就是公式中的b1
    B2=(A1*A4-A2*A3)/(A3^2+A4^2);  %% 计算b2,B2就是公式中的b1
    A_r_f(kk,:)=sqrt(B1^2+B2^2);  %% 幅频特性
    phi_r_f(kk,:)=atan(B2/B1)*180/pi;  %% 相频特性
end
subplot(2,1,1)  %% 绘制相频和幅频特性曲线
plot(f,A_r_f)
xlabel('频率f(Hz)')
ylabel('横摆角速度增益(1/s)')
grid on
subplot(2,1,2)
plot(f,phi_r_f)
xlabel('频率f(Hz)')
ylabel('相位差(deg)')
grid on

figure(2)
subplot(2,1,1)  %% 绘制对数坐标下的幅频和相频特性曲线
A_r_f=20*log10(A_r_f);
semilogx(f,A_r_f);
xlabel('频率f(Hz)')
ylabel('横摆角速度增益(dB)')
grid on
subplot(2,1,2)
semilogx(f,phi_r_f);
xlabel('频率f(Hz)')
ylabel('相位差(deg)')
grid on

mengxg1986 发表于 2012-4-4 18:06:24
谢谢分享!!!

zhongch2010 发表于 2012-7-18 00:34:58
我正在寻找,谢谢分享!

archy_mao 发表于 2013-3-29 15:56:36
楼主您好,我用您的程序进行计算,计算结果感觉有点问题,可否请教一下。我仔细推算了程序,没有问题。我的仿真的结果感觉也没有什么问题,就是频域特性曲线有问题。求教。

pptao 发表于 2014-8-24 10:29:44
用了自己的数据进行处理,幅频和相频特性图还是有问题

dcyan0907 发表于 2015-2-2 16:42:13
看看,谢谢。。。

four_D 发表于 2015-8-6 22:06:51
最近正在做这个

2006zhongch2006 发表于 2017-12-2 10:49:02
感谢分享资料

anka119 发表于 2021-11-24 09:21:59
感谢分享资料:loveliness:
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表