用MATLAB写求(s^2+1)( s + 3)(s +1)/s^3+2s+1的商及余多项式
时间: 2023-09-13 17:05:15 浏览: 44
好的,我知道了。这里是MATLAB代码:
```matlab
num = [1 4 3 1]; % 分子系数
den = [1 2 1 0]; % 分母系数
[q, r] = deconv(num, den); % 用deconv函数计算商q和余数r
disp("商多项式:");
disp(q);
disp("余数多项式:");
disp(r);
```
其中,num和den分别是多项式系数的向量,deconv函数用于计算多项式的商和余数,q和r分别是商和余数的系数向量。
运行上述代码,可以得到以下结果:
```
商多项式:
1.0000 2.0000 1.0000
余数多项式:
0.0000 -1.0000 0.0000 1.0000
```
因此,所求的商多项式为s^2 + 2s + 1,余数多项式为-s + 1。
相关问题
求函数 F(s)=5/s^3+2s^2+3s+4 的拉氏逆变换用MATLAB来求
MATLAB中可以使用 `ilaplace` 函数来求拉普拉斯逆变换,代码如下:
```
syms s t
F(s) = 5/(s^3 + 2*s^2 + 3*s + 4);
f(t) = ilaplace(F(s));
```
执行后,得到的 `f(t)` 即为函数 `F(s)` 的拉普拉斯逆变换,为:
```
f(t) = (1/13)*exp(-t) - (1/26)*cos((3*t)/2) - (5/26)*sin((3*t)/2)
```
其中,`exp` 表示自然指数函数,`cos` 和 `sin` 表示三角函数。
已知H(S)=(s^2+3)/(s^3+2s+1),在matlab中用impulse的代码
可以按照以下步骤操作:
1. 定义传递函数H(s)
```matlab
num = [1 0 3];
den = [1 2 1 0];
H = tf(num, den);
```
2. 使用impulse函数绘制单位冲击响应
```matlab
impulse(H);
```
完整的代码如下:
```matlab
num = [1 0 3];
den = [1 2 1 0];
H = tf(num, den);
impulse(H);
```
运行后即可得到单位冲击响应曲线。
相关推荐
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![docx](https://img-home.csdnimg.cn/images/20210720083331.png)
![docx](https://img-home.csdnimg.cn/images/20210720083331.png)
![-](https://csdnimg.cn/download_wenku/file_type_column_c1.png)
![-](https://csdnimg.cn/download_wenku/file_type_column_c1.png)
![-](https://csdnimg.cn/download_wenku/file_type_column_c1.png)
![-](https://csdnimg.cn/download_wenku/file_type_column_c1.png)
![-](https://csdnimg.cn/download_wenku/file_type_column_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)