结果要这样子的。看起来很复杂的程序，我matlab的水平，直接跪了

或者这样子的结果

你自己贴吧，自己不来啊。混帐！
我没看明白那个机械功率是怎么求的。。。你看明白了吗？
看来还是需要 @sasa_mj

我可是本吧排名第4的元老级人物啊。当心我把你这破吧主给弹劾下来。。。哼哼

已经弄到matlab的文件了。但是。问题是老师要我自己做。别用别人的，弟弟，这玩意太复杂了。要不你帮我弄弄。改下？比如2个文件并一个，或者能弄的简单点啥的，改名字啊啥的。反正看起来像自己弄的，不像别人的就行。还有这个只能弄出3楼的图，怎么搞出2楼的图，你晓得么？ 下面2个是matlab的文件。
http://pan.baidu.com/share/link?shareid=461349&uk=1124266812

这东西从没学过怎么会懂……
matlab顶多改改变量名字，
他名字起的都是有实际意义的，怎么改也是这个样子

改了个pg64.m
其中题目里阻尼0.05，程序里原本是0

看了半天不知道哪个量对应哪个量，不懂这东西没办法

问题就是，这是个电力系统，有64台发电机构成，导线有些阻抗啊，啥的，都是初始值，不变的。然后呢，因为一些随机事件，导致了一些不平衡，然后这个不平衡呢，会对电力系统造成扰动，以波的形式，在64台发电机间传播。所以你看到3楼图，其实是64台发电机在0-100时刻的情况，高出来的地方，说明扰动传播到这里，平的地方就是正常情况，可以很清楚从3楼图中，看到扰动波的传播。像波浪一样。大致理解了吗？ 看看横纵坐标。一个是时间0-100，一个是发电机编号1-64. 我想得到2楼图，就是特殊几个时间内，各台发电机的情况。其实就是3楼图的特殊几个时刻，不知道你能不能通过那2个file，搞出2楼的图呢？ 谢谢

(.5*exp(-.1*(s-16).^2))
只找到这个式子

正是因为不懂公式所以画不出来。
你自己研究下程序画出的图会发现，有64条线，对应64个发动机，每条线有自己的相对波形，合在一起模拟了那个三维的图的外观。
但事实上，每条线的纵坐标并不是实际2楼中的y轴或者3楼中的z轴，也就是那个扰动/弧度，而是经过变换的其他值（64号发动机纵坐标都超过6了，根本不是0.几）
按照t=0s也就是矩阵的第一列来绘制[扰动/弧度-发动机编号]的图，画出来不是2楼的样子，而是这个样子的

1号正常，64号纵坐标已经超过6了
程序目的应该是模拟3楼的图，所以不是真实的扰动/弧度值
不知道你理解我的意思了没
因为不懂公式的意义，只能尝试改其中各个步骤，依然不知道怎么做出真正的扰动/弧度值

擦。弗格森退休了。。。。。。。。
你温格来不来接班？

我现在做了一点修改，
把最后plot那句，改了下，
plot(q,delta_n(q,:)-2*pi*q/n);
然后能出这个图，看起来还不错。

但问题来了。我没法改时间，我在plot前定义事件t=10没用，图像不变化，我把tspan=[0:Td:Tend] 里的0 改成10，20啥的，图像还是不变。。。。

温格的时代已经到来