chd_lkl的博客

julia> using Plotsjulia> gr()Plots.GRBackend()julia> plot(rand(10),rand(10))julia> savefig("1.png")julia> 

s1 = "The quick brown fox jumps over the lazy dog α,β,γ"# returns the first index of a chari = search(s1, 'b')println(i)#> 11# the second argument is equivalent to the second argument of spli...

# String: 转化和格式化# strings can be converted using float and Inte_str1 = "2.718"e = float(e_str1)println(5*e)#> 13.5914num_15 = parse(Int, "15") # 对于整数的字符串，需要Int与parse配合使用println(3*num_15)#&g...

# strings are defined with double quotess1 = "The quick brown fox jumps over the lazy dog α,β,γ"println(s1)# println: adds a new line to the end of output# print: can be used if u dont want adds ...

# 本文将简单地介绍julia如何其他sub文件中的函数# 1. 先创建一个sub.jl文件。内容如下：function fun1(x::Float64, y::Float64) z = x + yendfunction fun2(x::Float64, y::Float64) z = x * yend# 2. 建立主程序文件main.jl。内容如下：includ...

# 楼主原来一直用fortran语言写代码，最近想学学julia语言。在数组下标这块感觉和fortran有点不一样。特此记录# 在fortran中，比如定义如下数组# real(kind=8) :: a(10) !// 表示有10个元素，默认下标从1开始，即a(1),a(2),...,a(10)# 当然你也可以指定下标的起始位置与终止位置，比如# real(kind=8) :: a(-5...

# 有如下数据1,23,45,67,89,10# 说明，本例代码在julia1.0上测试通过# 直接用数组读取using DelimitedFilesa = readdlm("123.txt", ',') # 数据之间用英文逗号隔开，注意readdlm中是单引号# a = readdlm("123.txt") # 读取用空格分割的数据println(a)for i i...

在开始讲解之前，我们先建立一个简单的线性系统：A = rand(3,3)x = fill(1, (3,1))b = A * x**Factorizations:**1. LU分解PA = LU此处P为置换矩阵，L是下三角矩阵（对角线元素为1），U是上三角矩阵julia使用函数lufact可以获得矩阵A,L,UAlu = lufact(A) # 在julia1.0中，Alu =...

本文将简要介绍julia1.0发布以后，有关文件的函数介绍下面的代码为批量产生10个文本文件，并在每个文本文件中写入其文件名的示例代码using DelimitedFilesn = 10filelist = Vector{String}(undef,n) # 用来存储文件名的字符串数组for i in 1:n filelist[i] = string(i, ".txt") ...

function printsum(a) # summary generates a summary of an object println(summary(a), ": ", repr(a))end# arrays can be initialised directly:a1 = [1, 2, 3]printsum(a1)# or initialised emp...

function printsum(a) # summary generates a summary of an object println(summary(a), ": ", repr(a))end# repeat can be useful to expand a gridm1 = hcat(repeat([1,2], inner=[1], outer=[3*2])...

function printsum(a) println(summary(a), ": ", repr(a))end# dicts can be initialised directly:a1 = Dict(1=>"one", 2=>"two")printsum(a1)#> Dict{Int64,String} with 2 entries: Dict(2=...

千呼万唤始出来，juliaPro1.0.1终于出来了。在使用Plots画图时，与以前有点不一样，现在记录如下。juliaPro1.0.1首先要安装Plots。本人安装步骤如下:1. using Pkg2. Pkg.add("Plots")耐心等待安装完成......3. using GR测试代码如下:using GRx = rand(100);y = rand(100);p...

在处理数据的过程中，我们有时候要对旧文件中的数据进行挑选，写入新的数据文件本博文给出julia0.6下运行的程序n = 10# 打开要写入数据的文件名newfile = open("newfile.txt", "w")for i = 1:n filename = string(i, ".txt") fileio = open(filename, "r") # 打开每一...

