function [ hp] = juanji( f,g )
%此函数只为计算两任意二维矩阵的卷积
%命令格式：C=juanji(A,B) C 为 A B 的卷积;
%若 A 为 m*n 矩阵，B 为 p*q 矩阵，则 C 为(m+p-1)*(n+q-1)矩阵;
[m1 n1]=size(f);
[m2 n2]=size(g);
M=m1+m2-1;
N=n1+n2-1;
fp=zeros(M,N);
gp=zeros(M,N);
fp(1:m1,1:n1)=f;
gp(1:m2,1:n2)=g;
Fp=reshape(fp',[M*N 1]); %fp 必须转置 reshape 把原矩阵的各列连接为列矩阵
for i=1:M
G(:,:,i)=gallery('circul',gp(i,:)); %生成循环矩阵
end
temp=gallery('circul',1:M)';
for i=1:M
for j=1:M
a1=(i-1)*N+1;
a2=i*N;
b1=(j-1)*N+1;
b2=j*N;
c=temp(i,j);
Gb(a1:a2,b1:b2)=G(:,:,temp(i,j))';
end
end
Hp=Gb*Fp;
hp=reshape(Hp,[M N]);
end
wmbuaa
- 粉丝: 0
- 资源: 2
我的内容管理
收起
- 我的资源
快来上传第一个资源
我的收益 登录查看自己的收益
我的积分
登录查看自己的积分
我的C币
登录后查看C币余额
我的收藏 
我的下载 
下载帮助
会员权益专享
最新资源
- VMP技术解析：Handle块优化与壳模板初始化
- C++ Primer 第四版更新：现代编程风格与标准库
- 计算机系统基础实验：缓冲区溢出攻击(Lab3)
- 中国结算网上业务平台：证券登记操作详解与常见问题
- FPGA驱动的五子棋博弈系统：加速与创新娱乐体验
- 多旋翼飞行器定点位置控制器设计实验
- 基于流量预测与潮汐效应的动态载频优化策略
- SQL练习：查询分析与高级操作
- 海底数据中心散热优化：从MATLAB到动态模拟
- 移动应用作业：MyDiaryBook - Google Material Design 日记APP
- Linux提权技术详解：从内核漏洞到Sudo配置错误
- 93分钟快速入门 LaTeX：从入门到实践
- 5G测试新挑战与罗德与施瓦茨解决方案
- EAS系统性能优化与故障诊断指南
- Java并发编程：JUC核心概念解析与应用
- 数据结构实验报告：基于不同存储结构的线性表和树实现
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
