tic; clear clc M=[, ];%空间步数 N=*M;%时间步数 :length(M) h=/M(k);%空间步长 tau=/N(k);%时间步长 s=tau/h;%步长比 x=:h:; t=:tau:; y=inline('exp(x+t)','x','t');%真解函数 :length(x) :length(t) exact(i,j)=y(x(i),t(j));%真解 end end u=zeros(M(k)+,N(k)+);%数值解内存单元 :M(k) u(i,)=exp(x(…

双曲圆锥螺线 #http://www.mathcurve.com/courbes3d/spiralehyperbolique/spiralehyperbolique.shtml vertices = t = *PI) to (*PI) r = a = rand2(PI*0.2, PI*0.8) x = r*cos(t)/t z = r*sin(t)/t y = r*cot(a)/t 看其脚本代码与截图,当t为无穷大时,上下两端的曲线才会相交.可以逆向思维一下,将t取倒数,脚本代码变为: #htt…

syms x y z=exp(2*x+y)+cos(3*x*y)-exp(1)-1; zz=subs(z,x,1) solve(zz)…

%%各时段电量需求 clc close all clear all cost_gd = [2200 1800 3800 4600]; cost_bj = [2.7 2.2 1.8 3.6]; cost_qd = [4000 1500 2500 1000]; fontsize = 24; axes('Parent',figure,'FontSize',fontsize); t=1:length(cost_gd); yyaxis left ; ylim([0,5000]); xlim = ([-2,…

clc; clear all; close all; x1 = 0:.1:40; y1 = 4 * cos(x1)./(x1 + 2); x2 = 1:.2:20; y2 = x2.^2 ./ x2.^3; h11 = line(x1,y1,'Color','r'); ax1 = gca; aa = get(ax1,'position'); aa(1:2) = aa(1:2) + .1; aa(3:end) = aa(3:end) - .2; % aa(1:2) = aa(1:2)+1 set(…

原文链接 多重网格法是一种用于求解方程组的方法,可用于插值.解微分方程等. 从专业角度讲多重网格法实际上是一种多分辨率的算法,由于直接在高分辨率(用于求解的间隔小)上进行求解时对于低频部分收敛较慢,与间隔的平方成反比.就想到先在低分辨率(间隔较大)上进行求解,因为此时,间隔小,数据量小,进行松弛时的时空耗费小,而且收敛快,而且一个很重要的优点是在低分辨率上对初值的敏感度显然要低于对高分辨率的初值的要求.这一点是显而易见的,例如我们平时看一个很复杂的物体,在很远的地方,你可能就觉得它是一个点或一个…

1.老老实实把课本上的题目做完.其实说科大的课本难,我以为这话不完整.科大的教材,就数学系而言还是讲得挺清楚的,难的是后面的习题.事实上做1道难题的收获是做10道简单题所不能比的. 2.每门数学必修课至少要看一本参考书,尽量做一本习题集. 3.数学分析别做吉米,除非你太无聊,推荐北大方企勤的习题集.此外注意一下有套波兰的数学分析习题集,是不是搞得到中文或英文版. 4.线性代数推荐普罗斯库列科夫的<<线性代数习题集>>和法捷耶夫的<<高等代数习题集>>.莫斯科…

本文是笔者在线看Lektorium上John Morgan在圣彼得堡国立大学欧拉研究所的讲座做的笔记.第一讲以如下内容组成 1. 黎曼曲面上的联络 黎曼流形$(M^n,g)$中,$M$为$n$维流形,而$g$为正定的黎曼度量,即$g_{ij}(x^1,x^2,\cdots,x^n)dx^i\otimes dx^j$,而$(g_{ij})$是对称正定的. $\nabla$是联络(我们可以把它看成“方向导数”($\nabla_X$为求$X$方向)),它的定义域与值域为$\nabla:Vect(M)\…

matlab快速入门 1矩阵 生成矩阵 ​ % 直接法 a = [1,2,3;4,5,6;7,8,9]; % 冒号一维矩阵 a = 开始:步长:结束,步长为1可省略 b = 1:1:10; % 1,2,...10 b = 1:10; %与上一个等价 % 函数生成 c = linspace(0,15,5); % linspace(开始,结束,元素个数),等差生成指定元素数的一维矩阵,省略个数则生成100个 % 特殊矩阵 e = eye(4); % eye(维数)单位阵 z = zeros(1,4)…

1.坐标轴的视点(viewpoint):从哪个方向看整个坐标系统,这决定了坐标轴的方向和位置,通过view函数实现视点的设置:view([z y ]):(将坐标系统想象为一座房子,而自己是个会飞的天使,最初你在房子的北,海拔角度为0的位置,这是时你的视角是view[0 0],之后可以根据需要设置z和y的值,分别是绕房子转动的方向,和沿空间轴向上妃的高度)z表示视角沿z轴(正方向)旋转的角度,y表示视角沿y轴(水平面)旋转的角度2.坐标轴显示.关闭.样式axis off; % 去掉坐标轴axis…