- 博客(13)
- 收藏
- 关注
RSS订阅原创 使用word中zotero插入参考文献教程
zotero这个插入参考文献的功能可以自动根据引用在文章的位置更新序号，可以自动生成引用的参考文献，可以统一设置引用格式会自动批量修改。
2023-12-14 17:17:52
6870
2
原创 java第一章练习题
5．编写一个Java应用程序，输出全部的希腊字母。提示：找到第一个希腊字母α在Unicode中的编码，最后一个希腊字母ω在Unicode中的编码。然后循环输出这个范围内所有的字符，就得到了希腊字母表public class helloworld { public static void main(String[] args){ for (int i = 945; i <= 969; i ++){ System.out.print((char)i); } }}.
2021-11-13 00:59:24
950
原创 数值计算方法第三章上机实践题三3.1
a=[3 59.14 3 1;5.291 -6.130 -1 2;11.2 9 5 2;1 2 1 1];b=[59.17;46.78;1;2];n=4;a(:,n+1)=b;for i=1:n for k=i+1:n lb=a(k,i)/a(i,i) a(k,:)=a(k,:)-lb*a(i,i) endend
2021-10-16 15:33:46
253
原创 matlab用雅可比（Jacobi）迭代法求解方程组
function x=myjacob(A,b,x0)myeps=0.0001;n=length(b);err0=1;while err0>myeps x1=zeros(n,1); for i=1:n s1=A(i,i:i-1)*x0(1:i-1); s2=A(i,i+1:n)*x0(i+1:n); x1(i)=(b(i)-s1-s2)/A(i,i); end err0=norm(x1-x0);x0=x1;end.
2021-10-16 15:32:51
1710
1
原创 matlab用三角分解法解函数
function x=sanjiao(A,b)he1=0;he2=0;sizes = size(A)for i=1:sizes(1) A(i,1)=A(i,1)/A(1,1)end%将矩阵A改成L和U矩阵%理论上L和U是两个单独的矩阵，为了减少内存，直接在A矩阵上修改，用到的时候直接取A矩阵的值for r=2:sizes(1) for i=r:sizes(2) for k=1:r-1 he1=he1+A(r,k)*A(k,i).
2021-10-16 15:16:51
2430
1
原创 数值计算方法第三章例题3.1代码实现
第一种：三角形方程组的解a=[1 1 1;0 4 -1;2 -2 1];b=[6;5;1];n=3;x=zeros(n,1);for i=n:-1:1% s=0;% for j=i+1:n% s=s+a(i,j)*x(j);% end s=a(i,i+1:n)*x(i+1:n); x(i)=(b(i)-s)/a(i,i);endx第二种：高斯消去法求解A=[1 1 1;0 4 -1;2 -2 1];b=[6;5;1];n...
2021-09-25 14:58:10
149
原创 matlab二分法求根
s1=1;s2=2;answer=0;f=inline('x^3+4*x^2-10','x')%别忘记计算需要迭代的次数n%这个n加一减一是啥意思n=fix((log(s2-s1)-log(eps))/log(2)-1)+1;%fix()函数是取整数for j=1:n if((f(s1+s2)/2)*f(s2)<0) s1=(s1+s2)/2; answer=s1; elseif ((f(s1+s2)/2)*f(s2)>0) ...
2021-09-18 15:04:36
544
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人