using DelimitedFilesn = 10for i = 1:n filename = string(i, ".txt") fileio = open(filename, "w") # 注意这里是双引号 # writedlm(fileio, [1 2 3]) 按行写 # writedlm(fileio, [1, 2, 3]) 按列写 writ...

# function to calculate the volume of a spherefunction sphere_vol( r ) return 4/3*pi*r^3end# 如果函数只有一个表达式，可以用下面简洁的形式quadratic( a, sqr_term, b ) = ( -b + sqr_term ) / 2*a# 计算x 0 = * x ^ 2 + b...

using DataFramesshowln(x) = (show(x); println()) # 这句代码很重要# A DataFrame is an in-memory databasedf = DataFrame(A = [2, 5], B = [e, pi], C = ["x", "y"])showln(df)#> 2x3 DataFrame#> |----...

fname = "test.txt"# using [do] means the file is closed automaticallyopen(fname,"r") do f for line in eachline(f) print(line, "\n") endend#> this is a simple file containing#> text ...

# 这一节讲解type的使用function printsum(a) println(summary(a), ": ", repr(a))endtype Person name::AbstractString male::Bool age::Float64 children::Int64endperson1 = Person("luk", f...

for i in 1:5 print(i, ", ")end#> 1, 2, 3, 4, 5,# =与in作用相同print("\n")for i = 1:5 print(i, ", ")endprint("\n")# array can also be looped over directly:a1 = [1,2,3,4,5]for i in a1...

Multiple dispatch是julia语言的一大特色，也叫多重分派。可以使Julia更通用更快捷。为了理解Multiple dispatch，我们看下面的例子。f(x) = x^2julia语言会自动判断接收参数的类型并对其进行运算f(10) = 100而对于f([1,2,3,4])则会报错，因为在数学里没有定义两个向量相乘，只定义了两个向量的内积**指定输入参数的类型**...

# 本人为了使用Plots简直快吐血了，因为是新手，对整体流程拿捏很不清楚。现在把我的经历写成log，供大家参考# 本人一开始安装了PyPlot，但是用得很不舒服。在网上重新寻找别的package，感觉Plots更适合我的口味。再安装了PyPlot之后add Plots各种报错。# 最终不得以把.julia下的lib文件夹以及cache文件夹删除。重新Pkg.add("Plots").这次安装...

# 假设文本文件test.txt中有如下数据100,200,300 # this is the number# 如何只读取前三个数据呢？注意这里的英文逗号以及300后面的空格都是要有的（阵对本文的方法）str = readuntil("test.txt"," ")list = Vector{String}(3) # 这里的3是文件中所要读取的数据个数list = split(str...

# 本文将简要介绍怎样向文本文件中写入数据# 代码如下x = linspace(1,10,10)y = linspace(11,20,10)f = open( "dat.txt", "a" ) # 由于一开始可能并没有dat.txt,所以这里用"a"for i = 1:length(x) writecsv( f,[x[i],y[i]] )end close( f )#...

# 本代码将展示如何批量产生文本文件# 第一种方法：利用外部程序实现[可参考我写的这一篇博文](https://blog.youkuaiyun.com/chd_lkl/article/details/81430166)n = 10for i = 1:n str = string( i, ".txt" ) run( `touch $str` )end # 第二种方法，利用Julia语...

# 在写代码的过程中，有时候会用到外部程序# julia提供了run函数可以实现, run的函数格式为run( `command line` ): 此处的command line以linux/unix命令行为准# 比如要创建10个文件夹（win10环境）n = 10for i = 1:n str = string( i ) run( `mkdir $str` ) # 此...

# 本博文将简单讲述julia如何从屏幕输入数据（包括字符与数字）str = readline() # 在屏幕上输入数据a,b时，用英文逗号分隔开a, b = split( str, "," )println( "a : ", a )println( "b : ", b )...

