clear all; %close all

N=20; %number of points in x-axis

L_Lake=2500; %length of lake (m)

L=10000; %wavelength of atmospheric wave fronts (m)

T=1.25*60^2; %period of atmospheric wave fronts

(seconds)

rho=1000; �nsity of water (kg/m^3)

d_x=L_Lake/N;

d=2; �pth of lake (m)

a=0.01*d; %amplitude of seiche

g=9.81; % Acceleration of gravity

M_n=0; % Manning's n

c=sqrt(g*(d+abs(a)));

alpha=1;% Factor to show effect of exceeding Courant

limit;

d_t=alpha*d_x/c;

Courant=c*d_t/d_x;

M=fix(10000*(1/d_t));

iplot=2; % plots every "iplot" timestep

% Initial conditions;

t=d_t*(0:M-1);

x=-L_Lake/2+d_x*(0:N-1)+d_x/2;% This puts x at eta

points.......

for m=

一维水动力matlab程序,matlab 一维湖泊浅水方程,pressure driven seiche code相关推荐

  1. 影像重采样matlab程序,matlab重采样代码

    本设计通过编 写 matlab 程序,借助对声音音色和音调的双重复合改变,实现输出声音的改变.共振 峰频 率的改变是基 于重采样 实现的,从重 采样原理 知道,这也同...... 本设计通过编 写 m ...

  2. 标号法matlab程序,matlab程序大全

    上例 Matlab 程序如下: %使用表8.4 给出的小波基 ld=[0 0.0... matlab经典代码大全_计算机软件及应用_IT/计算机_专业资料.哈哈哈 MATLAB 显示正炫余炫图:plo ...

  3. 插值拟合matlab程序,Matlab插值和拟合

    数值插值 一维插值 YI=interp1(X,Y,XI,METHOD) X.Y为采样点和样本值,XI为表示要插值的点,METHOD默认是线性插值'linear' METHOD可以取 'nearest' ...

  4. 解抛物型微分方程c语言程序,偏微分方程的matlab代码_偏微分方程matlab程序,matlab偏微分代码...

    matlab偏微分方程含源代码,大家可以下载运行,里面的都能用 Pie Eais ae h t lat BLep 话k气飞"目基回 古班比比,一生热件身的些懂 1身定 回量1 古典显式格式稳 ...

  5. 龙格库塔法解微分方程组的matlab程序,MATLAB实例源码教程:龙格库塔法求解微分方程组源代码实例.doc...

    MATLAB实例源码教程:龙格库塔法求解微分方程组源代码实例.doc MATLAB实例源码教程龙格库塔法求解微分方程组源代码实例题目用经典 Runge-Kutta方法求下列一阶微分方程组的近似解y1 ...

  6. dmc matlab程序,matlab编的DMC程序.doc

    matlab编的DMC程序 clear all; % close all; %系统模型建立 num=[0.8]; den=[225 1]; [a,b,c,d]=tf2ss(num,den); % st ...

  7. 哈密顿圈 MATLAB程序,matlab求最佳哈密顿圈遇到错误

    问题描述: matlab求最佳哈密顿圈遇到错误 求最佳哈密顿圈遇到错误 程序如下 %用矩阵翻转方法来实现二边逐次修**过程,求最佳哈密尔顿圈(H圈) clc clear load('zd.mat'); ...

  8. 约束优化的matlab程序,Matlab多约束最优化问题请高手求助,我编的程序

    问题补充,这是我自己编的, 我刚学matlab没有多长时间,请大家帮忙看下程序有什么不对的地方! x1==csvread('x1.txt');    %x1  是300*10的矩阵 [x1,y,z,v ...

  9. 导出matlab程序,Matlab数据导入导出

    Matlab提供了从磁盘文件或剪贴簿转载数据至工作区(数据导入)和将工作区变量存入磁盘文件(数据导出)的多种途径. 最简单的办法是使用界面导入向导,打开文件菜单中的导入数据而后按提示操作. 一.导入文 ...

最新文章

  1. 使用OpenCV进行直播(附代码)
  2. 站长新手入门:从0开始搭建微信小程序商城,不会代码也能开商城(附带源码)
  3. 金碧辉煌!皇城定制5月22日正式对外运营开业!
  4. 【Ubuntu】Ubuntu16.04安装 搜狗输入法 史上最详细
  5. 车联网服务non-RESTful架构改造实践
  6. 产品经理最痛苦的事情
  7. [转载] python数学计算模块之math常用函数学习使用
  8. Chapter 4 Invitations——18
  9. springcloud架构特点_董事长挖来一位京东T9架构师,送我们两份微服务文档,实在太香了...
  10. leach算法 matlab,求助LEACH的MATLAB仿真代码
  11. 特斯拉开始发布其汽车的开源Linux软件代码
  12. android webview 文件下载,Android编程使用WebView实现文件下载功能的两种方法
  13. 大型连锁百货运维审计用什么软件好?有哪些功能?
  14. VsCode文件屏蔽
  15. JS基础——多张图片循环切换
  16. 全球开发者各出奇招:我们想这样适配iPhone X
  17. kinect for Unity开发
  18. ubuntu 16.04下chrome浏览器Flash过期解决方法
  19. 慕课网 前端JS面试技巧 笔记
  20. 金蝶eas系统服务器怎么启动,EAS应用服务器启动失败

热门文章

  1. 打印棱形--进阶(Java)
  2. 二本计算机专业校招进银行,你是普通二本能够校招进银行吗?
  3. 关于win7 系统屏幕扩展的操作方法
  4. 买就送!你与新年Flag之间只差一个¥888京东卡!
  5. python json数组解析_Python基础系列(一)搞懂json数据解析与字典之间的关系
  6. 风变编程的python课程多少钱-风变编程Python小课最近很火,大家学完感受如何?...
  7. 【图形】交互式绘图软件开发设计
  8. 简单计算器(YZOJ-1047)
  9. Nodejs 快速入门教程
  10. 软件测试职业发展规划,带你摆脱年龄的困扰!不能不看