matlab画三维图像例题_matlab画三维曲面

matlab画三维图像例题_matlab画三维曲面我可能真的太闲了吧···在网上看到这图觉得挺好看:t=0:pi/360:2*pi;x=sin(t);y=cos(t);z=2*x.^2+y.^2;figureplot3(x,y,z,’Color’,’r’,’LineWidth’,2);xlabel(‘x’);ylabel(‘y’);zlabel(‘z’);title(‘z=(cost)^2+2*(sint)^2’);ax…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE稳定放心使用

二维函数的图像必须用三维坐标系画,很多函数的图像长得很好看,很有意思。会可视化也会有利于我们理解更深刻。

(1)plot3,三维曲线图

在网上看到这图觉得挺好看:

t=0:pi/360:2*pi;
x=sin(t);
y=cos(t);
z=2*x.^2+y.^2;
figure
plot3(x,y,z,'Color','r','LineWidth',2);
xlabel('x');
ylabel('y');
zlabel('z');
title('z=(cos t)^2+2*(sin t)^2');
axis([-1.2 1.2 -1.2 1.2 1 2.2]);

在这里插入图片描述

(2)mesh,网面图

clear all,clc;
t=-0.5:.01:0.5;
[x,y]=meshgrid(t);%形成格点矩阵
z=sin(4*pi*x)+cos(6*pi*y);
figure(1)
mesh(x,y,z);
axis([-0.5 0.5 -0.5 0.5 -2 2]);
title('z=sin(4*pi*x)+cos(6*pi*y); mesh')
colormap cool%cool是一种配色方案,还有其他方案如winter,summer····见help colormap
colorbar

在这里插入图片描述

画个二元函数,观察采样后频谱的混叠现象

在这里插入图片描述

二元函数图像的绘制还是把变量分为多个离散点绘制,只是增大点数。

clear all,clc;
t=-1:.01:1;
[x,y]=meshgrid(t);%形成格点矩阵
f=cos(2*pi*(2*x-y));
figure(1)
mesh(x,y,f);grid on
xlabel('x');
ylabel('y');
zlabel('f');
axis([-0.5 0.5 -1 1 -1.5 1.5]);
title('f=cos[2*pi(2x-y)]; mesh')
colormap winter
colorbar

在这里插入图片描述
先手动计算一下它的二位傅里叶变换:
不想打公式(逃)

在这里插入图片描述


F= fft2(f);%做二维DFT变换,把连续函数取样得到一个矩阵,把这个矩阵当图像一样处理
figure(2)
subplot(121),imshow(abs(F),[]);title('函数f=cos(2*pi(2x-y))二维频谱')
F=fftshift(F);%使用matlab的fftshift函数把低频分量移到频谱中心
subplot(122),imshow(abs(F),[]);title('中心化频谱')

在这里插入图片描述

t=-1.1:.333:1.1;%以频率3采样(题目要求)
[x,y]=meshgrid(t);%形成格点矩阵
f=cos(2*pi*(2*x-y));
figure(1)
mesh(x,y,f);grid on
xlabel('x');
ylabel('y');
zlabel('f');
axis([-1 1 -1 1 -1.5 1.5]);
title('采样后的f=cos[2*pi(2x-y)]; mesh')
colormap winter
colorbar

F= fft2(f);%做二维DFT变换,把连续函数取样得到一个矩阵,把这个矩阵当图像一样处理
figure(2)
subplot(121),imshow(abs(F),[]);title('函数f=cos(2*pi(2x-y))采样后的二维频谱')
F=fftshift(F);%使用matlab的fftshift函数把低频分量移到频谱中心
subplot(122),imshow(abs(F),[]);title('中心化频谱')

采样后信号:
在这里插入图片描述
采样后信号频谱:
函数的x方向频率是2,y方向频率是1,所以一定会有混叠失真
可以看出频谱的能量中心的相对位置都不对,这就是混叠失真引起的

在这里插入图片描述

t=-1.1:.17:1.1;%以更高频率采样
[x,y]=meshgrid(t);%形成格点矩阵
f=cos(2*pi*(2*x-y));
figure(1)
mesh(x,y,f);grid on
xlabel('x');
ylabel('y');
zlabel('f');
axis([-1 1 -1 1 -1.5 1.5]);
title('采样后的f=cos[2*pi(2x-y)]; mesh')
colormap winter
colorbar

F= fft2(f);%做二维DFT变换,把连续函数取样得到一个矩阵,把这个矩阵当图像一样处理
figure(2)
subplot(121),imshow(abs(F),[]);title('函数f=cos(2*pi(2x-y))采样后的二维频谱')
F=fftshift(F);%使用matlab的fftshift函数把低频分量移到频谱中心
subplot(122),imshow(abs(F),[]);title('中心化频谱')

以更大的采样频率采样,避免不满足奈奎斯特定理产生混叠:
可见得到的频谱和连续函数的频谱更相似。采样频率越大,得到的频谱和连续函数的频谱越相似
在这里插入图片描述

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/183032.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • 【转载】一些VS2013的使用技巧

    【转载】一些VS2013的使用技巧

    2021年11月18日
  • 证书注册_农行提示证书未申请

    证书注册_农行提示证书未申请https证书是互联网传输中很重要的数字证书,为自己的网站申请https证书似乎成为了一件十分有必要的事情,https证书不仅仅可以很好地维护访问网站的用户信息安全,还能为网站树立一个安全可信的形象,维持一个网址的信誉。为https申请证书也变成了一件迫在眉睫的任务,但是如何进行https申请证书?https证书用处的原理是怎么样的?​如何进行https申请证书?这是很多刚接触互联网的网友会带有的疑问。但是我们需要先了解一下何为https证书。https证书是一种有http网络协议和ssl证书.

  • 查看gcc的版本

    查看gcc的版本查看gcc版本 命令:debian:dpkg-lgccredhat:rpm-qa|grepgcc

  • CSRF攻击与防御(写得非常好)「建议收藏」

    转载地址:http://www.phpddt.com/reprint/csrf.htmlCSRF概念:CSRF跨站点请求伪造(Cross—SiteRequestForgery),跟XSS攻击一样,存在巨大的危害性,你可以这样来理解:攻击者盗用了你的身份,以你的名义发送恶意请求,对服务器来说这个请求是完全合法的,但是却完成了攻击者所期望的一个操作,

  • bootstrap 合并单元格之mergeCells属性合并

    bootstrap 合并单元格之mergeCells属性合并bootstrap合并单元格之mergeCells属性合并合并单元格有多种实现方式本文是根据bootstrap自带的mergeCells属性实现的单元格合并,原理是根据有规律的排序数据然后在展示层根据数据行数合并,这样的劣势是数据一定要有规律事先要排序还有一种是比较灵活的是事先处理好数据合并成一行,然后自己去展示层设计怎么展示,而且不用考虑分页问题,详情在我的…

  • html一个汉字空格占位_html空格字符

    html一个汉字空格占位_html空格字符1. (常用)不换行空格,全称No-BreakSpace,它是按下space键产生的空格。空格不会累加(只显示一个)。使用html表示才会累加,该空格占据宽度受字体影响。2. 半角空格,全称EnSpace,en为em宽度的一半(em类似于px受设置不同为20px=1em或其他自定义大小)。占据0.5个中文宽度,不受字体影响。3、&em…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号