参考:http://www.ilovematlab.cn/thread-101148-1-1.html
这个ok:
clc;clear; h = [10 9 8 7 6 5 4 3 2 1];
size=length(h);
t=zeros(1,size);
t(1)=h(1);
t(1,2:size)=h(size:-1:2);
H=toeplitz(h,t)
这个也ok:
clc;clear;
h = [10 9 8 7 6 5 4 3 2 1];
size = length(h);
newh = zeros(size,size);
newh(1,1) = h(1);
newh(1,2:size) = h(size:-1:2);
for i = 2:size
newh(i,1) = newh(i-1,size);
newh(i,2:size) = newh(i-1,1:9);
end
这个也ok:
clear;
clc;
a = [1 1 1 2 2 2 3 3 3 4 4 4 5 5 5 6 6 6 7 7 7 8 8 8];
h = [a;a;a];
%function producematrix(h,m)
[height width] = size(h);
m = 3;
c = cell(1,width/m);
for i = 1:m:width
c{(i+m-1)/m} = b(1:m,i:i+m-1);
end
newc = cell(length(c),length(c));
newc{1,1} = c{1};
for i = 2:length(c)
newc{1,i} = c{length(c)-i+2};
end
for i = 2:length(c)
newc{i,1} = newc{i-1,length(c)};
for j = 2:length(c)
newc{i,j} = newc{i-1,j-1};
end
end
%h和m是你的输入,m是h中构成方阵的大小;
对角化:
[V, D]= eig(newh) D=inv(V)*newh*V
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。
发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/184947.html原文链接:https://javaforall.cn
扫码关注腾讯云开发者
领取腾讯云代金券
Copyright © 2013 - 2025 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号:粤B2-20090059 深公网安备号 44030502008569
腾讯云计算(北京)有限责任公司 京ICP证150476号 | 京ICP备11018762号 | 京公网安备号11010802020287
Copyright © 2013 - 2025 Tencent Cloud.
All Rights Reserved. 腾讯云 版权所有