景兰-优快云博客

原创 MATLAB画图详细教程

本文将详细介绍如何用matlab绘图并美化。关于figure()创建图窗窗口：figure()figure()的属性：'Name'：在标题栏显示的名称，接字符串，如'Test''Position'：在电脑屏幕上的位置和大小，后接向量[left,bottom,width,height]也就是说指定了图窗的左下角位置，再向右+width、向上+height延伸。'unit'：单位，可以选择'normalized'，使得位置参数为屏幕的归一化位置，屏幕归一化坐标(x,y)，(0,0)表示左

2020-12-10 21:14:57 17124

原创 MATLAB并行实现的简单方法

此方法只是利用了matlab的设定，不需要额外知识。众所周知，matlab是单线程的，但matlab的每个应用窗口都是一个线程，因此可以同时开启多个MATLAB跑程序，占满CPU的所有core。例如处理多个文件的计算，单线程会很慢，所占内存却不多，这时可以写4个m文件，在四个MATLAB应用里同时跑，使得速度翻倍。每个MATLAB应用里的变量都是独立的，因此计算完全不会互相影响。这就类似于开多个窗口刷视频播放。如果并行需要简单的交互，可以写一个通信文件，让所有的M文件都读写这个文件，这样互相之间

2021-03-19 16:00:07 10915 5

转载【SAC】笔记

读写：https://seisman.github.io/SAC_Docs_zh/commands/write/保存文件：https://seisman.github.io/SAC_Docs_zh/graphics/save-image/

2019-09-03 10:52:02 359

原创【linux】设置镜像源

首先找到默认的源$cd /etcapt$cp source.list source.list_default然后替换掉source.list文件为国内源。清华源

2019-09-02 16:34:45 614

原创【matlab】零相位延迟滤波器

使用designfilt命令自动生成滤波器，然后my_filt=‘生成的滤波器代码’，最后filtfilt(my_filt,signal)

2019-09-01 14:53:46 3052

原创【服务器】sbatch 提交作业脚本

在已经写好的脚本的#!bin/bash下面加上：#SBATCH --get-user-env#SBATCH --mail-type=end#SBATCH -J simple_module#SBATCH --nodes=4#SBATCH --ntasks-per-node=24```然后执行脚本的时候用sbatch ./name_of_bash即可...

2019-07-08 16:07:12 23922

原创【matlab】画图的文字调整大小

h=xlabel('x');set(h,'Fontsize',14);h=ylabel('df(x)/dx');set(h,'Fontsize',14);h=title('精确解和二阶差分对比');set(h,'Fontsize',14);h=legend('精确的一阶偏导','二阶差分得到的偏导');set(h,'Fontsize',14);这样就行了...

2019-05-06 00:38:38 2400

原创【Matlab】滤波器常用命令

在命令行中输入：designfilt然后就会弹出滤波器的种类，选定后会让你填参数。参数填好之后点确定，会在命令行窗口生成这个滤波器，复制到编辑器里直接用就好。%读取音频文件%[x,Fs]=audioread('C:\Users\16000\Desktop\testfile.wav');%x是一个序列%Fs是这个音频的采样率%复制过来的滤波器%High_Filt=designfi...

2018-12-10 17:19:57 2372

原创 [matlab]斜阶梯式长条引起的重力异常计算

画的是《重力学与固体潮》那本书的斜阶梯的重力异常function [gravity]=draw_gravity_dixing(h1,h2,alpha,xigma,x0,z,direction)%h1=100;%上表面%h2=200;%下表面%x0=0;%移动原点位置，向右为正%xigma=3*10^3;%剩余密度pi=3.1415926535%alpha=alpha*pi;%倾角...

2018-12-03 13:02:06 1097

原创立方体引起的引力异常计算＆画图

地球重力学需要我们计算立方体引起的重力异常，公式见《重力学与固体潮》。这个程序取的z方向是竖直向下的，也就是说地面向下为正，地面向上为负%定义一个立方体函数%function [gravity]=draw_square(a,b,c,x0,y0,H,ph,z)%长方体模型参数说明%%a=2000;%长%b=200;%宽%c=100;%高%%质心坐标x0,y0,z0 %H=1000立方体深...

2018-12-03 11:54:16 993

原创以二维振动为例展示使用matlab画圆形三维图的偷懒方法

画图主要是用极坐标，我们使用画复数的函数来画。clc;%清空上次的%draw_Bessel%画这个东西a=4;%圆形半径%%定义角向和径向的步长和范围fai=[-4:0.1:4];r=[0:0.02:a];%将fai和r组合起来并生成复平面[R,T] = meshgrid(fai,r);Z = R.*exp(1i*T);%和画具体图有关的m=1;%阶数%k_11=4.6...

2018-12-03 11:40:45 1652

转载使用matlab播放特定频率的声音

A=2;%振幅f_0=397%声音频率fs=10000; %采样频率N=3000; % 信号样点数，播放时长y=A*sin(2*pi*f_0*(0:N-1)/fs); %单频信号sound(y,fs); %通过声卡放音

2018-12-03 11:23:31 9809 3

原创 Matlab自带的Bessel函数及如何扩展其应用

第一类虚宗量Bessel函数I_v:help besseli第二类虚宗量Bessel函数K_v:help besselk第一类Bessel函数J_vhelp besselj第二类Bessel函数 Y_v(诺伊曼函数N_v)help bessely汉克尔函数H_vhelp besselh以上函数的问题是输出结果都是实数。如果想算复数域，以上函数都有级数表达式，用级数求和...

2018-11-28 09:40:17 3358

原创【用word快速输入公式】希腊字母名称和大小写

用word按住alt和加号，先输入\ 然后后面跟代码，输完以后用空格转换。输入大写就首字母大写，小写就首字母小写。代码基本就是发音。大写小写发音AαAlphaBβBetaΓγGammaΔδDeltaΕεEpsilonΖζZetaΝνNuΞξXiOοOmicronΠπPi...

2018-11-28 09:40:01 14709

原创 [Matlab]函数表达式运算

MATLAB对函数表达式进行求和有时候要进行函数表达式的计算syms m; %m是符号，声明符号f1=exp(-im*1);f2=exp(-im*2);f=f1+f2;%可以以符号的形式存储下来%计算%subs(f,'m',1);%将f中的m替换为1总结一下，syms用于声明符号，subs用于将符号替换为数字，进行计算。syms了以后，还有很多函数用于求导积分什么的。详情可...

2018-11-21 20:41:16 12190

qq_37486958的博客