数学建模美赛集训的时候要用到一个海面模拟,分享一下海面模拟的MATLAB代码 先贴一下结果图: 下面是源代码~~~ function waterwave n = 64; % grid size g = 9.8; % gravitational constant dt = 0.01; % hardwired timestep dx = 1.0; dy = 1.0; nplotstep = 8; % plot interval ndrops = 5; % maximum number of drop

颜色和样式 八种内建默认颜色缩写b:blue g:green r:red c:cyan m:magenta y:yellow k:black w:white其它颜色表示方法可以参照百度给的值https://baike.baidu.com/item/%E5%8D%81%E5%85%AD%E8%BF%9B%E5%88%B6%E9%A2%9C%E8%89%B2%E7%A0%81/10894232?fromtitle=%E9%A2%9C%E8%89%B2%E5%8D%81%E5%85%AD%E8%BF%

MATLAB模拟布丰投针实验 标签(空格分隔): 算法 Buffon's Needle 桌面上有距离为a的若干平行线,将长度为L的针随机丢在桌面上,则这根针与平行线相交的概率是多少?假定L < a. 思路:从针据横线的距离与夹角得出. 解决: 假设针的中点到最近横线的距离为y,则\(y\in[0,\frac{a}{2}]\); 因为投针是随机的,所以y服从均匀分布: \[ f(y) = \begin{cases} \frac{2}{a}, & \text{$0 \leq y \leq \fr

多普勒失真信号采样Matlab模拟分析 方案 水声通信指的是使用声信号在水中数据传输. 相对而言.电磁信号在水中吸收严重衰减过快,光信号受水中悬浮颗粒的影响,也无法完毕远距离传输. 这两种信号的传播距离约为数百米,而声信号在水中传播距离达到数十千米.成为水下通信的首选方式. 水声通信的一个重大挑战就是运动带来的多普勒. 多普勒造成信号的频率漂移,影响系统性能,因此很有必要消除多普勒带来的影响.而多普勒在时域的表现为信号的长度压缩和扩展.即'y(t) = x((1+a) * t)',当中a为多普勒

