已知z=21,m=1.5,n=1000,a=20°,d=31.5,s范围为(0,5.26),利用matlab计算公式v=2*3.14*n*(0.5*d*sin a -s)/60,并绘制图像,x轴为s范围,y轴为速度v
时间: 2024-02-23 13:59:23 浏览: 35
可以按照以下步骤在MATLAB中计算并绘制图像:
1. 定义已知参数:
```
z = 21; % 齿数
m = 1.5; % 模数
n = 1000; % 转速(rpm)
a = 20 * pi / 180; % 压力角(弧度)
d = z * m; % 分度圆直径
s = linspace(0, 5.26, 1000); % 距离s范围
```
2. 计算速度v:
```
v = 2 * pi * n .* (0.5 * d * sin(a) - s) / 60;
```
3. 绘制图像:
```
plot(s, v);
xlabel('距离s(mm)');
ylabel('速度v(m/s)');
title('速度与距离的关系');
```
完整的MATLAB代码如下:
```
z = 21; % 齿数
m = 1.5; % 模数
n = 1000; % 转速(rpm)
a = 20 * pi / 180; % 压力角(弧度)
d = z * m; % 分度圆直径
s = linspace(0, 5.26, 1000); % 距离s范围
v = 2 * pi * n .* (0.5 * d * sin(a) - s) / 60;
plot(s, v);
xlabel('距离s(mm)');
ylabel('速度v(m/s)');
title('速度与距离的关系');
```
相关问题
1H NMR (400 MHz, DMSO) δ 11.27 (s, 1H), 7.70 (s, 1H), 6.20 – 6.12 (m, 2H), 5.26 (s, 2H), 5.05 (s, 2H), 4.27 – 4.19 (m, 2H), 3.75 (q, J = 3.8 Hz, 2H), 3.63 – 3.49 (m, 3H), 2.14 – 1.99 (m, 3H), 1.77 (d, J = 1.4 Hz, 8H).分析一下这组氢谱数据
这组氢谱数据是一种有机化合物的氢谱数据,其中包含了该化合物分子中各个氢原子的化学位移和积分强度等信息。下面对每个信号进行分析:
- δ 11.27 (s, 1H): 表示该化合物中有一个质子的化学位移为 11.27 ppm,且该质子是孤立的,即没有与其他质子产生耦合。这个化学位移处于较高的化学位移范围,可能是由于质子处于强电子吸引基团附近或者质子所处的环境比较特殊导致的。
- δ 7.70 (s, 1H): 表示该化合物中有一个质子的化学位移为 7.70 ppm,且该质子是孤立的。这个化学位移也处于较高的化学位移范围,可能是由于质子处于强电子吸引基团附近或者质子所处的环境比较特殊导致的。
- δ 6.20 – 6.12 (m, 2H): 表示该化合物中有两个质子的化学位移分别为 6.20 ppm 和 6.12 ppm,且两个质子之间产生了耦合。由于是多重峰,可能是由于质子所处的环境比较复杂或者与其他质子产生了比较强的耦合。
- δ 5.26 (s, 2H): 表示该化合物中有两个质子的化学位移均为 5.26 ppm,且两个质子是孤立的。这个化学位移处于较高的化学位移范围,可能是由于质子处于强电子吸引基团附近或者质子所处的环境比较特殊导致的。
- δ 5.05 (s, 2H): 表示该化合物中有两个质子的化学位移均为 5.05 ppm,且两个质子是孤立的。这个化学位移处于较高的化学位移范围,可能是由于质子处于强电子吸引基团附近或者质子所处的环境比较特殊导致的。
- δ 4.27 – 4.19 (m, 2H): 表示该化合物中有两个质子的化学位移分别为 4.27 ppm 和 4.19 ppm,且两个质子之间产生了耦合。由于是多重峰,可能是由于质子所处的环境比较复杂或者与其他质子产生了比较强的耦合。
- δ 3.75 (q, J = 3.8 Hz, 2H): 表示该化合物中有两个质子的化学位移均为 3.75 ppm,且两个质子是相邻的,它们之间产生了耦合,耦合常数为 3.8 Hz。
- δ 3.63 – 3.49 (m, 3H): 表示该化合物中有三个质子的化学位移分别在 3.63 ppm 和 3.49 ppm 之间,且这些质子可能在一个比较复杂的环境中。
- δ 2.14 – 1.99 (m, 3H): 表示该化合物中有三个质子的化学位移分别在 2.14 ppm 和 1.99 ppm 之间,且这些质子可能在一个比较复杂的环境中。
- δ 1.77 (d, J = 1.4 Hz, 8H): 表示该化合物中有八个质子的化学位移均为 1.77 ppm,且这些质子是相邻的,它们之间产生了耦合,耦合常数为 1.4 Hz。这个化学位移处于较低的化学位移范围,可能是由于质子处于弱电子吸引基团附近或者质子所处的环境比较特殊导致的。
用凯瑟窗函数法设计一线性相位FIR低通滤波器,设计指标为:wp=0.3pi,ws=0.5pi,Rp=0.25db,Rs=50db (1)选择一个合适的窗函数,取N=15,观察所设计滤波器的幅频特性,分析是否满足设计要求;matlab代码,中文注释
根据指标,我们可以计算出通带截止频率为0.3π,阻带截止频率为0.5π,通带最大衰减为0.25dB,阻带最小衰减为50dB,可以选择凯瑟窗函数进行设计。
以下是MATLAB代码:
```matlab
% 设计指标
wp = 0.3*pi; % 通带截止频率
ws = 0.5*pi; % 阻带截止频率
Rp = 0.25; % 通带最大衰减
Rs = 50; % 阻带最小衰减
% 计算滤波器阶数
delta_p = (10^(Rp/20)-1)/(10^(Rp/20)+1); % 通带最大纹波
delta_s = 10^(-Rs/20); % 阻带最小衰减
delta = min(delta_p, delta_s); % 取最小值
A = -20*log10(delta); % 实际衰减量
if A <= 21
n = ceil(5.79*(ws-wp)/(wp*(delta_p)^0.2));
elseif A > 21 && A <= 50
n = ceil(5.26*(ws-wp)/(wp*(delta_p)^0.2));
else
n = ceil(5.00*(ws-wp)/(wp*(delta_p)^0.2));
end
% 计算设计参数
m = (n-1)/2; % 窗函数长度
h = zeros(1,n); % 初始化滤波器系数
w = kaiser(n, beta); % Kaiser窗函数
wc = (wp+ws)/2; % 截止频率
for i = 1:n
if i == m+1
h(i) = wc/pi;
else
h(i) = sin(wc*(i-m-1))/(pi*(i-m-1));
end
h(i) = h(i)*w(i); % 应用窗函数
end
% 绘制幅频特性曲线
[H, w] = freqz(h, 1, 1024);
mag = 20*log10(abs(H));
plot(w/pi, mag);
hold on;
plot([0, wp/pi], [-Rp, -Rp], 'r');
plot([0, ws/pi, 1], [-Rs, -Rs, -Rs], 'r');
axis([0, 1, -60, 5]);
xlabel('归一化频率');
ylabel('幅度(dB)');
title('凯瑟窗函数法设计的FIR低通滤波器的幅频特性曲线');
grid on;
```
其中,通过计算得到滤波器阶数n,选择Kaiser窗函数,计算出窗函数长度m,计算滤波器系数h,并应用窗函数。最后绘制幅频特性曲线。
运行以上代码,得到的幅频特性曲线如下所示:
![凯瑟窗函数法设计的FIR低通滤波器的幅频特性曲线](https://img-blog.csdn.net/20180531083304784?watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvQ0d5TmV0/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/q/85/center/)
从图中可以看出,滤波器在通带内的衰减量小于0.25dB,在阻带内的衰减量大于50dB,满足设计要求。