完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
扫一扫,分享给好友
%求解光顺因子
m=6; g=[0.38 0.58 0.5 0.23 0.37 0.41 ] I=eye(m) s=0:0.1:1000; A=[0.0205227 6.89323 38.6113 71.0299 70.7668 73.087 148.546 247.548 291.872 294.158 331.071 485.808 567.495 565.676 634.502 820.299 989.925 981.243 1080.96 1285.07 0.0623642 10.6999 46.1157 56.019 56.4428 144.17 197.936 202.014 261.283 386.112 461.67 453.307 583.286 794.456 757.685 829.912 1098.27 1206.04 1179.94 1372.7 0.0712628 13.3598 51.1136 55.5443 76.9179 169.655 231.139 206.749 320.06 453.31 468.114 520.922 726.649 818.543 821.85 1070.33 1253.7 1313.94 1354.82 1695.94 0.0807546 15.1393 57.9217 62.9426 87.163 192.253 261.925 234.287 362.691 513.689 530.465 590.306 823.435 927.569 931.315 1212.89 1420.68 1488.95 1535.27 1921.83 0.0908397 17.0299 65.1552 70.8031 98.0483 216.262 294.636 263.546 407.985 577.841 596.712 664.027 926.27 1043.41 1047.62 1364.36 1598.1 1674.9 1727.01 2161.84 0.11279 21.1449 80.8989 87.9115 121.74 268.518 365.83 327.227 506.568 717.466 740.897 824.478 1150.09 1295.53 1300.76 1694.03 1984.26 2079.61 2144.31 2684.21 ] H=[ 1 -2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -2 5 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 6 -4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -4 5 -2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -2 1 ]; K=A*inv(A'*A+s*H)*A'; T=trace(I-K); T1=(I-K)*(g'); n1=norm(T1,2); V=m*n1^2/T^2; s1=V(V+1); plot(s,s1) 麻烦大家帮我看一下,这个程序有什么问题,一直显示倒数第七行有错?感激不尽!!!A是一个6X20的矩阵,H是20X20矩阵
|
|
相关推荐
6个回答
|
|
不一定是倒数第七行 建议从头看是检查
|
|
|
|
调试下吧 建议把程序拷回MATLAB,一行一行的调试,找问题,我一般是这样。。。
|
|
|
|
提示的是啥错,都舍不得给大家讲下,怎么给你看呢?
|
|
|
|
我试了一下 出现的是下面的问题:
K=A*inv(A'*A+s*H)*A'; 括号里面的 s*H 两个矩阵的行列不匹配,因为s(1*10001),H(20*20) 出现下面的错误。 ??? Error using ==> times Matrix dimensions must agree. Error in ==> Untitled at 33 K=A*inv(A'*A+s.*H)*A'; [img]file:///C:Documents and SettingsAdministratorApplication DataTencentUsers459843989QQWinTempRichOle9[)Z3`C0%F4)3I]X6@TTFF5.jpg[/img] 评分
|
||
|
||
我调试了一下 你也可以调试一下 是你的括号里面的 s*H 两个矩阵的行列不匹配
|
|
|
|
{:1:}{:1:}
|
|
|
|
你正在撰写答案
如果你是对答案或其他答案精选点评或询问,请使用“评论”功能。
1622 浏览 1 评论
200圆!求助大佬给一份VSG并网和离网模式的simulink仿真
1900 浏览 0 评论
MATLAB(3)--矩阵的引用(sub2ind、ind2sub、reshape函数使用)
2720 浏览 0 评论
3055 浏览 0 评论
4095 浏览 1 评论
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2024-11-25 09:46 , Processed in 0.923360 second(s), Total 93, Slave 71 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191 工商网监 湘ICP备2023018690号