Reckonist-CSDN博客

原创 Matlab绘制双波干涉

直接上图**源波****双波叠加****干涉条纹**代码%波的干涉 clcclear allf=10;%波源频率(两波同频)A=2;%幅值x=-5:0.01:5;y=-5:0.01:5;L_x=length(x);z=zeros(1,L_x);[X1,Y1,Z1]=griddata(x,y,z,linspace(min(x),max(x),150)',... linspace(min(y),max(y),150),'cubic');%插值[X2,Y2,Z2]=g

2021-05-17 10:40:27 3792 2

原创 Makefile出错 make: *** No rule to make target ‘.c‘，needed by ‘.o’.Stop

今天在写makefile的时候遇到了这个问题然后搜索了一下，非常多的相关文章描述了这个问题，指出问题的可能原因是在makefile中多了空格遂逐行检查，反复几遍都没有发现问题（烦躁ing）想了想是不是文件名有错误，又来回对了几遍然后扫了一眼文件夹，发现文件的名字没有对齐，那个空格竟然是在文件名里（可算找到了），再次make正常所以遇到这个问题时候可能原因：1）makefile 文件中出现多余空格2）可能是文件名出错导致（我这里是多了一个末尾空格）...

2020-12-03 12:40:41 5425 1

原创 while循环中scanf()函数被跳过的问题

while循环中scanf()函数被跳过的问题在写一个循环输入时遇到的情况，在一个scanf()之后第二个scanf()被跳过；于是写了段测试代码找找问题。#include "stdio.h"int main(){ char s; while(1) { printf("enter the value of s\n"); scanf("%c", &s); printf("thr value of s is %c\n", s); } return 0;}上面这段代

2020-11-07 15:56:53 2326 2

原创 Matlab 编程《计算流体力学基础及应用（约翰D安德森）》全亚声速等熵喷管流动CFD解法拉瓦尔喷管非守恒形式方程解法

Matlab 编程《计算流体力学基础及应用（约翰D安德森）》全亚声速等熵喷管流动CFD解法拉瓦尔喷管非守恒形式方程解法问题之全亚声速等熵喷管流动的CFD解法初始化参数迭代过程绘图结果问题之全亚声速等熵喷管流动的CFD解法问题的提法见《计算流体力学基础及应用（约翰D安德森）》中的7.4节，下面直接给出写的程序。初始化参数%全亚声速的喷管流动非守恒形式方程解法%MacCormack方法预估-校正具有二阶精度%初始参数（无量纲）clc;clear;clfx1=0:0.1:1.5;

2020-05-12 20:45:23 5038

原创 Matlab 编程《计算流体力学基础及应用（约翰D安德森）》亚声速-超声速等熵喷管流动CFD解法拉瓦尔喷管守恒形式方程解法

《计算流体力学基础及应用（约翰D安德森）》亚声速-超声速等熵喷管流动CFD解法，拉瓦尔喷管守恒形式方程解法 Matlab 编程问题之亚声速一超声速等熵喷管流动的CFD解法初始化参数迭代过程绘图结果最近学习经典的计算流体力学入门书籍《计算流体力学基础及应用（约翰D安德森）》，跟着作者给出的过程写了计算程序，程序写的比较初级。问题之亚声速一超声速等熵喷管流动的CFD解法问题的提法见《...

2019-12-17 16:32:43 5732 4

原创 Matlab广义追赶法（Thomas法）

矩阵分解方法（广义追赶法）对于一般形式的线性方程组Ax=bA x=bAx=b将系数矩阵A分解成下三角阵L和上三角阵UA=LU1A=LU_1A=LU1化归为两个三角方程组Ly=bLy=bLy=b，Ux=yUx=yUx=y对于3阶矩阵：A=[a11a12a13a11a12a13a11a12a13]=[l1100l21l220l31l32l33][1u12u1301u23001]=LU...

2019-06-21 11:28:46 4744

原创 Matlab 模拟"行波"动画

在[0,4π]区间，根据y(t,x)=e−0.2xsin(πt24−x)y(t,x)=e^{-0.2x}sin(\frac{πt}{24}-x)y(t,x)=e−0.2xsin(24πt−x)表现“行波”matlab程序：t=pi*(0:0.01:20);%时间控制l=length(t);%获取时间数组长度x=pi*(0:0.01:4);%[0,4π]区间y=exp(-0.2*x)....

2019-05-26 12:31:19 5829

原创 Matlab 数值计算遇到的问题 linspace与定步长冒号赋值

有没有大佬解答一下定步长冒号赋值存在的问题问题如下：二阶线性系统归一化冲激响应clf,clc,clear;t=0:0.1:18;zt=0.2:0.2:1.4;%阻尼系数m=length(zt);n=length(t);y=zeros(1,n);%初始化hold onfor k=1:m bt=sqrt(abs(1-zt(k).^2)); if zt(k)<1...

2019-05-24 11:32:12 1373

转载 Matlab condeig函数与eig函数

特征值、特征向量及与特征值有关的条件数eig语法形式e = eig(A)[V,D] = eig(A)[V,D,W] = eig(A)e = eig(A,B)[V,D] = eig(A,B)[V,D,W] = eig(A,B)[___] = eig(A,balanceOption)[___] = eig(A,B,algorithm)[___] = eig(___,eigvalOp...

2019-05-20 21:53:53 1713

原创 matlab的有限元编程练习 4nodes

matlab的有限元编程练习求各个节点的位移主函数%单元节点编号clear allELEM=[1 2 5 4; 2 3 6 5; 4 5 8 7; 5 6 9 8];global C;%节点坐标矩阵x1=-1;y1=0;x2=0;y2=0;x3=1;y3=0;x4=-1;y4=1;x5=0;y5=1;x6=1;y6=1;x7=-1;y...

2019-05-16 15:07:29 9511 16

马儿不吃草