还记得原来写过一个对图像进行波纹扭曲操作的博文. 这次实现的是水波特效,其实就是通过正余弦函数表示波纹中心位置慢慢向外扩散,通过叠加衰减因子使振幅不断减小,进而产生水波的效果. 效果如下: 原图: 波纹特效: matlab代码如下:   clear all; close all; clc; cenH=128; %波纹中心 cenW=128; R=80; %扩散半径 amplitude =3; %幅度 wavelength=30; %波长 phase =20; %相位 img = double(i

2012-06-08 21:26:42 用matlab来建模,仿真不同时刻os task在队列中的装载情况.输入参数如下 作为初学者,M文件写的有点长.能实现功能就算学以致用了. clear;clc %os每隔0.5ms进入一次中断,进行taskpush.设置初始值,间隔值,及结束时间. StartTime=0;      %第0ms开始push task入队列. TimeInterval=0.5; %每隔0.5ms push一次task. EndTime=100;       %计算到第10m

.katex { display: block; text-align: center; white-space: nowrap; } .katex-display > .katex > .katex-html { display: block; } .katex-display > .katex > .katex-html > .tag { position: absolute; right: 0px; } .katex { font-style: normal; font

1.初版,边缘未处理. % % In----near % If----far % In=imread('D:\文件及下载相关\桌面\模拟变焦拼接\Matlab_code\nearframe\frame1.bmp');%近图 If=imread('D:\文件及下载相关\桌面\模拟变焦拼接\Matlab_code\farframe\frame1.bmp');%远图 If=rgb2gray(If); In=rgb2gray(In); I_temp=If; I_need=zeros(,); I_midd

原文地址 https://blog.csdn.net/lircsszz/article/details/80249017 最近在研究图像校正,现将镜头失真中常见的径向畸变(radial distortion)和切向畸变(tangential distortion)使用Matlab模拟出来. 镜头畸变的数学模型如下: 将此模型进一步简化为两部分 由于opencv中畸变系数只与k1,k2,k3,p1,p2有关,所以在此我也只考虑这五个系数对畸变的影响. 使用矢量场绘制径向畸变: function r

1. 卷积编码增益性能.BER 信道环境:AWGN 信噪比SNR :0:0.1:6 MATALB仿真架构:源比特 +卷积码 +BPSK +AWGN +Viterbi +BER 说明:卷积编码,不同的R码率.不同的约束长度 2. Turbo增益性能.误码率 信道环境:AWGN 信噪比SNR :-1:0.1:2 MATALB仿真架构:源比特 +Turbo编码 +BPSK +AWGN +Turbo译码+BER 说明:turbo译码.不同的译码算法那,不同的交织长度,不同的迭代次数. LTE 标准的tu

本文是利用蒙特卡罗算法对马氏链过程的模拟.假设有10个状态,从每个状态到与之相邻状态的概率是相同的,仿真次数为1000,及进行了1000次状态转移.我们以动画的形式再现了状态转移的过程,并记录了到达每个状态的次数,具体实现如下: close all;clc;clear; figure; s=1; n=1000; r=1; % 圆圈的半径 title('等概率情况的计算机模拟') set(gcf,'doublebuffer','on'); % 设置图形渲染效果 xlabel('Please pre

  %%%% Tutorial on the basic structure of using a planar decision boundary %%%% to divide a collection of data-points into two classes. %%%% by Rajeev Raizada, Jan.2010 %%%% %%%% Please mail any comments or suggestions to: raizada at cornell dot edu

原文地址:http://blog.sina.com.cn/s/blog_79ecf6980100vcrf.html 低通滤波器参数:Fs=8000,fp=2500,fs=3500,Rp=1dB,As=30dB,其他滤波器可以通过与低通之间的映射关系实现. %%模拟滤波器 %巴特沃斯——滤波器设计 wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30; [N,wc]=buttord(wp,ws,Rp,As,'s')%计算率波器的阶数和3dB截止频率 [B,A]=butter(N

低通滤波器参数:Fs=8000,fp=2500,fs=3500,Rp=1dB,As=30dB,其他滤波器可以通过与低通之间的映射关系实现. %%模拟滤波器 %巴特沃斯——滤波器设计 wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30; [N,wc]=buttord(wp,ws,Rp,As,'s')%计算率波器的阶数和3dB截止频率 [B,A]=butter(N,wc,'s');%计算滤波器系统函数分子分母多项式 fk=0:800/512:8000;wk=2*pi*fk; H

理论基础:大数定理,当频数足够多时,频率可以逼近概率,从而依靠概率与$\pi$的关系,求出$\pi$ 所以,rand在Monte Carlo中是必不可少的,必须保证测试数据的随机性. 用蒙特卡洛方法进行计算机模拟的步骤:[1] 设计一个逻辑框图,即模拟模型．[2] 根据流程图编写程序,模拟随机现象．可通过具有各种概率分布的模拟随机数来模拟随机现象．[3] 分析模拟结果,计算所需要结果. ex1.投针试验求$\pi$ %蒲丰投针实验的计算机模拟 format long; %设置15位显示精度 a=

在无风情况下的喷泉模拟 我的python代码 import numpy as np import random import matplotlib matplotlib.rcParams['font.sans-serif']=[u'simHei'] matplotlib.rcParams['axes.unicode_minus']=False import matplotlib.pyplot as plt tt = np.linspace(0,4,1000) #y = -tt**2+4*tt #

大家好,我是小鸭酱,博客地址为:http://www.cnblogs.com/xiaoyajiang %单服务员排队模型模拟100天 clear clc day = 100 ;s = zeros(1,day) ;wait =zeros(1,day) ; for i = 1 :day %首先,赋初值 tg = exprnd(10) ;     %第一个顾客到达时间 ts = tg ;             %售货员已经工作的时间,这里的初值是等待第一个顾客到达的时间 nd = 1 ;      

Github地址:https://github.com/nuptboyzhb/WaterWaveView 欢迎Fork.欢迎Star 1.先看效果 watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvbnVwdDEyMzQ1Njc4OQ==/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/Center" alt="" /> 2.再看关键代码 描绘函数y

1. 一维信号 构造离散时间向量: Fs = 1000; % sampling frequency,采样频率 T = 1/Fs; % sampling period,采样周期 L = 1000; % Length of Signal t = (0:L-1)*T; % time vector 构造信号 S = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t); % 信号包含一个振幅为 0.7 频率为 50 hz 的正弦信号, % 振幅为 1 频率为 120 的正弦信号 加噪

http://www.matlabsky.com/thread-43937-1-1.html   <量化投资:以MATLAB为工具>连载(3)基础篇-N分钟学会MATLAB(下)     <量化投资:以MATLAB为工具>简介          <量化投资:以MATLAB为工具>是由电子工业出版社(PHEI)下属旗舰级子公司——北京博文视点资讯有限公司出版的<量化投资与对冲基金丛书>之一,丛书主编为丁鹏博士,<量化投资:以MATLAB为工具>由李