- 博客(8)
- 资源 (1)
- 收藏
- 关注
RSS订阅原创 matlab行阶梯变换，LU
lu函数厂商代码 1 function [L,U,x]=lux(A,b) 2 %LU 分解法解线性方程组（列主元LU分解） 3 [n,n]=size(A); 4 p=eye(n);%p记录了选择主元时候所进行的行变换 5 for k=1:n-1 6 [r,m]=max(abs(A(k:n,k))); %选列主元 7 m=m+k-1; 8 if(A(m,k)~=0) 9 if(m~=...
2021-12-31 10:23:13
486
原创 matlab笔记
Function：max、min,column的最值，可以max（max（…））从而获取更精准的sum，column的和mean，整个column的平均sort，每个列从小排到大（具有独立性）sortrows，按每行首元素进行对所有行进行从小到大排序（绑定在一起）size，矩阵大小，如输出（3 3）find,从第一个竖着排序后寻找元素，如a=[1 2 3;4 5 6;7 8 9]中find(a==4)得到2rem,求余数，如rem(10,5)=0disp，打印函数（p.
2021-12-30 20:21:13
544
原创 matlab画图持续更新
示例1.[X,Y] = meshgrid(-3:.125:3);Z = peaks(X,Y);figure(1); %此句可以省略subplot(2,1,1);surf(X,Y,Z);subplot(2,1,2);surfc(X,Y,Z);axis([-3 3 -3 3 -10 5]);示例2.[X,Y] = meshgrid(-8:.5:8);R = sqrt(X.^2 + Y.^2) + eps;Z = sin(R)./R;mesh(X,Y,Z)col...
2021-12-07 21:10:12
677
原创 python三位水仙花数笔记
shuixianhua=[]for i in range(1,10): for j in range(10): for k in range(10): if i*i*i + j*j*j + k*k*k == 100*i + 10*j + k: shuixianhua.append(100*i + 10*j + k)for i in shuixianhua: if i == shuixianhua[-1]: print(i) else: pr...
2021-12-05 23:25:26
79
原创 python些许笔记-动态进度条、等边三角形、数码管
# TextProBarV3.pyimport timescale = 50print("执行开始".center(scale // 2, "-"))start = time.perf_counter()for i in range(scale + 1): a = '*' * i b = '.' * (scale - i) c = (i / scale) * 100 dur = time.perf_counter() - start print("\r{.
2021-12-05 11:24:50
77
原创 matlab-绘图-心形线初探
我的心形线初探t=-10*pi:0.01:10*pi;p=3*(1+cos(t));x=p.*cos(t);y=p.*sin(t);plot(x,y,'-ro')grid ontitle('❤');xlabel('x');ylabel('y');大佬：用MATLAB实现心形线 - 知乎 (zhihu.com)
2021-12-04 23:40:23
343
原创 matlab字符串练习（一）
>> ch='awewqe213AEW'%定义hc=[ch(end:-1:1)]%倒置ch14=[ch(1:4)]%取1~4的字符s=find(ch>'a'&ch<'z');%寻找小写字母ch(s)=ch(s)-32%小写变大写ans=length(s)%共多少小写？
2021-12-02 15:27:51
200
原创 matlab水仙花数
>> a=[100:900];m1=rem(a,10);m2=rem(fix(a/10),10);m3=fix(a/100);f=find(a==m1.*m1.*m1+m2.*m2.*m2+m3.*m3.*m3);k=[a(f)]k = 153 370 371 407
2021-12-02 15:18:33
1755
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人