xjzhang 发表于 2009-9-30 11:13

怎么用符号推倒不定积分。

例如我知道传递函数H（S），里面有很多变量（符号退出来的），现在求H（Jw）在整个频域上的积分，积分变量是w，请问怎么办？
顺便问一下我会把传递函数中的分子分母分开，但是怎么能把这个含有很多符号变量的关于s的多详式的系数得到？

sogooda 发表于 2009-9-30 14:26

回复 楼主 xjzhang 的帖子

你提问很笼统，别人很难准确的帮助你。我觉得你还是把问题说更清楚点儿。

xjzhang 发表于 2009-9-30 21:39

请教，这个程序，我怎么用不了啊？
clear
clc
syms ks ms muku C_off C_on alpha s w

A=;
B=';
C=;
D=0;

I=eye(4,4);
G=C*inv(s*I-A)*B+D;
% =ss2tf(A,B,C,D);
% =ss(A,B,C,D);
C2=;
C3=;
G2=C2*inv(s*I-A)*B+D;
G3=C3*inv(s*I-A)*B+D;
output=;
=numden(output);
num=collect(num,'s');
den=collect(den,'s');
我的问题就是把G的分子分母的系数提出来,还有就是
G1=subs(G,'s',j*w);
% Fn=int(G1,w,-inf,inf)
我要对G1求积分，上下限分别是-inf，inf。用int感觉太慢了，怎么办？

ChaChing 发表于 2009-10-1 00:10

感觉上楼主的东东好像就是状态方程求解
考虑是否可直接利用control toolbox现有的函数

xjzhang 发表于 2009-10-1 11:46

我的确是在用状态矩阵推导传递函数H，但是我推导完传递函数后，还要对传递函数|H|^2在整个频域上积分（S——〉jw），
这个我就不会了，都是符号的，我要推公式，controlbox里面好像大多数不支持符号计算阿。

xjzhang 发表于 2009-10-2 04:46

多项式系数问题解决，但是这个积分还是没有搞定，好似H2优化里面经常用到。再次请教阿！

多项式系数如下程序：

syms a b c x
y = a + b*sin(x) + c*sin(2*x)
coeffs(y,sin(x)) =
coeffs(expand(y),sin(x)) =

syms x y
z = 3*x^2*y^2 + 5*x*y^3
coeffs(z) =
coeffs(z,x) =
= coeffs(z,y) returns c = , t =
页: [1]
查看完整版本: 怎么用符号推倒不定积分。