# 假设有如下的文本文件dat.txta, 1.0b, 2.0c, 3.0d, 4.0e, 5.0f, 6.0g, 7.0# 目标：只存储第二列数据到数组f = open( "dat.txt", "r" )n = countlines( f )seekstart( f )x = zeros(n)for i = 1:n str1, str2 = split( r...

# 假设有文本文件中存储中英文字符如下chineseenglish我爱周杰伦# -----------------------------# 代码如下f = open( "str.txt", "r" )n = countlines( f )seekstart( f )for i = 1:n str = readline( f ) println( str )...

在julia中优化代码的一种方法是确保类型的稳定性。如果函数中一些变量的类型不确定，编译器无法对这些变量进行确认，执行效率就会降低。当函数内部变量的类型是确认时，执行效率会提高我们来看下面的两个函数func1和func2下面的代码就是典型的函数内部变量类型“模糊”现象：因为输出时的数据类型可能是Float，Int，Rationalfunction func1() s = ran...

# 假设有文本文件dat.txt,里面存储的数据如下1,22,33,44,55,6# **用第一种方式读取，上面数据末尾的空行一定要有**# 下面本文将演示如何读取# 第一种方法f = open( "dat.txt", "r" )n = countlines( f ) # 计算数据行数seekstart( f ) # 返回文件起始位置for i = 1:n ...

# 本文将简要介绍一下如何同时打开多个文本文件并对其写入数据# 代码如下n = 10 # 文件数filelist = Vector{String}(n) # 创建10个空字符串for i = 1:n filelist[i] = string( i, ".txt" ) # 批量产生文件名endx = linspace(1,10,10)y = linspace(11,20,1...

# 本文将简单介绍julia如何使用module，[进一步的学习请参考文档](https://docs.julialang.org/en/stable/manual/modules/)或自己琢磨# 首先建立一个module文件，假设叫module.jl,内容如下module MyModule export my, he, f my = 123 he = 24 function...

# 下面这个代码演示输出到文本的数据不换行输出x = linspace(1,10,10)f = open( "dat.txt", "a" )for i = 1:length(x) print( f, x[i], " " ) # print默认不换行endclose( f )# 下面的代码演示换行输出x = linspace(1,10,10)f = open( &

# 本文将演示如何使用Julia语言实现三次样条插值# 如果没有安装Dierckx包，首先要Pkg.add("Dierckx"),然后Pkg.update()# 当然也可以用Interpolations.jl,但对其语法不太熟悉# 用下面代码做示例x = linspace(0,pi)v = sin.(x)xq = linspace(pi/2,pi)# 计算三次样条插值spl = ...

# 本文将简单描述使用PyPlot包在一个画布上画多张图using PyPlotx = -10:0.001:10subplot(221)plot(x, x)subplot(222)plot(x, x.^2)subplot(223)plot(x, x.^3)subplot(224)plot(x, x.^4)...

# function的声明方法# 第一种function myadd(x, y) x + yend function# myadd(2,3) = 5# 第二种f(x, y) = x * y# f(2, 3) = 6# 第三种f = x -> x / 10 # 这种形式的函数貌似只能是单变量函数# f(1) = 1 / 10 = 0.1# Duck-typ...

# if语句的结构为if *condition 1* *option 1*elseif *condition 2* *option 2*else *option 3*end# 判断一个数的范围a = 51if ( a > 50 ) && ( a <= 100 ) println( "a" )elseif ( a > ...

# while循环：格式如下while *condition* *loop body*end# 例如用while计算1-100的和s = 0i = 1while i <= 100 s += i i += 1endprint( s ) # for循环：结构如下for *var* in *loop iterable* *loop body*...

julia1.0发布之后，有的函数被废弃，有的语法被更改。不过不要抱怨因为1.0发布之后，julia的语法会趋于稳定，而且修改语法也是为以后更好地发展。1. 与文件有关的语法julia1.0发布之后，需要using DelimitedFiles之后，相关的writecsv等函数才能使用2. 线型代数相关的函数首先需要using LinearAlgebra。而且相关的lufract, ...